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I . INTRODUCTION 


The Georgia Institute of Technology Engineering Experiment Station 
(i:i:s) has been funded since 1973 by NASA Marshall. Space Flight Center 
(MSFC) to assist the State of Georgia in utilizing Landsat digital 
analysis. In support of the Georgia state and local agencies and the 
usefulness of Landsat applications, the purpose of this project was to 
facilitate a transfer of technology, especially as related to low-cost 
Landsat data analysis systems. The major tasks accomplished during 
this project are discussed below. 

The initial task involved a comprehensive survey of capabilities 
for digital processing of Landsat data in the Southern United States. 
Included was a review of the software and hardware currently being used 
by the southeastern states. Also, information was obtained on the 
commercial services used by those desiring to employ a Landsat data 
processing service. 

As a second task, EES developed a set of minimum requirements 
(hardware and software) for a stand-alone Earth Resources Data Analysis 
System (ERDAS). The requirements developed in this phase represent a 
system that can be constructed at a minimum cost to the prospective 
user and can also satisfy the user’s basic needs. The hardware for 
this system includes a minicomputer, line printer, disk drive, and 
tape drive while the software includes basic Landsat analysis techni- 
ques — level slicing, maximum likelihood pattern recognition, clustering, 
and software scaling of the data. The primary products from such a sys- 
tem are computer printouts and statistics. 

Since the user groups supporting an ERDAS system may vary greatly 
in size and available resources, a comprehensive plan was developed for 
expansion of the basic ERDAS system. This plan details a set of options, 
including approximate costs, and discusses the advantages of special 
equipment as well as the way that each equipment type complements the 
ERDAS concept. 

One of the prime concerns in setting up an operational system is 
the hardware and softxjare support needed to sustain a satisfactory 
operating schedule. Therefore, as the third task for this project, EES, 
through experience with its own ERDAS system, developed minimum re- 


1 


quirements for continuing support for such a facility ^software and 
hardware). A system such as ERDAS is versatile and is likely to be 
used for purposes other than earth resources analysis. It therefore 
will be in use most of the time and will necessarily require the 
availability of a skilled electrical technician. Software support 
may be obtained "in-house” or by contract. Some equipment manu- 
facturers hav i extensive service contracts which prove cost effective. 

As the fourth task of the project, EES has assembled a software 
library for the digital processing of Landsat and other multispectral 
scanner data using, whenever possible, existing techniques from NASA, 
the University System, and other public agencies. Software techniques 
for maximum likelihood classification, linear classification, clus- 
tering, level slicing, registration and rectification, and table look- 
up classification are available for implementation into an ERDAS system. 
Each module is designed on a user keyword structure so that no detailed 
knowledge of programming is needed. Documentation is available on pro- 
grams written or changed exclusively at EES with references provided 
for other programs. Even though the basic software was developed on 
a Data General minicomputer at Georgia Tech, the software library was 
designed entirely in Fortran IV and the routines are transferable to 
other 16 bit minicomputers. 

The final task of the project was the implementation of a limited 
version of the IMGRID Geographic Analysis program on the ERDAS system. 

By providing this system to state and local users, the wide range of 
applications for rectified Landsat data becomes far more evident than 
if the system were geared to produce only land cover maps. IMGRID 
produces a dynamic modeling tool for site planning, erosion control, 
environmental impact, and many other uses. 

Report Organization 

The remainder of this report discusses each of the five tasks in 
detail. Section II discusses the results of the comprehensive survey 
of capabilities for digital rrocessing of Landsat data in the South- 
eastern United States. Section III presents various options for low 
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cost earth resources processing systems and the minimum requirements 
for support of such facilities. Section IV discusses the IMGRID 
geographic analysis program. Appendix A contains descriptions of 
some of the major software available at EES, Appendix B contains 
a listing of the available Fortran software, and Appendix C presents 
flowcharts of some of the main analysis programs. 


II. CAPABILITIES FOR AND USE OF 
DIGITAL LANDSAT DATA IN THE SOUTHEAST 


Numerous individuals and/or organizations who are involved in 
the use of remote sensing data in the southeastern states were con- 
tacted in order to obtain information on their states 1 capabilities 
for digital processing of Landsat data. Included is a review of the 
software and hardware which are currently being used by these states 
as well as information on the commercial services used by those who 
wish to employ a Landsat data processing service. 

The states involved are: Alabama, Georgia, Kentucky, Missouri, 

North Carolina, South Carolina and Tennessee. A brief summary of 
each state's capabilities also includes information on state funded 
projects utilizing Landsat data digitally processed by companies or 
organizations outside the state. 

Alabama 

NASA/Marshall Spare Flight Center (MSFC) has the only current 
capabilities in Alabama for digital processing of Landsat data. 

Their hardware system consists of an IBM 360/75 with two megabytes of 
memory and several discs and tapes and a PDP11/45 with two 250-mega- 
byte discs, three tape drives, and various displays and terminals. 
Approximately 375 computer software routines are available including 
all algorithms developed by any of the NASA centers. However, very 
little analysis of Landsat tapes is actually done with this system at 
present. 

Georgia 

In order to assist the State of Georgia agencies with their desire 
to incorporate digital Landsat data into their planning activities on 
an operational basis, the Georgia Tech Engineering Experiment Station 
(EES) approved the design and acquisition of the Earth Resources Data 


4 





Analysis System (ERDAS) . EES was responsible for selecting and inte- 
grating both the hardware and software components of the system. 

Hardware . ERDAS was designed and constructed to allow true inter- 
active digital processing of all types of remote sensing data. The 
system consists of a set of four modules: (1) minicomputer subsystem, 

(2) input medium, (3) hardcopy output medium, and (4) display sub- 
system. 

The minicomputer subsystem consists of a N0VA-2/10 minicomputer 
with 64K bytes of core memory and a dual Diablo disk system with 5.0 
megabytes of storage for programs or data. 

The input medium for the ERDAS system is a set of two nine track 
dual density (phase encoded/NRZI selectable) magnetic tape drives and 
controller — both drives with a capacity for 10-1/2 inch reel tape. 

One hardcopy output device is a twenty inch electrostatic dot 

matrix printer/plotter. Scaled maps of Earth Resources data can be 

(R) 

made using this device. A CROMALIN' photographic process may then 
be used to generate a color coded outptit hardcopy product. Another 
output method currently in use consists of storing images on a magnetic 
tape and sending these tapes to be made into images by the use of a dig- 
ital film writer. This method is currently very inexpensive. 

The display subsystem consists of a high quality color video image 
analysis system that is interfaced to the minicomputer for complete user 

interaction in the choice of training samples for earth resources classi- 
fication. 

Software . Initially, EES implemented a basic Landsat digital 
analysis program called ASTEP (Algorithm Simulation Test and Evalua- 
tion Program) which was obtained from NASA/JSC and has been extensively 
modified by EES personnel. The ASTEP system was designed as a modular 
program whereby various classification algorithms may be tested against 
one another using a standard input/output system. 
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Software for the rectification of Landsat data to map coordinates 
using a least squares fit of Landsat data to control points was ob- 
tained from NASA/MSFC and transferred to the Georgia Tech computer. 

Other software, such as various spatial clustering algorithms, was 
studied but has not yet been transferred. Software for a table lookup 
formulation for Landsat classification (ELLTAB) was obtained from 
NASA/ERL along with software for rectification, destriping, and polygon 
location of Landsat data. A fast combination table lookup and maximum 
likelihood classifier from NASA/ERL has been implemented which signi- 
ficantly decreases the length of time needed for a scene classifica- 
tion. A fast clustering algorithm from NASA/ERL has also been imple- 
mented. In addition, many software algorithms for image manipulation, 
spatial filtering, rectification, training field selection, and high 
speed classification have been developed at EES. EES software exists 
for image analysis on UNIVAC and CDC large computers and Data General 
minicomputers . 

Projects . EES has been deeply involved in the formulation, plan- 
ning, and implementation of a Georgia Natural Resources Inventory since 
its conception in 1972. During 1972 and 1973, EES presented various 
State of Georgia agencies with the background information needed to 
make an initial assessment of the usefulness of digital Landsat infor- 
mation. A trial project was initiated between the Georgia Department 
of Natural Resources (DNR) and EES in 1973 to test the capability for 
using digitally processed Landsat data to determine land use in the 
Atlanta area. 

EES has been funded since 1973 by NASA/MSFC to assist the State 
of Georgia in utilizing Landsat digital analysis for various resource 
problems within the state. In this multi-year effort, several related 
tasks have been performed in conjunction with numerous local and state 
agencies within Georgia, including the Department of Natural Resources, 
the -Department of Transportation (DOT), and the Office of Planning and 

u 

Budget (OPB) . 
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In 1975 Georgia Tech EES was designated as the technical interface 
with NASA/EEL for the transfer of NASA software to the State of Georgia. 
Under this technology transfer project EES purchased the necessary digi- 
tal processing equipment required for operational processing of Landsat 
data. EES then coordinated a project with the Georgia DNR office of 
Planning and Research for classifying and mapping land cover for die 
State of Georgia. Products of the effort assured further acceptance 
of digital processing of Landsat data as an operational tool for en- 
vironmental analysis. 

Probably the best indicator of the success of the technology trans- 
fer efforts of this project is the number of agencies which have commit- 
ted funds and/or personnel time to a project to map the entire State of 
Georgia using Landsat data. This project is concerned with mapping land 
cover using Landsat data processed with ERDAS and, where appropriate, 
inferring land use. 

The agencies which have committed funds to the mapping project 
include : 

Georgia Department of Natural Resources 
Environmental Protection Division 
Game and Fish Division 
Office of Planning and Research 

Georgia Forestry Commission 

Georgia Office of Planning and Budget 
Bureau of Community Affairs 

U.S, Department of Agriculture 
Soil Conservation Service 
Forest Service 

United States Army Corps of Engineers 
Fort Benning 

Savannah Engineer District 

North Georgia Area Planning and Development Commission 

Coosa Valley Area Planning and Development Commission 

Other organizations which are interested but as yet have supplied 
no funds include: 

Georgia Department of Natural Resources 
Earth and Water Division 

Georgia Department of Transportation 

Five other area planning and development commissions. 
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EES is currently assisting the Georgia DNR in a geographic data 
base demonstration project. DNR manually obtained 30 different data 
variables on a 10 acre cell basis for North Fulton County, Georgia. 

EES is assisting with transfer of the data to the ERDAS system and 
analysis of the data using IMGRID and NIMGRID. 

In addition to the activities at Georgia Tech EES, the Atlanta 
Regional Commission (ARC) was involved in a project to test the 
feasibility of using automatic processing of Landsat data to detect land 
use changes in the seven county planning area of the ARC, and thereby 
update the USGS/LUDA land use map of the Atlanta Region. 

The technique chosen to accomplish the land use change detection 
was the ratioing of two different dates of Landsat data. This was 
accomplished at the EROS Data Center utilizing the Image 100 System. 

The results of the project indicated that ratioing Landsat data was a 
feasible technique for ARC to use in updating the USGS land use maps 
of the area. Accuracy evaluations showed that 91% of the change theme 
was accurate to within about 3 pixels (accuracy sufficient for ARC’s 
purposes). Additional manual analysis was required to identify 
accurately the types of land use changes. 

Ke ntucky 

Kentucky is currently developing the capability for in-house 
processing of Landsat digital tapes. A few state agencies and 
universities have small interactive graphics systems which are capable 
of digital processing, but these are being discarded in favor of a 
central state computer system. The central system includes an IBM 
370/168 MP with 9 megabytes of main storage and various discs. The 
Bureau of Computer Services has been established as the central 
hardware agency. Their personnel have recently visited NASA/ERL to 
obtain compatible software for the system. 
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A representative from Eastern Kentucky University recently 
attended a workshop at GSFC on a joint Appalachian Regional Commission- 
NASA project to apply Landsat to the study of geological lineaments 
using GE's Image 100 System. An ORSA package was ordered from GSFC 
which provides an offline printout of similar pixels vs. significant 
changes, but this is not adequate for final classifications. 

Five professors from Murray State University attended a short 
course at ERL and had plans to return the end of October for hands-on 
experience with digital tapes for strip mining studies in Western Kentucky. 

Previous projects in Kentucky include a survey in 1973 by Earthsat 
of water impoundments larger than two acres. The results were not 
beneficial because the imagery was taken after a flood and strip mines 
appeared as impoundments. However, this led to a Landsat-2 follow-Dn 
study in 1976 for which ERIM provided the digital processing. The 
objective of the project was to look at the operations of inspection and 
enforcement to detect significant violations of surface mining laws. 

Several projects have been contracted with Bendix Corporation for 
Area Development (AD) Districts in conjunction with EPA 208 land use 
planning programs. Color-coded maps and overlays were provided in 1976 
for the Kentuckiana AD (seven Kentucky counties, two Indiana counties) 
at 1:48,000 scale. Area tabulations were also provided for each county. 
Processing for the Big Sandy AD (Prestonsburg, Kentucky) is in progress 
and maps are scheduled for delivery this fall. A contract has also been 
signed for the Green River AD (Owensboro, Kentucky) to be delivered in 
December, 1977. This project will provide information to meet HUD 701 
requirements for land use planning. 

Two counties of Kentucky were included in a digital processing 
project by Bendix for the Ohio-Kentucky-Indiana Council of Government 
(OKI COG). Color-coded maps at 1:62,500 scale and computer tabulations 
for each of 229 drainage basins were provided. Also, five counties 
centered around Fayette County (Lexington, Kentucky) known as the 
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Central Blue Grass Region (not an AD district) were mapped for the 
Army Corps of Engineers. 

The LARS program at Purdue was used to produce a land cover map 
of Henderson County in Western Kentucky (part of Southeastern Indiana 
COG). 

Missouri 

According to the final report of a project on Earth Observation 
Data Mangement Systems in December, 1976, "few agencies (in the five 
state Midwestern region of Illinois, Iowa, Minnesota, Missouri, and 
Wisconsin) now have the staff or computer capabilities to handle 
digital satellite data."^ At the present time, Missouri still has no 
capabilities for digital processing of Landsat data. However, the 
University of Missouri at Rolla recently received a grant from the 
National Science Foundation and has ordered a Comtal Interactive 
Image Analysis System which will be used for such processing. Delivery 
occurred in early 1978. Software will be obtained from NASA/ERL. 

Several Missouri agencies have funded projects for digital process- 
ing in the past. A demonstration study of land use in the Ozarks 

2 

Planning Region of Southern Missouri was conducted by NASA/ERL in 1975. 

The Soil Conservation Service in Missouri has digital tapes processed by 
LARS/Purdue for water analysis and study of soil patterns. Also, the Missouri 
Geological Survey had the University of Kansas in Lawrence process some digital 
tapes for an area around Kansas City, Missouri. 

North Carolina 

Currently North Carolina has no facilities for processing digital 
Landsat data although considerable interest exists for establishing a 
centralized state system. Several projects have been completed by 
Bendix Corporation for EPA 208 planning regions in North Carolina. 

Color-coded land use overlays at 1:96,000 scale were produced using the 
Multispectral Data Analysis System (M-DAS) at Bendix in 1975 for 


Eastwood, et al., "Project on Earth Observation Data Management Systems," 
Final Report, Washington University, St. Louis, Missouri. Prepared for 
GSFC, December 31, 1976. 

2 

"A Computer Implemented Land Use Classification Technique Applied With ERTS 
Digital Data Acquired Over Southern Missouri," Report number 143, April 1975, 
ERL/JSV, A. T. Joyce and J, D. Derbonne*. 
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Planning Region J (also known as Triangle J — five counties in central 
North Carolina — Raleigh, Durham, Chapel Hill areas) . An additional 
product of the analysis process was the generation of statistical data 
by 50x50 meter grid cells in data sets corresponding to 54 7-1/2 
minute USGS quadrangle maps. Overall classification accuracy of the 
land use categories was judged to be around 90% and the cost of 
processing the Landsat scene and generating the products was approxi- 
mately $4.00 per square mile. 

A similar analysis was done in 1976 by Bendix for Planning Region D 
in the Northwestern part of the state. Also, analysis of the Dan River 
Sub-Basin (Roanoke River Basin) was completed this year for the Corps 
of Engineers. NASA had some involvement in the most recent project. 

South Carolina 

South Carolina has no capabilities at present for digital process- 
ing of Landsat data. A proposal has been submitted to NASA/ERL for 
projects which may involve digital processing. Previous activities 
include a contract with General Electric Company for Image 100 process- 
ing of Landsat data to produce color-coded land use maps and area calcu- 
lations for three Council of Government (COG) regions comprising approxi- 
mately 25% of the state. Landsat tapes were processed by NASA to produce a 
map of the Congaree Swamp area for the Wildlife and Marine Resources 
Department. The Land Reserve Conservation Commission and the Bureau of 
Mines visited the EROS Data Center for processing of digital tapes for 
a study of mining areas in South Carolina. 

Tennessee 

The only current capabilities in Tennessee for digital processing 
of Landsat tapes are at the University of Tennessee, Knoxville. The 
facilities are presently being used for image processing projects other 
than earth resources , although previously they have been used for Landsat 
analysis. Oak Ridge National Laboratory has accomplished strip mining 
surveys using some Landsat data. 
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Bendix Corporation completed a digital processing project on water 
resources and strip mining in the New River Drainage Basin in North- 
central Tennessee for the Soil Conservation Service in early 1976. They 
provided a color-coded map of the Basin at 1:62,500 scale, fifteen 
color-coded 7-1/2 minute quads, computer tabulations, and rescanned- 
resampled tapes for the area. 

Summary 

Of all the Southeastern states, Georgia is the most advanced in 
the use of digital Landsat data. The University of Tennessee at 
Knoxville has the capability for digital processing of Landsat data but 
no such projects are currently underway. All other states have relied 
on NASA or commercial facilities. The most extensive use of digital 
Landsat data among these states has been in fulfillment of EPA 208 and 
HUD 701 planning requirements. 

The results of this survey indicate that users of remote sensing 
data in the Southeastern U.S. are increasingly turning to digital 
processing techniques. All the states surveyed have had some involve- 
ment in projects using digitally processed data. Even those states 
which do not yet have in-house capabilities for digital processing are 
extremely interested in and are planning to develop such capabilities. 
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III. DESIGN OF LOW COST 


EARTH RESOURCES DATA PROCESSING SYSTEMS 

In the design of an earth resources data processing system there 
are many factors to be considered. In some cases potential users 
should buy turn-key systems that are currently on the market; in other 
cases, users might consider the design and implementation of their own 
systems by buying components and assembling the systems using their 
own technical expertise. This section deals with alternative systems 
that might be considered by users having the technical expertise for 
assembly of such a system within their agencies or support groups. 

With the appropriate technical personnel available, a significant 
cost savings often may be realized by user design and implementation 
of systems. At least one part- or full-time computer hardware technician, 
one applications software analyst, and several applications programmers 
are desirable for all phases of system design. These requirements may be 
lessened, however, if sufficient support is available from the various 
equipment manufacturers. 

Computer software for digital processing of earth resources and 
other geographically based data is currently becoming available for 
minicomputers at minimal or no cost. Since a significant amount of this 
computer software has been developed under government contracts, it is 
in the public domain and readily available. Thus the costs involved in 
the acquisition of the systems listed here are primaril> the actual 
hardware costs. If they are needed, personnel training costs and sys- 
tems interface costs are extra 

Figure 2 indicates five alternative configurations for low cost 
earth resources data processing systems. At the low end is a nominal 
system consisting of a minicomputer, floppy disk, magnetic tape unit, 
color terminal, and line printer. The estimated cost range for this 
system (depending on the exact components selected) is $22,000 to $45,000. 
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At the upper end of the range (still, however, at a price signi- 
ficantly less than many systems on the market) are systems costing an 
estimated $165,000 to $220,000. This configuration, with substantially 
increased capability over the lower cost system, consists of a mini- 
computer with an array processor, a 96 megabyte disk, dual magnetic 
tape drives, a digitizer, a color display, a line printer, and a film 
recorder. 

In general, as the systems increase in cost, the processing sophis- 
tication is improved and the speed with which a data set can be analyzed 
increases rapidly. Thus, for users requiring only a low volume of pro- 
cessed data, a system in the low or middle cost range might be suitable. 
For users desiring a faster processing speed and an increased through 
put, a more expensive system might be in order. 

The breakdown of the total costs for each system is shown in Table I. 
For each system, the high and the low cost estimate for each component 
is given. Thus, by selecting particular components with a greater cr 
lesser capability, systems could be configured that cost anywhere within 
the range of the cost extremes given. 

To facilitate the estimation of costs for system configurations 
other than those listed in Figure 2, the costs of the individual 
components are given in Table II. Using these data, the approximate 
costs of many more low cost system configurations could be derived. 

Typical vendors of such components are given in Table III. 

Estimated Cost for ERDAS Support 

A minimum cost estimate for the support of the Georgia Tech ERDAS 
system over a one year period totaled $2,500. This estimate includes 
approximately 400 hours of a resident part-time technician along with 
all electrostatic line printer supplies and any additional maintenance 
charges incurred when repairs exceeded in-house capabilities. 

This estimate Is approximately 3.3% of the total cost of the ERDAS 
system. As maintenance contracts usually run about 10% of system costs 
per year, a significant savings was realized using in-house methods of 
maintaining the equipment. 
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K 

L7.NE PRINTER 

$ 

5 

- 

8 

K 

COLOR DISPLAY 

$ 

19 

- 

23 

K 

(256 X 256) 







T 

53 

- 

74 

"k 

TAPE DRIVE 

$ 

12 

— 

14 

K 

CPU 

$ 

8 

- 

IS 

K 

DISK 

$ 

9 

- 

11 

K 

ELECTROSTATIC 






PRINTER/PLOTTER 

$ 

10 

- 

20 

K 

COLOR DISPLAY 

$ 

30 

- 

35 

K 

(3 images 256 X 256) 







T 

69 

- 

98 

”k 

TAPE DRIVE 

$ 

12 

_ 

14 

K 

CPU 

$ 

8 

- 

18 

K 

DISK 

$ 

20 

- 

35 

K 

FILM WRITER 

$ 

40 

_ 

50 

K 

ARRAY PROCESSOR 

$ 

40 

- 

45 

K 

LINE PRINTER 

$ 

5 

- 

8 

K 

DIGITIZER 

$ 

5 

- 

15 

K 

COLOR DISPLAY 

$ 

30 

- 

35 

K 


(3 images 256 X 256) 

$165 -.220 K 
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TABLE II. INDIVIDUAL COMPONENT COST ESTIMATE 


1 

I 


1) 

1 Mag Tape Drive & Controller 

$ 5 

- 10 

K 

2) 

2 Mag Tape Drives & Controller 

12 

- 14 

K 

3) 

16 Bit CPU with 32 K Memory- 

8 

- 18 

K 

4) 

Line Printer 

5 

- 8 

K 

5) 

Floppy Disk 

2 

- 4 

K 

6) 

Dual 2.5 Megabyte Drive & Controller 

9 

- 11 

K 

7) 

10.0 Megabyte Cartridge Disk 

9 

~ 11 

K 

8) 

Electrostatic Printer/Plotter 

10 

- 20 

K 

9) 

Array Processor 

40 

- 45 

K 

10) 

96 Megabyte Disk 

25 

- 35 

K 

IX) 

Color Terminal 

2 

- 5 

K 

12) 

Color Display Image 256 x 256 

19 

- 23 

K 

13) 

Color Display 3 Image 256 x 256 

28 

- 35 

K 

14) 

Film Writer 

40 

- 50 

K 
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TABLE III. TYPICAL VENDORS 


COMPUTERS 


Data General Corporation 
Digital Equipment Corporation 
Hewlett Packard Corporation 

MAGNETIC TAPE DRIVES 

Digi-Data Corporation 
Pertec Computer Corporation 
Kennedy, C. J. Company 


DISKS 


Control Data Corporation 
Data General Corporation 
Digital Equipment Corporation 

COLOR DISPLAY SYSTEMS 

Aydin Corporation 
Comtal Systems Corporation 
ITT Grinnell Corporation 
Ramtek 

ELECTROSTATIC TRINTERS/PLOTTERS 

Versatec Incorporated 
Varian Data Machines 
Gould Incorporated 

DIGITAL FILM RECORDERS 

Dicomed Corporation 

Optronics International Incorporated 


DIGITIZERS 


Summagraphics 
Bendix Corporation 
Talos System Incorporated 
Aristo Graphics Corporation 

LINE PRINTERS 

General Electric Company 
Okidata Corporation 
IBM 

Varian Data Machines 
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IV. THE IMG RID GEOGRAPHIC ANALYSIS PROGRAM 


IMGRID, an information manipulation system for grid cell data 
structures, is a package of computer programs designed for the analy- 
sis of natural resource and land planning data which is qualitative 
in its thematic content and varies over geographic space. Originally 
developed by David Sinton of the Harvard University Department of 
Landscape Architecture, the IMGRID system has been adapted for use on 
mini-computers at Georgia Tech. 

The IMGRID system has been designed for people having no pre- 
vious experience with computers. The basic operations are controlled 
with simple keyword commands which may be used with a basic knowledge 
of planning principles but without any knowledge of programming. Thus, 
a link is provided for easy access and manipulation of digital data 
bases. Some typical project applications of the IMGRID system are: 

River basin planning 

Siting of facilities such as airports 
or sanitary landfills 

Environmental impact statement review 

Visual Analysis 

Project review by regional planning agencies 
KEYWORD STRUCTURE 

The basic structure for an IMGRID keyword command involves three 
processes : 

1. Retrieve one or more data elements from the da-a file. 

2. Transform or manipulate the values for each grid cell 

in the data elements retrieved. 

3. Store the new data element created in the data file. 

On the Georgia Tech ERDAS, keyword commands are entered on a CRT 
terminal, executed by the NOVA II mini-computer, and output either in 
color on the Comtal video display, or in black & white on a dot matrix 
printer. The data files are usually stored on tape then read into a 
disk file. The IMGRID program itself also resides on a disk. 
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The keyword commands fall basically into 9 groups: 

1. ) Data Entry & Management 

Two keywords can be used to enter new information into 

the data base: 

• STORE - operates on an element (data variable) by 
basis. When a new variable, such as a soil type or 
slope category is added to the data base for all cells, 
STORE can be used. 

• UPDATE - operates on a cell by cell basis. If a single 
cell changes characteristics for a particular variable, 
such as a land use change from agricultural to residen- 
tial, UPDATE is used along with the row and column loca- 
tion, plus the new value of the cell. 

Three keywords which are used for data management are: 

■ RELOC - permits a data element to be moved to a new loca- 
tion within the data file, such as relocating the results 
of an analysis as a new data element in another location. 

• RENAME - allows the name of a data element to be changed 
without affecting the contents of the data. 

• LIST - allows a user to list the names associated with 
the contents of part or all of the data base. 

2 . ) Delimiter Keywords 

• MODEL - used as a first keyword in a sequence of keywords 
defining an analysis. The primary function of MODEL is 

to assign a title to subsequent analyses. It also includes 
the function of the CLEAR keyword. 

• CLEAR - clears the results of operations performed by pre- 
vious keywords, 

• END - identifies the completion of IMGRID input. 
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3. ) Display Keywords 


* SYMBOL - allows a user to enter a specific set of charac- 
ter symbols to be used when making a map display. 

• MAP - makes a graphic display of the contents of a data 
element in the file. This can be a data variable such as 
a MAP of a slope, or of the results of an analysis such as 
vulnerability to soil erosion. 

• TEXT - premits user to insert textural descriptions of 
the procedures being undertaken. 

4. ) Spatial Analyses 

* SEARCH - generates a set of values which identify the proxi- 
mity of each cell in the study area to a specified condition, 
such as roads, rivers, or airports. The analysis determines 
how far every other location is from the preselected data 
items . 

► ASEARCH - operates on a cell by cell basis examining for a 
prespecified radius, the conditions around every cell, or a 
defined subset of cells in the study area, such as how many 
cells of wetlands are there within a five cell radius of a 
landfill. 

5. ) Rescaling or Restructuring of Data Values 

• RECODE - assigns new values to an old set of values for a 
data element, such as recoding an old set of values for 
element "land market value" to reflect land use or tax policy 
changes. Recode assigns values in the range of 0-9. 

* X.REC0DE - an extended recode which assigns values in the 
range of 0-19. The RECODE keywords assign the last value 
specified in the case where multiple elements are being 
combined with different value scales. 
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' OVERLAY - assigns the highest value, rather than the 
last value, when multiple elements are being combined. 

The system will over-ride the values assigned on the 
first element if the new values generated by the rescaling 
are greater than the previously existing values. 

6 * ) Logical Combinations of Elements 

’ MATRIX - results in a series of values which identify 
each of the possible combinations of features, such as 
the highly erodible soils on steep slopes. 

7 . ) Reject Conditions 

* REJECT - permits the user to identify a group of cells 
which must be ignored (’’dropped out") in all analysis 
and display of the data, such as rejecting unstable soils 
in an analysis for siting a large industrial complex. 

8 . ) Mathematical Manipulation of the Data 

* MULTPLY - allows user to multiply the values in one data 
element by the values in a second data element. 

* INDEX - generates a weighted index of several data elements 
and also provides for addition and subtraction. 

* NORMAL - takes data values over large ranges and normalizes 
them to the range of 0-99. 

’ REDUCE - generates values in the range of 0-19. 

9. ) User Written Fortran Subroutines 

* USERSUB - allows user to write a Fortran subroutine to 
perform any set of computations that is desired and call 
that subroutine through the USERSUB keyword. 

Figure 3 shows an example of IMGRID output where elements of 
slope, depth to water table, soils, vegetation, proximity to roads, 
and travel time were weighted, overlaid, recoded and mapped to pro- 
duce an attractiveness model for industrial parks. The darkest cells 
are best for industrial sites usii.g the criteria specified. 
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Because of the large size of IMGRID and the limited address space 
of a 16 bit computer, a method was devised to partition the IMGRID pro- 
gram into independent subroutines. Each subroutine corresponds to a 
keyword in the above discussion. By making each major subroutine an 
overlay, IMGRID was implemented on the NOVA 2 minicomputer. A core 
part of the IMGRID program resides in the computer at all times, and 
when a particular keyword is selected, only programs that provide that 
function are loaded into the main memory of the minicomputer from disk 
storage. While some delay occurs in the program execution due to the 
overlay procedure, the time required is negligible compared to the 
execution time of the analysis modules. 

The structure of the minicomputer IMGRID program allows two data 
variables in a format of up to sixty elements by sixty elements to be 
analyzed at one time. In the Atlanta, Georgia area, this array size 
is sufficient to represent data for a USGS seven and one half minute 
quadrangle in approximately ten acre elements. For a large area, 
therefore, data for a data set of the whole area would have to be par- 
titioned into segments dependent on the selected cell size. 

All keyword functions with the exception of the search algorithms 
could easily be performed on such a segmented data set. Any application 
for which a search from some criterium is specified, however, would en- 
counter problems when needing to search beyond the boundaries of one 
data segment. In addition, the normal mode of operation for IMGRID is 
execution in a batch runs tr earn in which a set of input data is required 
in special formats. While this method is sometimes desirable when 
using cards as input to a large computer, an alternate interaction 
method should be considered when dealing with minicomputers where the 
user has direct contact with the computer itself and he may often be 
the only user using the system at one time. 

In an effort to alleviate some of the difficulties given above a 
new minicomputer version of IMGRID, called NIMGRID, has been developed 
at EES. NIMGRID reflects to a large extent the program philosophy of 
the original IMGRID program regarding types of functions performed on 
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the data base, but a new philosophy of data storage and retrieval and 
program interaction was implemented in early 1978. The data storage 
philosophy envisions each data variable as being represented by a 
raster data set which is only constrained by the disk space available 
to the user of a particular system. The data are processed on a line 
by line basis, and even for the search algorithm, only two lines of 
the data set are needed in the computer at one time. While this 
method involves many more input/output operations than the IMGRID 
method, very little degradation in preformance has been noted in com- 
paring N IMGRID and IMGRID. 

The prime advantage of this storage and retrieval philosophy is 
that the size of the data set that may be processed by a user is vir- 
tually unlimited. Also, searches may be made over large areas without 
problems associated with crossing data segment boundaries. Another 
virtue of the NIMGRID system is that the program has been made inter- 
active. Prompting of the user occurs, giving him all possible choices 
and asking him to select the desired command. No knowledge of FORTRAN 
and very little knowledge about the particular minicomputer is re- 
quired of the user. Although NIMGRID is still under development, 
applications personnel from the Georgia Department of Natural Re- 
sources are currently using the program at EES with ease. 
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APPENDIX A 


DESCRIPTION OF MAJOR SOFTWARE 
RELATED TO IMAGE PROCESSING AND 
EARTH RESOURCES AT EES 


b 


PROGRAM NAME: RECRAW 

LANGUAGE FORTRAN IV 

COMPUTER SEL 32/55, CDC Cyber 74 

SPECIAL PERIPHERALS: 2 Tape Drives 

Disk Files 


PURPOSE: RECRAW uses a first order transformation matrix from COORD 

to resample raw Landsat data by a selectable nearest neigh- 
bor or bilinear interpolation to format the data into a 
standard coordinate system. 


PROGRAM NAME : RECCL 

LANGUAGE: FORTRAN IV 

COMPUTER: Data General NOVA 2, CDC Cyber 74, SEL 32/55 

SPECIAL PERIPHERALS: 2 Tape Drives 

Disk Files 


PURPOSE: RECCL uses a first order transformation matrix computed by 

COORD to resample Landsat classified data using nearest 
neighbor and format the data into a Latitude-Longitude or 
UTM coordinate svstem. 


PROGRAM NAME: 3DPL0T 

LANGUAGE: FORTRAN IV 

COMPUTER: Data General NOVA 2, CDC Cyber 74 

SPECIAL PERIPHERALS: Carcomp plotter or printer/plotter 

Disk Files 


PURPOSE: 3DPL0T provides a perspective view of three dimensional 

data (r = f(x,y)). Viewing angle and scaling may be spe- 
cified by the user. This technique is especially useful 
in analysis of topographic related data. 
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PROGRAM NAME: 


MAXIS (NASA/ERL) 


LANGUAGE: FORTRAN IV 

COMPUTER: SEL 32/55, CDC Cyber 74 

SPECIAL PERIPHERALS: Two Tape Drives 

Disk Files 


PURPOSE: MAX18 was developed by Ronnie Pearson of NASA/ERL (Slidell, 

Louisiana) as a fast classifier of Landsat MSS data. The 
program is a combined table look-up-maximum likelihood type 
of classifier which uses the best points of each technique. 
Instead of creating a look-up table defining the boundaries 
of statistical distribution of signatures, this technique 
iterates quickly through the data building a table of where 
in channel space the majority of the data in a Landsat scene 
lie. Then, these vectors are classified using a maximum like- 
lihood scheme. A second iteration is made through the data 
for classification. Each pixel of MSS data is checked to see 
if the data vector associated with that pixel is in the already 
classified data table. If so, the classification is derived by 
simply indexing into the classified table. If the data vector 
has not already been classified a maximum likelihood decision 
rule is used. There is a distinct trade off between amount of 
storage and speed of classification. On the SEL machine a 
Landsat scene may be classified into 60 classes in 1.1 - 2 hours. 


PROGRAM NAME : 2DFFT 

LANGUAGE: FORTRAN IV 

COMPUTER: Data General NOVA 2, CDC Cyber 74 

SPECIAL PERIPHERALS: Magnetic lape Drives 

Disk Files 


PURPOSE: 2DFFT performs a two dimensional Fast Fourier Transform on 

image data. Options in the program include: 

a) Forward transform 

b) Inverse transform 

c) Image filtering by specification of 
one of many optional filters. 
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PROGRAM NAME : 


FILTER 


LANGUAGE: FORTRAN IV 

COMPUTER: Data General NOVA 2 

SPECIAL PERIPHERALS: Image display system and 1 tape drive 

or 2 magnetic tape drives 


PURPOSE: FILTER is a program whose concept was taken from a technique 

used at USGS, Flagstaff, Az. for high frequency enhancement 
of image data. The basic method involves taking a windowing 
approach to create a low pass filtered image and subtracting 
that image from the original image. The resultant is essen- 
tially a high pass filtered image. By adding the high pass 
filtered image to the original image the high frequency en- 
hancement is achieved. 


PROGRAM NAME: IMGRID (Harvard U.) 

LANGUAGE : FORTRAN TV 

COMPUTER: CDC Cyber 74, NOVA 2 (DaLa General) 

SPECIAL PERIPHERALS: Disk Files 

PURPOSE: IMGRID is a general purpose geographic data base manipulation 

pro cram developed by the Harvard University Graduate School of 
Landscape Design. This program provides for manipulation of 
multiple data variables related to the same geographic area 
and ir. a sridded format. Spatial searching, statistics genera- 
tion, and modeling via multivariable weighting parameters are 
key cc its analysis capability. Visual and environmental impact 
analysis are two uses of such a system. 

PROGRAM NAME: TOPO 

LANGUAGE : FORTRAN IV 

COMPUTER: Data General NOVA 2 

SPECIAL PERIPHERALS: lape Drive 

Image Display 


After breaking down the NCIG topographic tapes into NOVA 4096 
record blocks with the CDC Cyber 74, TOPO unpacks and displays 
the data on the imaging system. Various scales may be represented 
by selection of parameters in the program. 


PURPOSE: 


PROGRAM NAME: FILTER 

LANGUAGE: FORTRAN IV 

COMPUTER: Data General NOVA 2 

SPECIAL PERIPHERALS: Image display system and 1 tape drive 

or 2 magnetic tape drives 

PURPOSE: FILTER is a program whose concept was taken from a technique 

used at USGS, Flagstaff, Az. for high frequency enhancement 
of image data. The basic method involves caking a windowing 
approach to create a low pass filtered image and sub tree ting 
that image from the original image. The resultant is essen- 
tially a high pass filtered image. By adding the high pass 
filtered image to the original image the high frequency en- 
hancement is achieved. 


PROGRAM NAME: IMGRID (Harvard U.) 

LANGUAGE: FORTRAN IV 

COMPUTER: CDC Cyber 74, NOVA 2 (Data General) 

SPECIAL PERIPHERALS: Disk Files 


PURPOSE:: IMGRID is a general purpose geographic data base manipulation 

program developed by the Harvard University Graduate School of 
Landscape Design, This program provides for manipulation of 
multiple data variables related to the same geographic area 
and in a gridded format. Spatial searching, statistics genera- 
tion, and modeling via multivariable weighting parameters are 
key co its analysis capability. Visual and environmental impact 
analysis are two uses of such a system, 

PROGRAM NAME: TOPO 

LANGUAGE: FORTRAN IV 

COMPUTER: Data General NOVA 2 

SPECIAL PERIPHERALS: Tape Drive 

Image Display 


PURPOSE: After breaking down the NCIC topographic tapes into NOVA 4096 

record blocks with the CDC Cyber 74, TOPO unpacks and displays 
the data on the imaging system. Various scales may be represented 
by selection of parameters in the program. 


PROGRAM NAME: 


SCORECARD 


LANGUAGE: FORTRAN IV 

COMPUTER: Data General NOVA 2 

SPECIAL PERIPHERALS: Tape Drive 

PURPOSE: SCORECARD performs a maximum likelihood classification on 

specific polygons within a data set. This is used to eva- 
luate accuracy of classification by comparing the classi- 
fied data for test fields to known ground truth. 


PROGRAM NAME: THERMAL 

LANGUAGE : FORTRAN 

COMPUTER: Data General NOVA 2 

SPECIAL PERIPHERALS: Tape Drive 

Image Display 


PURPOSE: THERMAL is a program designed to unpack and display digital 

thermal data from a NASA owned thermal scanner - RS18. 


PROGRAM NAME: CHAN 24 

LANGUAGE: FORTRAN IV 

COMPUTER: Data General NOVA 2 

SPECIAL PERIPHERALS: Image Display 

Tape Drives 

PURPOSE: CHAN24 unpacks data from the Bendix 24 channel aircraft 

scanner and reformats the data such that single channels 
may be accessed and written out to either tape or an 
image display. 
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PROGRAM NAME: 
LANGUAGE: 

COMPUTER (S) : 

SPECIAL PERIPHERALS: 


M60CI- 

FORTRAN IV 

Data General NOVA 2 

2 Tape Drives 
Disk 


PURPOSE: M60CL uses a maximum likelihood decision rule to classify 

Landsat data into one of up to sixty classes for which 
means and covariances are available on a disk file. Thresh- 
olds (probability of correct classification) are output for 
each record of data. This program is a record by record 
classifier. One record is read from tape, classified, and 
then written to an output tape before the next record is 
processed. 


PROGRAM NAME: 
LANGUAGE : 
COMPUTER: 


ASTEP 
FORTRAN IV 

CDC Cyber 74, UNI VAC 1108 


SPECIAL PERIPHERALS: Tape Drives 

Disk Files 


PURPOSE: ASTEP is a general purpose earth resources analysis program 

developed by TRW Systems for NASA Johnson Space Center. The 
acronym ASTEP stands for Algorithm Simulation. Test and 
Evaluation Program. The program is of modular construction 
with standardized input-output such that they are essentially 
transparent to the user. Different classification, clustering, 
statistics generating, or feature selection algorithms may be 
tested against one another with a minimum of programming change 
to the whole system. In addition to its usefulness as an al- 
gorithm test bed, ASTEP has been used effectively as an opera- 
tional, interactive classification system for Landsat data. 
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PROGRAM NAME: 
LANGUAGE: 

COMPUTER: 

SPECICAL PERIPHERALS: 


SEARCH (NASA/EKL - Ronnie Pearson) 

FORTRAN IV 

SEL 32/55 

1 Tape Drive 
Disk Files 


PURPOSE/DESCRIPTION: SEARCH is a program developed by Ronnie Pearson of 

NASA/ ERL for unsupervised development of signatures 
for use in a maximum likelihood classification scheme. 

A 3 k 3 or 6 x 6 pixel moving window is used in a 
single iteration through the raw Landsat data to form 
candidate signatures. A maximum number of acceptable 
signatures is specified and a divergence criterion is 
used for merging, splitting, and selection of signa- 
tures. This program normally takes approximately one 
hour for development of signatures for one Landsat 
scene. Auxiliary programs are available for intui- 
tively assigning color values for each class for use on 
a color display based on a two dimensional plot of the 
signature means for channels 2 and 4 of Landsat MSS 
data. 


PROGRAM NAME: 
LANGUAGE: 
COMPUTER (S) : 


COORD 

FORTRAN IV 
CDC Cyber 74 


SPECIAL PERIPHERALS: Disk Files 


PURPOSE: COORD accepts pairs of Latitude-Longitude or UTM coordinates 

and Landsat pixel coordinates for Ground Cortrol Points (GCP) 
and computes a least squares fit of the transformation matrix 
needed to map Landsat data into a standard coordinate system. 
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PROGRAM NAME: 


CLUSTER 


5v 


LANGUAGE: FORTRAN IV 

COMPUTER (S): Data General NOVA 2, CDC Cyber 74 

SPECIAL PERIPHERALS: 2 Tape Drives 


PURPOSE: CLUSTER is a sequential clustering algorithm which creates 

an unsupervised classification of Landsat or aircraft multi- 
spectral scanner (MSS) data using a Euclidean distance measure 
as a decision criterion. This system decides how many "different" 
types of land cover there are in a MSS scene. This system is 
dependent on user input parameters which specify the criteria 
for number of clusters, merging, creation of new clusters, and 
exclusion of clusters. This technique is often used to define 
training fields for supervised classification. 

PROGRAM NAME: CLUST (NASA/ERL - Ronnie Pearson) 

LANGUAGE: FORTRAN IV 

COMPUTER: SEL 32 /55 

SPECIAL PERIPHERALS: I Tape Drive 

Disk Files 


PURPOSE/DESCRIFTION: CLUST is structurally similar to the SEARCH program 

for unsupervised development of class signatures. 
CLUST, however, uses a Euclidean distance measure 
in its sort and merge control logic for clusters. 
CLUST is often more useful than SEARCH in very broken 
terrain where fields of 40 acres are not common. 
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PROGRAM NAME: 
LANGUAGE: 

COMPUTER: 

SPECIAL PERIPHERALS: 


SUPERG (NASA/ERL - Marcellus Graham) 

FORTRAN IV 

SEL 32/55 

2 Tape Drives 
Disk Files 


PURPOSE/DESCRIPTION: SUPERG is a rectification, program for Landsat MSS data 

which includes 2nd order mirror corrections into a 
least squares determination of a covariance matrix for 
conversion of Landsat pixels into a UTM coordinate 
system. The data are resampled along scan lines to 
satisfy scaling considerations and the output file 
contains new pixels which are directly related to 
the UTM system. Between 10 and 30 Ground Control 
Points are suggested for complete determination of 
the transformation. Rotation of the data to true 
North is not accomplished hy this program. 


PROGRAM NAME: 
LANGUAGE: 

COMPUTER: 

SPECIAL PERIPHERALS: 


TRAIN 

FORTRAN IV 

Data General NOVA II 

COMTAL Interactive Color Video Display (3 images) 
Magnetic Tape Drive 
Disk Files 


PURPOSE: TRAIN is an interactive training field selection and statistics 

generation program for Landsat digital analysis. In conjunction 
with a video display system with a cursor or joystick, a subset of 
a Landsat image may be selected by drawing an arbitrarily shaped 
polygon (up to 100 vertices) around an area on the display screen. 
The program keeps track of position on the input Goddard format 
GCT and calculates the normal statistics (mean and covariance) 
of the selected training fields. Histograms of the multivariate 
distribuitons are displayed on the display screen and the mean 
and variance, polygon vertices, and histograms may be saved on 
a disk file. The program acts in a question-answer mode which 
requires no knowledge of computer languages by the user. 

A parallelopiped classifier is also implemented in this pro- 
gram which will classify a 256 x 256 element scene for one class 
in near real time. 
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G jk 

C ADDS1GM * 

G * 

G THIS ROUT DTE ADDS SPECIFIED SIGNATURES * 

C TO FORM NEW MEM MD COVARIANCE * 

C * 

G SEQUENCE: ADDSIGM INPUT SIGFIL NEWS I G * 

C * 

G » 

C CREATED AT GEORGIA TECH EES 

C 

G PROGRAMMER: NICKOLAS L. FAUST 

C * 

G ^wjc^jc ^C'!- *£v^i' I* i|i *4* 3 Ji i *|* J fi if C 5)C ifi * 1 1 

DIMENSION COV1 (4, 4) ,C0V2(4,4> , NOTH ( 12) ,NUH2( 12) , I0RDER(60) 

DIMENSION IDC 6) ,AHEAN(4,60) , BCOV( 4, 4, 60) , 11(30) , 12(30) , 13(30) 
DIMENSION NP( 60) ,NAMI(60) ,NAM2(60> ,NAM3(60) ,NAM4(60) 

DIMENSION VMlt 12) f VM2( 12) ,VM(4> , ISW(2) 

COMMON/NPF/NP 

ND=4 

IP= 12 


C 

G COMARG AND OPEN STATEMENTS 

C 

CALL OPEN! 1 , "COM. CM" , 1 , 1 ERR) 

CALL COMARG( 1 , 11, ISW, IERR) 

CALL COMARG( 1,11, ISW, IERR) 

CALL COMARG( 1, 12, ISW, IERR) 

CALL COHARG( 1 , 13, ISW, IERR) 

CALL OPEN ( 2 , 11,0, IE) 

CALL FOPEN ( 4 , 13, "B”) 

G 

C 

TYPE » INPUT NEWSIG NAME " 

READ( 11,200)N1,N2,N3,N4 
READ(2) NSIG 
WRITE! IP) NSIG 
C 

G GET SIGNATURES LISTED IN INPUT FILE 

G 

CALL GSIGC AMEAN, BCOV, NS IG, 2 , 3 , 12, I ORDER, NAM1 , NAM2 , NAMQ , NAM4) 
NUMIC 1 ) ~NP( 1) 

DO ( J= 1,4). 

: VM1(J)=AMEANCJ, 1> 

: DO (L= 1 , 4) 

: : C0V1 ( J , L) =BCOV< J , L, 1 ) 

: : . . FIN 

: . .FIN 
NS1=NSIG~1 
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LOOP TO COMBINE NSIG SIGNATURES 

DO ( K- 1 , NS 1 ) 

: KP1=K+1 

: DO !J=1,4) 

: : VK2! J)=AHEMC J,KPI) 

: : DO !L=1,4> 

: : : C0V2C J,L) “BGOVC J,L,KP1) 

: : i . .FIN 

: : . .FIN 

: NUM2! 1)=NP!KP1) 


CALL ADDS IG! COVI , C0V2 , VM1 , VM2 , ND, NUJfl , NUMB) 


. . .FIN 
DO ( Is 1,4) 

.* VM! I) =VM1C I) 

- FIN 

WRITE! IP , 20 1 ) N l , N2 , N3 , R4 
WRITE! IP) NUUK 1) , VII, COVI 
DO ( L= 1 , 6 ) ID( L) =0 
WRITE BINARY! 4) ID 
WRITE BINARY! 4) N1,N2,N3,N4 
WRITE B INARY! 4) NUK1 ! 1 ) , VM, COVI 

200 FORMAT! 4A2) 

20 1 FORMAT! 2.X, 4A2 ) 

STOP 

END 
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C *{* 

C ALARM < SUBROUTINE) i 
C jft 
C THIS SUBROUTINE ALARMS 1 CLASS ON THE COMAL * 
C CORRESPONDING TO THE LAST SIGNATURE, * 
C IT ALARMS TO GRAPHICS OVERLAY #IOV * 
C * 


P »A» *.t* »J. »J> *±* vj* tl* «L* «J» *J • »L**L* i. *U < L.^1> -J f 

C 

C CREATED AT GEORGIA TECH EES 

C 

C PROGRAMMERS: NICKOLAS L. FAUST 

C ROBERT A. MADDOX 

C 


* 


X 

:lt 


X 

X 


cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


SUBROUTINE ALRM2C IOV) 

DIMENSION IMINC4) , IMAXC4) , 1X2(0:255,3) 
COMMON/DTRARS/ IMAG1 ( © : 255 > , I XDC 0 : 5 1 1 , 4) , IHAG( 0:511) 
COMKON/HXSTX/ICOUNT(4, 100) 

EQUIVALENCE ( IXD, 1X2) 

DO (1=0,15) IMAG( I) =0 

DO ( I=©,255) CALL GH( IOV, I, I MAG, 16) 

DO ( 1=1,4) 

: DO (J= 1,100) 

: : IF( ICOUNTC I , J> , GT. 0) GO TO 30 

* s * * f irr 

30 : CONTINUE 

: IKINC I) = J 

: DO (JJ=1, 100) 

: : N= 101-JJ 

: : IFC ICOUNTC I ,N) . GT. ©)G0 TO 50 

• * ^ p J]\J 

50 : CONTINUE 

: IMAXC I) =N 

: . .FIN 

IMINl=miN( 1) 

IMIN2=miNC2) 
miN4= IMIN( 4) 

I MAXI = IMAXC 1) 

IMAX2=IMAX(2) 

IMAX4= IMAX( 4) 

DO ( K=0, 255) 

: CALL IMRDC 0 , K, IMAGC 0) , 128) 

: CALL IMRD( 1 , K, IMAGC 128) , 128) 

: CALL IMRD(2,K, IMAGC256) , 128) 

: CALL UPAC8C IMAG, 1X2(0, 1), 384) 

: DO (L=0,255) 

: : IMAG1(L)=0 

: : IF ( IX2CL, 1) .LT. IMIN1) 

: : IF C IX2CL, 1) .GT. IMAX1) 

: : IF C IX2CL,2) .LT. IMIN2) 

: : IF ( IX2(L,2) .GT. IMAX2) 

: : IF ( IX2CL,3) .LT. IMIN4) 

: : IF ( IX2(L,3) .GT. 01AX4) 

: : IMAGK L) = 1 

70 * : CONTINUE 

: : . , F IN 

: CALL PACBC IMAG1 , IMAG, 16) 

: GALL GWRC IOV, K, IMAG, 16) 

‘...FIN 
RETURN 
END 


GOTO 70 
GOTO 70 
GOTO 70 
GOTO 70 
GOTO 70 
GOTO 70 


ORIGINAL 
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G 

G ANGDIS (SUBROUTINE) 

C 

G PURPOSE: 

C COMPUTES THE DISTANCES AND ANGLES BETWEEN A SET OF 

G VECTORS OF ARBITARY DIMENSION 

G 

G DESCRIPTION OF PARAMETERS 

G 

C INPUT: 

C CALL IKG SEQUENCE 

G VM - NVM VECTORS OF DIMENSION ND , STORED BY COLUMNS 

G NVM - NUMBER OF VECTORS 

C ND - DIMENSION OF EACH VECTOR 

G IDISF - . Ed. 1 COMPUTE EUCLIDEAN DISTANCE 

C . NE. 1 COMPUTE LI DISTANCE 

C OUTPUT 

C CALLING SEQUENCE 

C R - AN NVH BY NVM MATRIX WITH I TH - J TH ELEMENT 

G CORRESPONDING TO I TH - J TH VECTORS AND EQUALS ANGL 

C IF ABOVE DIAGONAL AND EQUALS DISTANCE IF BELOW 

C DIAGONAL 

C 


* 

* 

* 

* 

* 

* 

* 

a: 

* 

Jft 

* 

* 

* 

* 

* 

* 

kU 

Jfc 

* 

m 


C 

C CREATED AT NASA/ JSC (ASTEP) * 

C * 




SUBROUTINE ANGDIS( VM, NVM, ND, IDISF, R) 
D I MENS I OR VMC ND , NVM) , R( NVM, NVM) 

N = NVM - 1 
DO 20 J = 1,N 
R( J, J) = O. 

IX = J * 1 

DO 10 I = II, NVM 

IFC IDISF. EQ. 1) GO TO 6 


D“0. 


DO 4 10= 1 , ND 

4 D“D+ ABSC VMC K, J) -VM( K, I ) ) 

GO TO S 

6 CALL EDISTC VMC 1 , J) , VMC I , I) , ND, D) 
8 CONTINUE 

CALL ANGLEC VMC 1 , J) , VMC 1 , I ) , ND , A) 
R< I , J) = D 
RC J, I) = A 
10 CONTINUE 
20 CONTINUE 

RC NVH, NVM) = 0. 

RETURN 


END 
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ANGLE ( SUBROUTINE) 

PURPOSE: 

COMPUTES THE ANGLE BETWEEN TWO VECTORS OF ARB IT ARY 
DIMENSION 

DESCRIPTION OF PARAMETERS 
INPUT 

CALLING SEQUENCE 
VI - 1 ST VECTOR 
V2 - 2 ND VECTOR 
ND - DIMENSION OF VI AND V2 

OUTPUT 

CALLING SEQUENCE 

A - ANGLE BETWEEN VI AND V2 IN DEGREES 


* 

* 

& 

:j: 

h; 

*2# 

:ii 

* 

» 

H: 


it* 

«r* 


it* 

A* 

A* 


c 
c 
c 

G 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

C 

C CREATED AT NASA/ JSC CASTEP) ;/c 

C ^ 

SUBROUTINE ANGLE C VI, V2. ND, A ) 

DIMENSION Vl(ND) ,V2(ND) 

Dl = 0. 

D2 = 0. 

A - 0. 

DO 10* I = 1,ND 
Dl = DI + VI (I)**2 
D2 = D2 + V2(I)**2 
10 A = A + VIC X)*V2( I) 

IF ( DI.EQ.0.0 .OR. D2.EQ.0.0 ) GO TO 20 
A = A/( SQRTC D I ) #SQRT( D2) > 

A = 57.29578*AC0SC A ) 

15 RETURN 
20 A = 90.0 
GO TO 15 
END 
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/ 1 *i* *J'*>* , 4**J'' Q-* •>*]-* a<*UvU«ua<.j/ 

>h 

autos nz * 

<jj 

DESIGNED TO OUTPUT CLASSIFIED TAPE IN PATTERN FORMAT * 

TO DISPLAY ON THE VERSATEC * 

* 

SEQUENCE: AUTOS UZ MTU: F INPUTFILE * 

* 

CREATED AT GEORGIA TECH EES * 

Jj/ 

PROGRAMMER: MICHAEL D. FURMAN * 

* 


INTEGER IHAGC 180) , ITAPEC 10) , IOUTC360) , I0RDERC25, 15) , ISWC2) 
INTEGER NCC25) , IWORKC 4100) , IMAG3C 180) , IMEMC 0:255) , IMAG2C360) 
INTEGER IWORK1C 360) , IWORK2C360) , IWORK3C360) , IW0RK4C360) 
INTEGER IWORK5C36G) , IWORIC6C36G) , IW0RK7C360) , I WORKS ( 360) 
INTEGER IW0RK9C360) , IWORK10C360) , IWORK1 1(360) , IWGRK12C360) 
INTEGER IMES( 15 , 40) , IFLDC 10) 

EQUIVALENCE < IWORKC 1) , IWORK1C 1) ) 

EQUIVALENCE ( IWORKC 359 ) , IWORK2C 1) ) 

EQUIVALENCE ( IWORKC 7 17) , IW0RK3C 1) ) 

EQUIVALENCE ( IWORKC 1075) , IW0RK4C 1) ) 

EQUIVALENCE C IWORKC 1433) , I WORKS < 1) ) 

EQUIVALENCE ( IWORKC 1791) , IW0RK6C 1) ) 

EQUIVALENCE C IWORKC 2149) , IW0RK7C 1) ) 

EQUIVALENCE < IWORKC 2507) , IW0KK8C 1) ) 

EQUIVALENCE C IWORKC 2865) , IW0RK9C 1) ) 

EQUIVALENCE C IWORKC 3223) , IW0RK10C 1) ) 

EQUIVALENCE C IWORKC 3531) , IWORKUC 1) ) 

EQUIVALENCE C IWORKC 3939) , IW0RK12C 1) ) 

CQMMON/RUM/ ID( 165) , IWORK 
DATA ID/ 


10 , 0 , 0 , 0 , 0 , 6 , 0 , 0 , 0 , 0 , 0 , 

20, 176K, 102K, 102K, 102K, 102K, 102K, 102K, 102K, 176K,0, 

30,0, 176K, 176K, 176K, 176IC, 176K, 176K, 176K,0,0, 

4300K, 340K, 160K, 160K,6GK, 30K, 14K, 16K, 16IC,7,3, 

5 176IC, 74K, 30K, 0 . 20 IK, 3G3IC, 20 IK, 0 , 30K, 74K, 176IC, 

6303K, 303IC, 303IC, 303K, 0,0,0, 303K, 303IC, 303K, 303K, 

730IC, 30IC, 30K, 30K, 30IC, 30K, 30K, 30K, 30K, 30K, 30K, 

B30IC, 30K, 30K, 30K, 377K, 377K, 377K, 30K, 30K, 30K, 30K, 

93,7, 16K, 16K, 14K,30K,60K, 160K, 160K, 340K, 300K, 

1 377K, 30 IK, 24 IK, 23 IK, 23 IK, 23 IK, 23 IK, 23 IK, 205K, 203K, 377K, 
120 IK, 303IC, 347K, 176K, 74K, 30K, 74IC, 176K, 347K, 303K, 20 IK, 
2374IC, 37©K, 36 IK, 36 IK, 363K, 347K, 3 17K, 2 17K, 2 17K, 37K, 77K, 
377IC, 37K, 217K, 2 17IC, 317K, 347K, 363K, 36 IK, 36 IK, 370K, 374K, 
4377K, 377IC, 377IC, 303K, 303IC, 303IC, 303IC, 303K, 377IC, 377IC, 377K, 
5377K, 377K, 377K, 377K, 377K, 377K, 377K, 377K, 377K, 377K, 377K/ 
CALL OPENC l , "COM. CM" , 1 , I ERR) 

CALL COMARGC 1, ITAPE, ISW, I ERR) 

CALL COMARGC 1, I TAPE, ISW, I ERR) 

CALL COMARGC 1 , IFLD, ISW, I ERR) 

CALL F0PENC3, IFLD, "B"> 
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10=3 

CALL MTOPDC 2 , ITAPE , 0 , IE) 

DO (11=1,25) 

: DO (12=1,15) IORDERC II, 12) =0 

: . .FIN 

TYPE "DOES INPUT TAPE HAVE THRESHOLD? ( 2= YES, l=N0) " 

READ( 10) ITH 
TYPE "LINE" 

READ( 10) ILINE 

TYPE "MAP PORTION (1 - 12) ",IANS2 

READ( 10) IANS2 

TYPE "LENGTH OF DATA " 

READC 10) 1ST 
IENB= IEL+359 

IF C IEND. GT. 4100) IEND=4100 

TYPE "PRODUCE OVERALL MAP? ( 1=YES) '* 

READC 10) IANS 

TYPE "NUMBER OF GENERALIZED CLASSES " 

HEAD( 10) I GEN 
DO ( 1= 1 , IGEN) 

: TYPE "INPUT CLASS DESCRIPTION " 

*. KEAD( 10,101) C IMESC 1,10 ,K= 1,20) 

: TYPE "NUMBER OF SUBCLASSES " 

: READ( 10) NC( I) 

: NC1=NC( I) 

: TYPE "INPUT SUBCLASSES " 

: READ( 10) ( IORDERC 1,10 , IC= 1 , NCI) 

: IF (I.NE.IGEN) TYPE "NEST GROUP" 

: . .FIN 

50 DO (1=1,253) IMAG2(I) = 0 
DO (1=1,180) IKAG3( I) = -' 1 
DO (1=0,255) IMEPK I ) =0 
NRITEC 12) 

L2= 154 

DO ( 1=1, ILINE) GALL MTD IOC 2,0, IN0RK2, IS, IE, IG> 

IF (IANS.EQ. 1) 

: DO (1=1,180) IOUT( I) =0 

; DO ( M5= 1 , IGEN) 

: : L3=NCCM5) 

: : DO (J2=I,L3) IMEMC IORDERC M5 r J2) ) =L2 

: : URITEC 12, 100) ( IMES< 1*15,10 ,IC= 1,20) ( IORDERC MS, I) , 1=1, L3) 

: : DO (J=l,ll) 

: : : IOUT( 15)=ID(L2+J) 

: : : GALL MTXC IOUTy 180) 

: : : . . F IN 

: : L2=L2- 1 1 

: :..FIN 

: NRITE( 12) 

:..FIN 
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DO ( PI1“ 1 , I GEN) 

: IF ( IANS. KE. I> 

: : L3=HC( Ml ) 

: : URITEC 12, 100) C IMES(M1,K) ,K=1,20) C IOBDERCM1, Z) , I=1,L3> 

: : DO C I- 1 , 1G0> I OUT! I) =0 

: : DO (1=0,255) IMEPI(I)=0 

: : DO ( J= 1 , 1 1 ) 

: : : L2=154 

: : : DO (J2=1,L3> 

i : : : IMEHC I ORDER( Ml , J2) ) =L2 

: : : : IOUT( J2*4+9) = IDCL2+J) 

: : : : L2=L2-11 

: : : J..FIN 

: : : CALL MTX( IOIJT, 180) 

: : : . .FIN 

: : WRITE! 12) 

: : . . F IN 

: DO C 1=1,8) CALL MTXC IIIAG3, 180) 

: DO ( M2= 1 , 1ST) 

: : DO (1=1, ITH) 

: : : CALL MTDI0C2 , 0 , IWOKK, IS, IE, IC) 

: J i IF (IC.LT.5) GOTO 55 

: : :..FIN 

: : CONDITIONAL 

: : : Cl ANS2 . EO. 1 ) CALL AUTEC C IW0RK1 , IMEM, 1 1 ) 

: •' : C IANS2.E0. 2) GALL AUTEC( IW0RK2, IMEM, 11) 

: : : ( I ANS2.E€t. 3) CALL AUTECC I WORKS, IMEM, 1 1 ) 

: : : ( IANS2.EQ..4) CALL AUTECC IW0RK4, IMEM, 11) 

: : : (IANS2.E0.5) CALL AOTECC IW0RK5 , IMEM, 1 1) 

: : : (IANS2.EGL6) CALL AUTEC( IW0RK6, I MEN, 11) 

: : : ( IANS2.EQ..7) CALL ADTEC( IW0RK7, IMEM, 1 1) 

: : : CIANS2.E0.8) CALL AUTECC IW0RK8, IMEM, 11) 

: : : ( IANS2.EQ..9) CALL AUTECC IW0RK9, IMEM, 11) 

: : : ( IANS2.EGL. 10) CALL AUTECC IWOEKIO , IMEM, 1 1) 

: : : < IANS2 . EO. 1 1 ) CALL AUTECC IWORK1 1 , IMEM, 11) 

: : : C XANS2.E&. 12) CALL AUTEC( IW0RK12 , IMEM, 1 1) 

: : : . .FIN 


J..FIH 

DO C 1=1, 11) CALL MTXC IMAG3, 180) 

DO (1=1,700) 

: BO (J= 1,360) IOIJTC J) = IDC 1 ) 

: DO <N10=1,1100) N11=N10 

: CALL PAC8C IOUT, IMAG,360) 

: CALL MTX( IOUT* 180) 

: . .FIN 

CALL MTDI0C2, 100O0K, IWORK, IS, IE, IC) 

DO (1=1, ILINE) CALL HTDI0C2, 0, IWORK2, IS, IE, IC) 
IF CIANS.EQ.. 1) 

: IANS=0 


: : GOTO 50 

: :..FXN 

: . .FIN 

DO (1=1,10) WRITEC12) 

STOP FINISHED 

FORMAT! IX, / , IX, 20A1 , / , ,r CLASSES ", 2014) 
FORMAT! 20A1 ) 

END 
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C jjc 

G CHAN24 * 
C '•!• 
C UNPACKS AND REFORMATS BEND IX 24-CHANNEL DATA FOR DISPLAY * 
C ^ 
C SEQUENCE CHAN24 MTU:F MTU:F * 
C 


C 

c 

C 

C 

C 

c* 


CREATED AT GEORGIA TECH EES 
PROGRAMMER: MICHAEL D. FURMAN 

INTEGER IDATAC 4100) , IMAGC512) , IMAG2<£56) , ITAPEH5) , TTAPE2C5) 
INTEGER I WORK! 256) , I EE ADC 30) , IX(3) , 

COMMON IX, IY 

EQUIVALENCE C IXC 1) , 1X1) , C IX< 2) , 1X2) 

CALL OPENC I , " COM. CM" , 1 , I ERR) 

CALL COMARGC 1 , ITAPE1, ISW, I ERR) 

CALL COMARGC 1, ITAPE1, ISW, I ERR) 

CALL COMARGC 1 , ITAPE2 , ISW, I ERR) 

DO (1=1,30) IHEADC I) =0 


IYC3) , ISWC2) 

, C IY( 1) , IY1) , C IY(3) 


IY3) 


10 FORMATt IX, 12, "/", 12, V", 12) 

1 1 FORMATt IX, 11 INPUT TAPE NO . “ , Z) 

12 FORMAT! IX, 10A2) 

13 FORMAT! IX, "BLOCK SIZE= ", 13) 

14 FORMAT! 10A2) 

TYPE "DATE H,B,Y" 

ACCEPT IBEAD! 3) , IHEADC 4) , IHEADC 5) 

TYPE "INPUT TAPE NUMBER? XXXX" 

ACCEPT IHEA1X1) 

TYPE "COMMENTS ON RUN — 2© CHARACTERS" 
READ! 11, 14) ( I HEAD! I) , 1=10,20) 

CALL COLORS!® 

TYPE "FAST SCAN? (2= YES, 1=H0> " 

ACCEPT IFSCAN 

TYPE "DATA BLOCKS TO SKIP?" 

ACCEPT ISKB 

TYPE "INPUT BLOW-UP FACTOR” 

ACCEPT IBLUP 
IBLUPM1= IBLUP- 1 
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* * * w *X ■:< * 


OF rn’m L PAGB Is 

* P(x ® QUAUIY 


IELEM= 1 

IF C IFSCAN. EQ. 1) 

. TYPE "START WITH ELEMENT?" 

; ACCEPT IELEM 
: . .FIN 

CALL 5CALESUR 

TYPE " INPUT CHANNEL NUMBER! 1-24) " 

ACCEPT I CHAR 

I SKB= ( ICHAN+2)/3+( ISKB*9)+2 
I NUM= 5+ 1 CILAN- ( ( ( ICHAN+2 ) /3+ 1 ) *3) 

L= INUK*393+2+ IELEM 
LAST= L+256+ ( ( IFSCAN- 1 ) *92) 

IF C IELEM. GT. 94) LAST=35 l+( INUM&393) 

CALL MTOPDC 3 , ITAPEl , 0 , IER) 

CALL KT0PD(4, ITAFE2,0. IER) 

DO (1=1, ISKB) CALL MTDIQC 3,0, IDATA, IS, IER, ICNT) 

M2=512/IBLUP 

K3=0 

DO (111=1, M2) 

: DO ( N= 1 , IFSCAN) 

: : CALL MTD 10(3, 30010IC, IDATA, IS, IE, IC) 

: : CALL KTDI0(3,0, IDATA, IS, IE, IC) 

: : IF (IC.LT. 10) GOTO 150 

: : , . F IN 

: J=1 

: DO (IC=L, LAST, IFSCAN) 

: : IKAG( J) = IDATA( K) 

: J J=J+1 

: : . . F IN 

: BLOWUP-AREA 

: DO ( 1= 1, IBLUP) 

: : CALL IKWRITEt 0, M3, IMAG2, 256) 

: : M3= M3+ 1 

: : . . FIR 

: . .FIN 

150 DO ( Ll= 1,3) 

: PAUSE — POSITION CURSER AND HIT RETURN 

: CALL RTARG( IX(L1) , IY( LI) ) 

: . .FIN 
IXl=IKl/2 
1X2= 1X2/2 
DO ( 1=1 ,3,2) 

: CALL IMREAD(9, IY( I) , I WORK, 256) 

: DO ( 12=1X1, 1X2) IWOKKC 12) =255 

: CALL IMWRITE(0, IY( I) , IW0RK.256) 

: . .FIN 

BO ( 1= IY1 , IY3) 

: CALL IMRE ADC 0,1, IWORK, 256 ) 

: IWORXC IXI)=255 

: IWORICC 1X2) =255 

: CALL IMWRITE( 0,1,1 WORE, 256 ) 

: . .FIN 
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TYPE "PROPER RECTANGLE? ( 1=YES; 2= NO) " 

ACCEPT IANS 

IF CIANS.Ea.2) GOTO 150 

CALL KTDIOC 3 , 1000OK, IDATA, IS, IE, ICNT) 

I X 1 = IJU* I FSCAN+I ELEM- 1 
I X2= 1 252* IFSCAN+ IELEH- 1 
I Yl= IY1*IFSCAN 
IY3= IY3&IFSCAN 
IPAS=9*IY1 
IDIF- ( IY3-IY1) 

IDIF2=IK2-IX1 
WRITE- IIE ADER 

CALL MTD 10(3, 3000 1K+IPAS, IDATA, IS, IE, ICNT) 

DO (1=1, IDIF) 

: CALL MTD 10(3, 300O 1 K, IDATA , IS , IE , ICNT) 

: DO ( K5= 3 , 8) 

: : CALL MTDJOC3,0, IDATA, IS, IE, ICNT) 

: J DO (L6=3, 739,393) 

: : : J= 1 

: : : DO CL5=IX1,IK2) 

: : : : IWORK( J) = IDATA( L5+L6) 

: : : : J=J+I 

: : : i . .FIN 

: : : CALL KTDIO( 4, 5O00OK+ IDIF2, I WORK, IS, IE, ICNT) 

: : S..FIN 

: : . .FIN 

: . .FIN 

DO (1=1,10) CALL NTDIO( 4, 60000K, IWORK, IS, IE, ICNT) 
CALL MTDIO( 4, 10000K, IWORK, IS, IE, ICNT) 

CALL MTDI0C3, 10000K, IDATA, IS, IE, ICNT) 

STOP 

C 

TO WRITE- HEADER 
: WRITEC 12 , 11) 

: WRITEC 12) I HE AD ( 1) 

: NRITE(12,10) IHEAD(3) , IHSAD(4) , IHEAD(o) 

: WRITEt 12, 12) C IHEAD( I) , 1=10,20) 

: WRITEC 12, 13) IDIF2 

: CALL fITDlO( 4, 50036K, I HE AD, IS, IE, ICNT) 

J..FIN 

C 

TO BLOWUP- AREA 
: DO ( K7=0, 255) 

: : K9=256~KT 

: : KDB = K9 / 1 BLUP 

: i I IL4GC K9 ) = IMAG( KDB) 

: : . .FIN 

: CALL PAC8< I MG , IMG2 , 5 12) 

: . .FIN 

200 STOP FIELD ERROR — CHAN24 MTU:F MTU:F 
END 
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C 

G 

C 

G 

C 

G 

C 

G 


CLASIFY 

THIS ROUTINE CLASSIFIES POLYGONS WITH A MAXLIK GLASS IF IEE 
A TABLE LOOK-UP FEATURE IS INCLUDED 

SEQUENCE: CLASIFY INPUT DATAP MAXTAP SIGFIL 


C 

C CREATED AT GEORGIA TECH EES 

C 

G PROGRAMMERS: NICKOLAS L. FAUST 

C ROBERT A. MADDOX 

G 

D I MENS ION AMEANC4.60) ,ECOVC4,4,60) , IVX( 10 1> , IVY( 101) ,DETC60) 

1 , LA( 133) , GOUNTC 0:61) , IBUF( 1640) ,CLC0N(61) 

2, NAM1(60) .NAM2C60) ,NAM3(60) ,NAM4(60) ,B1(4,4) ,B2C4,4) 

DIMENSION JBUFC 4 ,810) , FIELD! 17) , ISWSC2) , IORD(0:61) 

DIMENSION DUM(6l) ,DUM1(4,4,6©) ,BUM2(4,60) , ITABLC5,300> 

COMMON JBUF 

EQUIVALENCE ( JBUF , BUM) , ( JBUF , DUM1 ) , C JBUF , DUM2) , ( JBUFC 1,20), DET) , 
1C JBUF, Bl) ,( JBUFC 1, 10) ,B2> 


G 

C 

G 


G 

C 


CQMARG CALLS AND OPEN STATEMENTS 


CALL OPEN! 1 , "COM. CM” , 1 , IE) 
GALL COMARGC 1, FIELD, ISWS, IE) 
C ALL COMARGC I , FIELD , ISWS , IE) 
GALL OPEN! 2, FIELD,©, IE) 

GALL COMARGC 1, FIELD, ISWS, IE) 
CALL KT0PDC3,FIELD,0, IE) 

CALL COMARGC 1 .FIELD, ISWS, IE) 
CALL HT0PDC5,FIELD,0, IE) 

CALL COMARGC 1 , F I ELD , ISWS , IE) 


GET PROGRAM NAME 
GET INPUT FILE NAME 
OPEN INPUT FILE 
GET DATA TAPE NAME 
OPEN TAPE FILE 
GET OUTPUT TAPE NAME 
OPEN TAPE OUTPUT FILE 
GET SIGNATURE FILE NAME 


IP- 12 
ND=4 

BNV=-1.E+12 ; BIG NEGATIVE NUMBER 

ITPTR=0 ; INITIAL TABLE POINTER 

READ! 2, 200) NSIG 
G 

C GET SIGNATURES FOR CLASSIFICATION 

C 

CALL GS IGC AKEAN , BCOV, NS IG, F I ELD, NAM1 , NAH2 , NAII3 , NAM4 , GLCON) 
DO C L= I , NSIG) 

: DO CJH=1,4) 

: : DO CJI-1,4) 

: : : B2C JH, JI) “BCOVC JH, JI,L) 

: : : . .FIN 

: : . .FIN 

: GALL SYKINVC B2 , Bl , DETP , 4) 

: DO C JH= 1 , 4) 

: : DO C JI= 1 , 4) 

: : : BCOVC JH, JI , L) =B1( JE, JI) 

: : : . . F IN 

: : . .FIN 

C : DIVIDE DIAGONAL ENTRIES OF COVARIANCE MATRIX BY TWO 

: DO ( J= 1 , 4) BCOVC J,J, L) = BCOVC J, J, L) P2. 

C : GET NATURAL LOG OF DETERMINANT 

: DETC L) = ALOGC DETP) 

: . .FIN 
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IYMIN= 10000 
IXMIN= 10000 
IYMAX=0 
WRITE! IP, 202) 

READ! 2, 200) NC 
DO ( I” 1 , NC) 

: WRITE! IP, 205) I 

: READ! 2, 206) lVKi I) , IVY( I) 

: WRITE! IP,206) IVX( I) , IVY< I) 

: IF! IVY( I) .LT. IYKIN) IYHIN= IVY( I> 

: IF! IVX! I) .LT. IXMIN) IXMIN=IVX( I) 

: IFC IVY< I) .GT. IYHAX) IYMAX= IVY( I) 

t..FIN 

I V5K NG+ 1 ) = FVX( I) 

IVY(NC+1) = IW( 1) 

NV=NC 

IDEL= IlT'IAIi- IYMIN 
ISICIP= IYIIIR+ 1 

SET UP THE AFRIORI INFORMATION EITHER FROM, 

A. THE INPUT FILE IF * IAPR’ EQUALS 1 IN INPUT FILE 
OR B. SET APRIORI FOR ALL CLASSES EQUAL TO 1 . 0 

READ (2,201) NTE ; NUMBER OF TABLE ENTRIES 

READ (2,200) ITH 
IF CITE.NE.0) 

: READ (2,310) THUS 

: THRS=ALQG( THRS) 

: . .FIN 

READ (2,200) IAPR 
IF ( IAPR. EQ. 0) 

: DO ( 1= 1 , NSIG) GLGONC I) = 1 . 

* T7 TTiT 

WRITE! IP, 410) 

DO (1=1, NSIG) 

: WRITE! IP, 420) I,NAM1( I) , NAM2( I> ,NAM3C I) ,NAM4( I) t GLGON( I) 

: . .FIN 

WRITE! IP, 400) 

CONFUTE GLASS CONSTANTS 
P ILOG= 2 . &ALOGC 2 . #3 . 141503) 

DO (1=1, NSIG) CLGONI I) =ALOG( CLCON( I))~.5*DET( D-PILOG 
CLGON( NS IG+ 1 ) = BNV 

SORT GLASS CONSTANTS FROM LARGEST TO SMALLEST 
ANB ALSO REORDER STATISTICS 

DO C J= 1 , NS I G) 

: BG = — 1 . E+ 10 

: DO ( 1=1, NSIG) 

: : IF (CLCONC I) .GT.BG) 

: : : BG=CLC0N(I> 

: : : IND= I 

: : J . .FIN 

i , FIN 
: IGHD( J) = IND 

: DBM! J) =GLGON( IND) 

: CLCGN< IND) =BNV 

: . .FIN 
IORD( 61) =6 1 
IORB<0)=0 
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REORDER GLASS CONSTANT ARRAY 

DO ( J= 1 , NSIG) CLCONC J)=DUMC J> 
CLCONC NS I G+l)= BNV 

REORDER COVARIANCE MATRICES 


SSStS ® 8 


DO ( 1=1 ,NSIG) 

: DO (J=1,ND) 

: : DO ( K> 1 , ND) 

i i . .FIN 
: , .PIN 

DO (1=1, NSIG) 

: DO CJ=1,ND) 

: : DO ( K= 1 , ND) 

: : . . PIN 


. .FIN 


DUMK J, K, I) = BCOV( J,K, IORDC I)) 


BCOV( J , K, I)=DUM1CJ,K, I) 


REORDER MEAN VECTORS 
DO ( 1=1, NSIG) 

: DO CJ=1,ND) DUM2C J, I)=AMEAN( J, IORDC I) ) 

: . .FIN 

DO (1=1, NSIG) 

: DO ( J= 1 , ND) AMEANCJ, I)=DUM3C J, I> 

: . . FIN 

SKIP RECORDS 

CALL MTDXQC3,30000K+ISKIP, IBDF, IS, I EE) 

DO (LL= 1,812) IBIIF(LL) =0 
IBBFC 1) = IDSL 
IBUFf 2) = IHMIN 

GALL MTDIO(5,50000K+810, IBIJF, IS, IER,NW) 

ISW=2 

JQ=810 

TOT=0 

T0THITS=0 

DO (1=0,61) C0UNT( I) =0 
DO ( 1=1,5) 

: DO (J=1,NTE) 

: : ITABLC I , J) =0 

: : . . FIN 

: . .FIN 

GALL FGTIMEt I HR, IMIN, ISEC) 

TYPE "LINE PROCESSING BEGAN <§" , IHR, IMIN, ISEC 
DO ( J= 1 , IDEL) 

: HITS=0 

: CALL mDIOC3,0, IBUF, IS, IE,NW) 

: CALL P0LY2C ISKIP+J-1 , ISW, IVX, IVY,NV,LA) 

: JD=0 

: DO ( K0= 1 , JO) 

: : DO ( KP= 1 ,4) 

: : t JBUF(KP,KO)=0 

: : i . .FIN 

: :..FIN 

: LA1=LA( 1) 

: DO ( Kl= 1 ,LA1) 

: : LSUB=2*K1 

: : L2=LSBB+ 1 

: : JS=LA( LSUB) 

: : JF=LA(L2) 

i : JD= JD+JF— JS+1 

: : IF(JF.GT.JQ) JF=JO 

: : DO ( K2= 1 , 4) 

: : : LL=C( JS-1) /2> *4+K2 

: : : LAST= LL+ ( JF- JS+ 1 ) *2 

: : : Jl=l 

: : : DO C II=LL,LAST, 4) 

: : : : JBUF(IC2, Jl) = ISHFTC IBDFC ID , -8) 

: : : : JBUFCK2, J1+1)=IANDC IBDF( II) ,377K) 

: : : : J1 =Jl+2 

i : : i . .FIN 

: : : . . F IN 

: : . .FIN 
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G 

C 

C 


1010 


1040 


1 

2 

2 

3 

4 


1020 


1030 


1000 

1050 

G 

G 


CLASSIFY THE POINTS BY MAXIMUM LIKLIHOOD ALGORITHM 
DO ( JPT= 1 , JQ.) 

: IF t JBUFt 1,JPT) .EQ.O) GOTO 1000 

: JB1= JBUFt 1 , JPT) 

: JB2=JBUF(2, JPT) 

: JB3=JBUF(3, JPT) 

: JB4=JBUF(4, JPT) 

: DO 1010 IC= 1 ,HTE 

: IF ( JB3.NE. ITABL13, IG) ) GOTO 1010 

: IF ( JB2.NE. ITABL(2, IC) ) GOTO 1010 

: IF C JB4.NE. ITABU4, IC) ) GOTO 1010 

: IF ( JB1.NE. ITABLt 1, IC) ) GOTO 1010 

i INDEX= ITABLt 5 , IC) 

: HITS=HIT8+1 

; GOTO 1O50 
: CONTINUE 

: M= 1 

: PR0B=BNV 

: VDl=JBl-AHEANt 1,M) 

: VB2- JB2-AMEAN( 2 , M) 

: VB3= JB3-AMEANC 3 , M) 

: VD4" JB4-AMEANC 4 , ID 

: CLTIHl- CLCONt II) - ( VB 1*VR 1*BC0V( 1 , I , M) 

: -t-VB2* C YD 1*BC0V( 2,1, II) + VB2*BC0V( 2 , 2 , K) ) 

: 4VD3*t VD1*BC0V( 3 , 1 , ID +VD2*BC0V< 3,2, M> 

: + VD3*BC0V ( 3 , 3 , ID ) 

: +VD4* C VD 1*BC0V( 4 , 1 , ID +VB2*BC0Vt 4,2, ID 

: + VD3*BC0Vt 4 , 3 , M) +VD4*BC0Vt 4 , 4 , ID ) ) 

i IF (CLTHR.LT.PROB) GOTO 1020 
: PROB=GLTHR 

i SNDEX=M 
: II=M+1 

: IF t PROB.LT. CLCONt M)) GOTO 1040 
: IF CITH.EQ.0) GOTO 1030 

: IF (PROB.LT.THRS) INDEX=61 

: ITPTR- ITPTB.+ 1 

: IF t ITPTR. GT. NTE) ITPTR- 1 

: ITABL( 1, ITPTR) = JB1 

: ITAELC 2 , ITPTR) “ JB2 

: ITABLt 3 , ITPTR) = JB3 

: ITABLt 4 , ITPTR) = JB4 

: ITABLt 5 , ITPTR) = INDEX 

t GOTO 1050 
: INBEX=0 

: IBUFt JPT) = IORDt INDEX) 

: . .FIN 


TGT-TOT+JD 

TOTH l TS = TOTH ITS+HI TS 

DO (K= 1 , JQ.) 

: IBK= IBUFt K) 

: COUNTt IBID = COUNTt IBK) + 1 

: . .FIN 

GALL FGTIMEt I HR, IMIN, ISEG) 
imiTEt 10,207) J, HITS 
TYPE ’'TIME = M HR, IMIN, ISEG 
GALL irTDIOC5,50000K+310 ( IBUF, IS, IER) 
.FIN 
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200 

201 

202 

205 

206 
207 
310 
400 
410 
420 
430 
440 
450 


WRITE! IP, 400) 

WRITE! IP. 430) 

DO ( IC= 1 , NS I G) 

: PER=! COUNT! K> *100) /TOT 

: IKK“IORD!K) 

: WRITE! IP , 440) K, KAMI ( K) , NAM2! ID , KAM3C K) , NAM4(K) , IKK 

1 : CLCON! IKK) , COUNT! K) , PER 

: . .FIN 

WRITE! IP, 450) COUNT! 61) 

FP ITE! I P ) TOTHITS , TOT 
WHITE! IP, 400) 

DO (KL= 1,810) IBUF! KL) =0 

BO ! ICL= 1 , 200) CALL MTDIO! 5 , 50000K+810 , IBUF, IS , IER) 

CALL MTDIO! 5 , 600O0K, IBUF, IS, IER) 

CALL MTDIO(5,60000K, IBUF, IS, IER) 

TYPE “TOTAL HITS, TOTAL POINTS TOTHITS, TOT 

FORMAT! 12) 

FORMAT! 13) 

FORMAT! 2X, " INPUT * OF CORNERS") 

FORMAT! 2X, “ INPUT CORNER # 11 , 13, "J, I ") 

FORMAT! 214) 

FORMAT! 2X, "LINE *“,I5, " PROCESSED THERE WERE \F8.0, " HITS") 
FORMAT! F 1 0 . 6 ) 

FORMAT! ///////////////) 

FORMAT! 2X, "NUMBER" , 2X, "NAME" , 6X, "APRIORI ") 

FORMAT! 4X, 12, 2X, 4A2 , 2X, F10 . 6) 

FORMAT! 2X, "NUMBER NAME RANKING CONSTANT # PIXELS PERCENT") 
FORMAT! 4X, I2,2X,4A2,3X, I2,2X,F10.6,2X,F8. 1 ,3X,F5. 1) 

FORMAT! //, 4X, “TOTAL POINTS NOT CLASSIFIED = 11 ,F10. 0,//) 

STOP 

END 
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C 

C CLEANUP 

C 

C USED TO CLEAN UP SIGNATURE , VERTEX , OR HIST FILES 

C 

C SEQUENCE: CLEANUP NAMES FILE IN FILEOUT 

C 

C * 

C CREATED AT GEORGIA TECH EES * 

C ?jt 

C PROGRAMMER: NICKOLAS L. FAUST * 

C * 

DIMENSION IFILC20) ,JFILC20) ,KFIL(20) ,SM<4) ,COV<4,4) , I<6> 

DIMENSION ICOUNTC 4 , 100) , IVXt 101) , IVYC 101) , ISWt2) 

CALL OPENC 1 , "COM. CM" , 1 , IERR) 

CALL COMARGt 1 * IF IL , ISW, IERR) 

CALL COMARGt 1, IFIL, ISW, IERR) 

CALL COMARGt 1 , JFIL, ISW, IERR) 

CALL COMARGt 1 , KF rL, ISW, IERR) 

CALL F0PEN(2, IFIL, "B") 

CALL FOPENt 4 , KFIL, "B") 

WRITEt 10,401) 

RE ADC 11) ICON 

CALL F0PEN(3, JFIL, "B") 

READ BINARYC 3) I 
WRITE BINARYt 4) I 
CALL FCL0SEC3) 

FOREVER 

CALL FOPENt 3, JFIL, ”B”) 

READ BINAPiY(3) I 
ISW=0 

READt 2, 100,ENB=121)NA1,NA2,NA3,NA4 
REPEAT WHILEC ISW.EQ.0) 

READ B INARYt 3 , END= 120 ) NS 1 , NS2 , NS3 , NS4 
IFt ICON. EQ. DREAD BINARYt 3) NP,SM,COV 
IFC IC0N.EQ.2) 

READ BINARYC 3) K, IYMIN, IYMAX 

DO ( IV= 1,K) READ BINARYC 3) IW, IVXC IV) , IVYC IV) 

FIN 

IFC ICON. EQ. 3) READ BINARYt 3) ICOUNT 
IFt ( NS 1 . EQ. NA1 ) . AND . t NS2 . EQ. NA2 > ) 

IFt C NS3 . EQ. NAS) . AND . C NS4 . EQ. NA4) ) 

WRITE BINARYt 4) NS 1 , NS2 , NS3 , NS4 
IFC ICON, EQ. 1) WRITE BINARYt 4) NP,SK,COV 
IFC ICON . EQ. 2) 

WRITE BINARYt 4) K, IYMIN, IYMAX 
DO C IV=1, K) WRITE BINARYt 4) IV, IVXt IV) , IVYt IV) 

.FIN 

IFC ICON, EQ. 3) WRITE BINARYt 4) ICOUNT 
ISW=1 


100 

401 

120 

121 


GALL 
.FIN 
.FIN 
.FIN 
.FIN 
STOP 

FORMATC 4A2) 
F0RMATC2X, " INPUT 


FCLOSEt 3) 


STOP - 
STOP - 
END 


S IGNATURE 


NORMAL EXIT 


SWITCH , 1 
NOT FOUND 


“SIG , 2 -VER , 3 -HIS”) 
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Jft 

CLDSTA ( SUBROUTINE) * 

5fc 

PURPOSE * 

ADAPTIVE CLUSTERING ALGORITHM, I -ST PASS THROUGH DEVELOP * 

CLUSTER MEANS USING PERIODIC ELIMINATION TESTS AND * 

MERGER TESTS, ALSO HAS FEATURES OF STRIP FORMULATION AND * 
SEQUENTIAL SEARCH FOR STRIP ASSIGNMENT, 2-ND PASS THROUG * 
DEVELOPS CLASSIFICATION MAP * 

* 

DESCRIPTION OF PARAMETERS X 


INPUT X 

GALLING SEQUENCE * 

V - VECTORS TO BE CLUSTERED x 

VM - INITIAL CLUSTER HEARS X 

ND - DIMENSION OF VECTORS X 

NV - NUMBER OF VECTORS IN V X 

NVM - NUMBER OF VECTORS IN VM X 

NVMMAX - MAXIMUM NUMBER OF VECTORS ALLOWED IN VM X 

NVG - WEIGHTS FOR CLUSTERS, I-TH VALUE IS NUMBER OF X 

POINTS IN I-TH CLUSTER * 


C, S, KP , R - CLUSTERING DISTANCE MEASURES FOR MERGERS, 
STRIP GENERATION, PRIORITY SEARCH, AND 
THRESHOLD FOR NEW CLUSTERS RESPECTIVELY 
NPC , NPT - COUNTER AND THRESHOLD FOR UPDATES TO PRIORITY 
LIST 

PL 1ST - PRIORITY LIST 

NEC.NET - COUNTER AND THRESHOLD FOR SMALL CLUSTER 
ELIMINATION TESTS 

NMIN - ELIMINATION THRESHOLD, NUMBER OF POINTS 
NMC.NMT - COUNTER AND THRESHOLD FOR MERGER TESTS 
IP ASS - PASS NUMBER OR ROUTING FLAG 

IP - PRINT FLAG, . EQ. O NO PRINT, .NE. 0 PRINT MERGERS 
AND ELIMINATION MESSAGES 

JFTP - NUMBER OF POINTS PROCESSED PRIOR TO THIS ENTRY 
RMINM, RMINV - MEAN AND VARIANCE FOR CLUSTER THRESHOLD 

DISTANCES (USED FOR IPASS=2 ONLY), CURRENT 
VALUES 

VHP, VAR - MEANS AND VARIANCES DEVELOPED DURING 2-ND PASS 
BASED UPON ACTUAL ASSIGNMENTS - CURRENT VALUES 


OUTPUT 

CALLING SEQUENCE 

VM - UPDATED MEANS 

NVM - NUMBER OF MEANS IN VM 

NVG - UPDATED WEIGHTS 

IMG - DEFINES CLASSIFICATION MAP 

TDIS - THRESHOLD ARRAY 

NPC , NEC , NMC - UPDATED VALUES 

RMINM, RMINV - UPDATED VALUES 

VHP, VAR - UPDATED VALUES 


•J* 

X 

X 

X 

X 

X 

X 

* 

X 

X 

X 

X 

X 

X 

* 

* 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

* 




CREATED AT NAS A/ JSC CASTEP) * 


X 

SUBROUTINE CLUSTA< V, VM, ND , NV, NVM, NVHMA X, NVG, C , S , RP , R, NPC , 

*NPT, PL I ST, NEC, NET, NMIN, NMC.NMT, IPASS, I MG, TDIS, 

* IP , JPTP , RMINM, RMINV , VMP , VAR) 

COMMON/ INOUT/NOUT, N IN 
INTEGER PL I ST 
INTEGER TDIS, V 

D IMENS ION VC ND , NV) , VH( ND , NVMMAX) , NVG( N VMMAXP , PL ISTC NVMMAX) , 

*VS C 24 ) , IHG( NV) , TD IS C NV) 

COMMON/D I ST/ ID I ST 

D I HENS I ON RMINMC NVMMAX) , RMINVC NVMMAX) 

DIMENSION VMP( ND , NVMMAX) , VAR( ND , NVMMAX) 

1111 FORMATC 1 ND ’,15) 

JPT=0 

NS=0 
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TESTS ON - NEST POINT, SMALL CLUSTER ELIMINATION, 
MERGING TIME, AND PLIST UPDATE TIME 


12 CONTINUE 

IFC IPASS.EGL2) GO TO 16 
NEC=NEC+NS 

IFCNEC.GE.NET) GO TO 70 
14 NMC=NMC+NS 

IF(NMC.GE.NMT) GO TO 80 
16 KPC=NPC+NS 

IF(NPC.GE.NPT) GO TO 90 

18 JPT= JPT+ 1 

IFC JPT.GT. CNV-1)) GO TO 100 
JSC= JPT 
JPT= JPT+ 1 

STRIP GENERATION 

19 1=1 

20 T= IABSC VC I, JPT) - VC I, JSC) ) 

IF C T.GT.S ) GO TO 22 
1 = 1+1 

IF C I.LE.ND ) GO TO 20 
JPT = JPT + 1 
IFC JPT.LE.NV) GO TO 19 
22 JPT = JPT - 1 

COMPUTE MEAN OF STRIP 

NS= JPT—JSC+ 1 
DO 26 1=1, ND 
26 VSC I)=0. 

DO 30 J= JSC, JPT 
DO 28 1=1, ND 
28 VSC I)=VSC I)+VC I, J) 

30 CONTINUE 
T= FLOATC NS) 

DO 32 1=1, ND 
32 VSC I) =VS( I) /T 

PRIORITY SEARCH FOR NEAREST CLUSTER 

34 KMIN= 1 . E+10 
DO 40 1=1, NVM 
L=PLIST( I) 

D=0. 

IFC ID 1ST. EO. 2) GO TO 3601 
DO 36 J= 1 , ND 

36 D=D+ABSC VMCJ,L) -VSC J) ) 

GO TO 3602 

3601 CONTINUE 

DO 3603 J= 1 , ND 
3603 D=D+C VMCJ,L)-VSC J) )**2 
D=SQRTC D) 

3602 CONTINUE 
1112 FORMATC ) 

IFCD.GT.RMIN) GO TO 40 

RMXN=D 

jl=L 

IFC KMIN. LT.RP) GO TO 42 
40 CONTINUE 
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ORIGINAL PAGE IS 
OP POOR QUALITY 


TEST DISTANCE TO CLUSTER FOR ASSIGNMENT OF STRIP 

42 IF(RMIN.GT.R) GO TO 50 

FIRST PASS » ASSIGN TO J1 

IFC IPASS . EGL, 2) GO TO 200 

CALL MODIFY! VMC 1, JI) , VS, NVGC Jl) ,NS,ND) 

GO TO 12 

FIRST PASS - NEW GROUP, POSSIBLE ELIMINATION OF SMALLEST 

50 IFC IPASS. EQ. 2) GO TO 210 
NVM=NVM+1 
J1=NVM 

IFC NVM. LE . NVMMAR) GO TO 60 
NVM=NVMMAX 
NVGP= 10000 
DO 52 1-2, NVM 

IFCNVGC I) .GT.NVGP) GO TO 52 
IMIN= I 
NVGP=NVGC I) 

52 CONTINUE 

JPTT= JPTP+ JPT 

900 F0RMATC9H CLUSTER , 13, 8H WEIGHT , 14, 19H ELIMINATED, JPT - , I4,7H 
*VM = 13) 

NVGC 1 ) =NVGC 1) +NVGC IMIN) 

J 1= IMIN 

60 DO 62 1=1, ND 
62 VMC I, J1)=VSC I) 

NVGC J1)=NS 
NPC= NPT 
GO TO 12 

SHALL CLUSTER TESTS AND POSSIBLE ELIMINATIONS 

70 NEC=0 
NPC=NPT 
1=1 

72 1=1+1 

74 IFC I. GT. NVM) GO TO 14 

IFCNVGC I) .GT.NMIN) GO TO 72 
JPTT= JPTP+ JPT 

IFCIP.NE.0) WRITEC NODT, 900) I, NVGC I) , JPTT,NVM 
NVGC 1 ) =NV'GC 1 ) +NVGC I ) 

GALL PACK! VM, ND , NVM, I ) 

CALL PACKC NVG, 1 , NVM, I) 

NVM=NVM-1 
GO TO 74 

TESTS FOR MERGING AND POSSIBLE MERGING 

80 NPC=NPT 
NMC=© 
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G MINIMUM DISTANCE BETWEEN CLUSTERS CALCULATION 

G 

G 

82 12- NVM- 1 

IF(I2.LE.O) GO TO 16 
RMIN= l.E+10 
DO 88 1=1,12 
11 = 1+1 

DO 88 J= 1 1 , NVM 
D=0* 

IF( ID 1ST. EQ.. 2) GO TO 8400 
DO 84 L=1,ND 

84 D=D+ABS!VM!L, I)-VM!L,J>) 

GO TO 8401 

8400 CONTINUE 

DO 8402 L= 1 , ND 

8402 D=D+!VM!L, D-VMCL, J>)**2 
D= SORT! D) 

8401 CONTINUE 
1007 FORMAT! ) 

IF!D.GT.RMIN) GO TO 86 

RMIN=D 

Jl= I 

J2=J 

86 CONTINUE 
88 CONTINUE 
C 

C THRESHOLD TEST 

C 

IF(RHIN.GT.C) GO TO 16 
C 

C CLUSTER MERGING, J2 INTO J1 

C 

JPTT= JPTF+ JPT 

IFCIF.NE.0) NRITE(NOUT,902> J2,NVG! J2) , J1,NVG( Jl) ,NVH, JPTT 
902 FORMAT! 7H MERGER, 6H J2 = ,I2,7H NJ2 = ,14, 

#6H Jl = , 12, 7H NJ 1 = , 14 , 7H NVM - ,I2,7HJPT= , I4> 

CALL MOD IFYC VM! 1 , J 1 ) , VMC 1 , J2> , NVG! Jl ) , NVG! J2) , ND) 

CALL PACK! VM, ND , NVM, J2) 

CALL PACK! NVG, 1 , NVM, J2) 

NVM- NVM- 1 
GO TO 82 
C 

C PLIST UPDATE 

C 

90 NPC=0 

CALL UPPLT! PLIST, NVG, NVM) 

GO TO 18 
C 

C POSSIBLE SPECIAL CASE FOR LAST POINT 

C 

100 IF! JPT.NE.NV) RETURN 
DO 102 1=1, ND 
102 VS! I)=V( I,NV) 

NS= 1 
JSC= JPT 
GO TO 34 
C 

C SECOND PASS - ASSIGN TO Jl AND UPDATE STATISTICS 

C 

200 CALL THRDST! RMINM! J 1 ) , RMINV! J 1 ) , NVG! J 1 ) , RMIN , NS ) 

DO 202 J= JSC, JPT 

CALL SEQST!VMP! 1, Jl) ,VAR! 1, Jl) ,NVG< Jl) ,ND,V(1, J) ) 

TD IS ! J ) = RMIN 
202 IMG! J) = Jl 
GO TO 12 
C 

C SECOND PASS - ASSGN TO UNASSIGNED 

C 

210 DO 212 J= JSC, JPT 
TDIS!J)= 10000 
212 IMG! J> = 1 

NVG! 1 ) =NVG! D+NS 

GO TO 12 

END 
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C * 
C COMBINE * 
C jj> 
C THIS PROGRAM COMBINES SIGNATURES # 
C WITH .4 GIVEN SPECIFIER FROM TWO FILES * 
C * 
C SEQUENCE: COMBINE FILE1 FILE2 FILE3 * 
C * 


C 

C 

C 

c 

c 


CREATED AT GEORGIA TECH EES 
PROGRAMMER: NICHOLAS L. FAUST 


* 

;j: 

:f: 

* 

* 




DIMENSION A( 4, 20) ,B!4,4,20) f Nil (20) ,N21C20) ,N31(20> ,N41!20) 
DIMENSION IFK20) , IF2C20) , IF3I20) , C! 4, 20) , D! 4, 4 , 20) 
DIMENSION E( 4, 30) ,F!4,4,30) ,N12(20) ,N22!20) ,N32(20) ,N42C20) 
DIMENSION Nl( 30) ,N2!30) ,N3(30) ,N4!S0) 

DIMENSION SM! 4> , COVC 4, 4) , ISW(2) 

CALL OPEN! i , "COM. CM" , 1 , I ERR) 

CALL GQHARG! 1 , IF 1 , ISW, IERR) 

CALL COMARG! 1 , IF I , ISW, IERR) 

CALL COMARG! 1 , IF2, ISW, IERR) 

CALL COMARG! 1, IF3, ISW, IERR) 

NS I G= 20 


WRITE! 10, 103) 
ND=4 


IP= 10 

READ! ID ISW 

CALL F0PENC2, IF1, "BD 

CALL F0PENC3, IF2, "B") 

CALL FOPENC 4 , IF3 , "B " ) 

READ BINARY! 2) ITAP1 , ITAP2, ITAP3, IDAT1 , IDAT2, IDAT3 

IF( ISW. NE. DREAD BINARY! 3) ITAP1, ITAP2, ITAP3, IDAT1 , IDAT2, IDAT3 

WRITE! 10,101) 

READ! 11, 100)NA1,NA2 
WRITE! 10, 102) 

NUN1-2 

READ! 11,104) ITAP 1 , ITAP2 , ITAP3 , ID ATI , IDAT2 , IDAT3 
NSIG=5 

CALL GSIG1! A, B, NA1 , NA2,NUH1 , NSIG,NUN1 , N1 1 , N21 , N31 , N41) 

WRITE! 12, 105) 

NUN2=3 

NSIG=20 

IF! ISW.EQ.2) 

: CALL GS I G1 ! C , D , NA1 , NA2 , NUM2 , NS I G, NUN2 , N12 , N22 , N32 , N42 ) 

: WRITE! 12, 106) 

: . .FIN 
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do ck=i,numi> 

: N1CID=K11(K) 

: N2CK)=N21UD 

: N3CIO -N3KIC) 

*- N4(IO=N41(K) 

: DO <K1=1,ND> 

: : ECICX , K) “ACK1 , K) 

: : DO ( E2- 1 , ND) 

: : : FCK1\IC2,K)=B(K1 ,K2,K> 

: : : , . FIN 

: J..FIN 

: . .FIN 
IF( ISTf.Ea. X) 

: DO ( M= I , NUM2) 

: : ML=NUM1+M 

: *. NIC ML) =N 12CH) 

: : N2CML) =N22( M) 

: : N3(ML)=N32(M) 

: : N4( ML) =N42( ID 

: : DO (.HI- 1 , ND) 

: : : E( Ml , ML) =C(M1 , M) 

: : : BO ( 112= X,ND) 

: : : i F(m,M2 f ML>-D(Ml,M2,M) 

: : : :..FIN 

: : t . . FIN 

: : ..FIN 

: . .FIN 

• IFC ISW.EQu 1)KL=NUMI 
NP= 100 

WRITE BINARY! 4) ITAPI, ITAP2, I TAP 3 , I DAT I , IDAT2, IDAT3 
DO ( 1=1, ML) 

: WRITE B INARYC 4) NH D , N2( I) , N3( I ) , N4< I > 

: DO (11=1, ND) 

: : SK( I1)=E< II, I) 

: DO C 12= 1 , ND) 

: : : COV( 1 1 , 12) =F( I X , 12, I) 

f : S..FIM 

: r . .FIN 

: WRITE B INARYC 4) HP, SM, COV 

: WRITEC IP, X04)N1C I) ,N2C D , N3( I) , N4C I) 

: WRITEC IP) NP , SK, COV 

i . . FIN 

X@0 F0RHATC2A2) 

X©1 FORMAT! 2X, " INPUT NAX , NA2 ” ) 

102 F0RMATC2X, “INPUT TAPE, DATE IN 216 FORMAT") 

103 F0RMATC2X, "INPUT OF FILES") 

104 FORMATC 6A2) 

105 FORMATC //2X, "ABOVE ARE SIGNATURES FROM FILE! ") 

106 FORMATC /Y2X, "ABOVE ARE SIGNATURES FROM F ILE2 " ) 

STOP 

END 
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G * 

C C0PY5 * 

G ;■; 

C THIS ROUTINE ADDS 5 SIGNATURE , VERTICES , OR HIST * 

C FILES AND PUTS THEM IN A GOPYFILE * 

G * 

G SEQUENCE: C0PY5 OUT INI IN2 IN3 IN4 IN5 * 

C * 

G x 

C CREATED AT GEORGIA TECII EES * 

C * 

C PROGRAMMER: NICKOLAS L. FAUST * 

C * 

DIMENSION IOUT( 20) , INIC20) , IN2C20) , IN3(20) , IN4C20) , IN5C20) 

DIMENSION IC0UNT( 4 , 100) , 1(6) ,SM(4) ,C0V<4,4) , IVXC 101) 

DIMENSION IVY( 101) , ISWC2) 

GALL GPENC 1 , "COM. CM" , I , I ERR) 

COMARG( 1 , TOUT, ISW, I ERR) 

G0MARG( 1 , I OUT, ISW, I ERR) 

INI, ISW, IERR) 

IN2, ISW, IERR) 

IN3, ISW, IERR) 

IN4, ISW, IERR) 

INS , ISW, IERR) 

ORIGINAL PAGE IS 
OF POOR QUALITY 


CALI. 

GALL 

GALL 

GALL 

GALL 

GALL 

GALL 


COMARG( 1 
COMARG( 1 
CGMARG( 1 
GOKARGC 1 , 

CQMARGC 1 
IPR= 10 

WRITE( IPR, 40 1 ) 

READt 11) ICON 
WRITE( IPR, 402) 

READ( 1 1)M 
L=M+2 

GALL FOPEN( 2, IOUT, "B") 

CALL F0PENC3, INI, "B") 

CALL FOPENC4, IN2, “B") 

CALL FOPENC 5 , IN3, "B") 

GALL FOPENC 6, IN4, "B“) 

CALL F0PENC7, IN5, "B") 

DO C K=3 , L) READ BINARYC K) I 
WRITE BINARYC 2) I 
TYPE " L = ",L 
DO CICL=3,L) 

FOREVER 

READ B I NARYC KL , END= 120 ) NS 1 , NS2 , NS3 , NS4 
WRITEC IPR, 403 ) NS 1 , NS2 , NS3 , NS4 
WRITE B INARYC 2) NS 1 , NS2 , NS3, NS4 
IFC ICON.EQ. 1) 

READ BINARYC EL) NP , SM, GOV 
TYPE 11 NP - ” , NP 

mi TE B I NARYC 2 ) NP , SM, COV 
FIN 

IFC ICON.EQ. 2) 

READ BINARYC KL) IC, IYMIN, IYMAX 

DO CK1= 1,10 READ BINARYC KL) K1 1 , IVXCK1) , IVYCIC1) 
WRITE BINARYC 2) K, IYI>IIN, IYMAX 
DO ( K1 = 1 , IQ WRITE B INARYC 2 ) K1 , IVXC IC1 ) , I VYC IC1 ) 
.FIN 

IFC ICON.EQ. 3) 

READ BINARYC KL) ICOUNT 
WRITE BINARYC 2) ICOUNT 
.FIN 


2 - VER , 3 - HIS") 


120 

: : . . FIN 

: CONTINUE 


401 

: . .FIN 
FORMAT! 2X, " 

INPUT SWITCH - 1 -SIG 

402 

FORMAT! 2X, 11 

INPUT * OF FILES '*) 

403 

FORMATC 2X, 4A2) 


STOP 

END 
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p *>* *i» *> «-V •> »>■>* tV't* kl+hl# \l< tl> il* ci» ti» *i*»U «'»«'J«^ »J-« «J-« lj^ vV- »V »J»«i» *4> ^1» «i» «fj «la ».£» *J* «.t« *L» 

C i*i 

G COUNTY * 

G # 

G DRAWS A POLYGON WHEN GIVEN A STRING OF COORDINATES. * 

C * 

C 5): 

C CREATED AT GEORGIA TECH EES * 

G * 

G PROGRAMMER : MICHAEL D, FURMAN * 

C * 


INTEGER ICC2, 256), IXYC2, 2) 

ACCEPT "# OF POINTS, STRING OF POINTS \ NP , < C ICC I , J) , 1= 1 , 2) , J= I , NP) 
DO C 1=1,2) 


: DO CJ=I,2) IXYC I, J)=ICC I, 1) 

:..FIN 


DO ( J= X , NP) 

: DO ( IC= 1 , 2) 

: : IF (C ICCK.J) .LT. IXYCK, 1)) .AND. I ICCK, J) .GE.0)) IXYC K, 1) = ICC K, J) 

: : IF ( ICCK, J) . GT. IXYCK.2) ) IKYC IC, 2) = ICC K, J> 

: : . . FIN 

: . .FIN 

FAR= ( IXYC 1,2)“ IXYC 1 , 1) ) /5 1 1 . 0 
FAY=C IXY<2,2)-IXYC2, 1) )/511.0 
IfflEN CFAX.GT. FAY) FACT= FAX 
ELSE FACT=FAY 
DO ( 1=1, NP) 

: DO CJ=1,2) ICCJ, I) = C ICCJ, D-IXYCJ, 1))/FACT 

: . . FIN 


NP1=NP- 1 

DO ( 1= 1 , NP1) CALL VECTORC0, ICC 1,1), IGC2, I) , ICC 1 , 1+1) , ICC2, 1+1) ,0,200) 
CALL VECTORCO, ICC 1, I) , ICC2, I) , ICC 1, 1) , ICC2, 1) ,0,200) 

STOP 

END 
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• <;» 
fc 

if" 




G X 

C CRDEM4 * 

G Hi 

C CREATE DEMO FROM DISPLAY IMAGE INCLUDING * 

G FUNCTION MEMORY, COLOR MEMORY, AND GRAPHICS * 

c ^ 

C SEQUENCE: CRDEM4 MTU: F * 

C * 


C Hi 

C CREATED AT GEORGIA TECH EES X 
G :}: 
C PROGRAMED: FRED L. THOMPSON X 
C * 




INTEGER ISW! 2) , INPUT! 512) , IARAY! 64) , IRRAY( 0:255) , IFLD( 10) 
CALL OPEN C 1 , "COM. CM" , I , IE) 

CALL COMARG C I , IFLD, ISW, IE) 

CALL COMARG (1, IFLD, ISW, IE) 

CALL MTOPD (3, IFLD, 0, IE) 

DO ! 1=1,5 12> INPUT! I) -0 
CALL RCM (0, INPUT! 3)) 

ACCEPT "512 OR 256 " , ISZ 

ACCEPT "TYPE 1 TO RECORD GRAPHICS 'MANS 

TYPE "TYPE A THIRTY CHARACTER DESCRIPTION » 

TYPE " 


READ !11, 100) ! INPUT! I) , 1-68,98) 

100 FORMAT !3©A1) 

INPUT! 1)=ISZ 
INPUT! 2) -IANS 

CALL MTDIO (3,50000K+ISZ, INPUT, 1ST, IE, ICNT) 

DO ( 1=0,2) 

: CALL RFUM ! I , INPUT) 

: CALL MTDIO (3, 500GOK+ISZ, INPUT, 1ST, IE, ICNT) 

•* . .FIN 


IEND= ISZ-1 


ISZl=ISZ/2 
DO C IY=0, I END) 

: CALL IKRD (0, I Y, INPUT, ISZ 1) 

: CALL MTDIO ! 3 , 50000K+ ISZ, INPUT, 1ST, IE, ICNT) 

: IF ! ISZ, EQ. 256) 

: : DO C 1= 1 , 2) 

: : : CALL IMRD ! I , IY, INPUT, ISZ1) 

: : : CAJX MTDIO ! 3, 50O00K+ ISZ, INPUT, 1ST, IE, ICNT) 

: : J..FIN 

: : IF C IANS.EQ. I) 

: : : DO ! 1=0,3) 

: : : : CALL GRD C I , IY, INPUT, 16) 

: : : : CALL MTDIO ! 3, 50000K4 ISZ, INPUT, 1ST, IE, ICNT) 

: : : :..FIN 

: : :..FIN 

: : . .FIN 


: . .FIN 

DO ( 1=1,2) CALL MTDIO C3.6000OK, INPUT, 1ST, IE, ICNT) 

STOP 

END 
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Q 

G CYTAPE * 

C # 

C MAKES FILE INTO A CYBER 74 COMPATABLE FORMAT * 

G * 

G SEQUENCE: CYTAPE FILE FILE * 

G * 

C * 

C CREATED AT GEORGIA TECH EES * 

C x 

G PROGRAMMER: MICHAEL D. FURMAN * 

G X 

DIMENSION SC30) , IFILEC 10) ,NSC4) , IFILE2( 10) , ITAPC3) , IDATC3) , ISWC2) 
CALL OPENC 1 , "COM. CM" , 1 , I ERR) 

CALL COMARGC 1 , IFILE, ISW, I ERR) 

CALL COMARGC 1, IFILE, ISW, I ERR) 

CALL COMARGC 1, IF ILE2, ISW, I ERR) 


CALL FOPENC 2, IFILE, "B") 
CALL FQPENC3, IFILE2, "B") 


READ B INARYC 2 , END=203 ) ITAP , IDAT 
FOREVER 

DO C 1“ 1 ,30) S( I)=0. 

READ B INARYC 2 , END- 205 ) NS 
READ B INARYC 2 , END= 205 ) NP, < S< I) , 1= 1 , 20) 
WRITE B INARYC 3) NS 
WRITEC3) 

WRITEC3, 100) NP, CSC I) , 1= 1,20) 

..FIN 
STOP 

205 STOP FINISHED 

1©0 FORMAT ( IK, I4,5C4F12.8,/) ) 

END 
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G # 

C DISDEM4 * 

G is*? 

C DISPLAYS DEMO TAPES TO COMTAL COLOR DISPLAY * 

G WITH GRAPHICS COLOR MEMORY AND FUNCTION MEMORY * 

C * 

C SEQUENCE: DISDEM4 MTU:F * 

C * 

G # 

G CREATED AT GEORGIA TECH EES * 

C * 

G PROGRAMMER: FRED L. THOMPSON # 

C Hi 


INTEGER ISW(2) , INPUT! 5 12) , IFLD! 10) 
CALL OPEN ( 1, “COM. CM", 1, IE) 

GALL COMARG ( 1 , IFLD, ISW, IE) 

GALL COMARG ( 1 , IFLD, ISW, IE) 

GALL KTOFD (3, IFLD, 0, IE) 

FOREVER 


: GALL MTDIO ( 3, 0, INPUT, 1ST, IE, ICNT) 

: ISZ= INPUT! 1) 

: I ANS= INPUT! 2) 

: IVRITE < 10, 10G) ! INPUT! I) , 1=68,98) 

100 : FORMAT t 1 X, 3 0 A 1 ) 

: IF ( ISZ.EQ.512) PAUSE SET BIG FOR 512 HIT RETURN 

: IF ! ISZ.EQ. 256) PAUSE SET BIG FOR 256 HIT RETURN 

: CALL WCM !0, INPUT! 3)) 

: DO ! 1=0,2) 

: : CALL MTDIO ! 3, 0, INPUT, 1ST, IE, ICNT) 

: : CALL NFUK ! I , INPUT) 

: : . .FIN 

: IEND= ISZ- 1 

: ISZ1= ISZ/2 

: DO ! I Y=0, I END) 

: : CALL MIDIO (3,0, INPUT, 1ST, IE, ICNT) 

: : CALL IKWR (0, IY, INPUT, ISZ1) 

: : IF (ISZ.EQ.256) 

: : : DO (1=1,2) 

: : : : CALL MTDIO (3,0, INPUT, 1ST, IE, ICNT) 

: : : : CALL IMWR C I , IY, INPUT, ISZl) 

: ’• : :..FrN 

: : : IF ! IANS. EQ. 1) 

s : : : DO ( 1=0,3) 

: : : : : CALL MTDIO (3, 0, INPUT, 1ST, IE, ICNT) 

: : : : : CALL GWR ( I, IY, INPUT, 16) 

: : : : :..FIN 

: : : :..FIN 

: : : , .FIN 

: : . .FIN 

: CALL MTDIO (3,0, INPUT, 1ST, IE, ICNT) 

: PAUSE HIT RETURN FOR NEXT FILE CTRL A TO END 

: . .FIN 

STOP 

END 
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C sfc 

C DRAWER * 

c m 

C USED WITH EL IP TO DRAW 2-D ELLIPSES OF SIGNATURES * 

C ON THE VERSATEC « 

C * 

C y 

C GREATER AT GEORGIA TECH EES * 

C ;j: 

G PROGRAMMER: G. DAVID GENTRY * 

C if; 

SUBROUTINE DRAWERC SM, E, COV, AX1 , AX2 , R, ID) 

DIMENSION RC2,2) ,SMC2> ,EC2,4) ,GOV(2,2) 

IF( ID. GE. I) GO TO 9 
CALL DRAW . 5 , . 5 , 1 , 100 1 ) 

CALL K0DEC7, 8. ,8. , 9999. ) 

GALL MQDEC 8, ISO. , 10. ,0.) 

GALL M0DEC9, 180. , 10. ,0.) 

GALL AXES ( 7 . 0 , " CHANNEL 11 , 7 . 0 , " CHANNEL " ) 

GALL NQTE(-.5,5. , AX2, 1000) 

GALL NOTEC 5, , ~.5,AX1, 1000) 

9 10= 10 

IF( ID.EQ.0) IDC=43 
IFC ID.EO. 1) IDG=31 
IFC ID. EQl. 2) IDC=42 
THETA=ATAN2CRC2, 1) , RC 1, 1)) 

IFC COV( 1 , 2) .LT. 0. )THETA=ATAN2CRC2,2) ,R( 1,2) > 

IFC COV( 1,2) .LT. 0. )WRITE( 10, 12) 

12 FORMAT! // 11 CGVCX,Y) IS NEGATIVE") 

THET=57 . 29573*THETA 
mi TEC 10, 1 1) THET 

1 1 FORI'IATC // 11 THE ROTAT I ON ANGLE IS" ,F6.2, ” DEGREES V/ ) 

DO C 1=1800,2600, 10) 

J=C C 1“ 1800) / 10) + 1 
X1=FL0ATC I)/ 10. 

X=X1-SMC 1) 

IFC COVC 1,2) .LT. 0. ) GO TO 6 

A= EC 2 , 1 ) * ( S INC THETA) ) **2+E< 1 , 1)*C COS ( THETA) ) **2 
B=2.*X*CE(2, 1) *CS INC THETA) )*< GOSC THETA) ) -EC 1, 1) *CS INC THETA) ) 

Is *( GOSC THETA) ) ) 

C=X**2*C EC 2 , 1 ) *( COSt THETA) ) **2+E( 1 , 1) *C S INC THETA) ) **2) 

Is -3.#EC 1, 1)*EC2, 1) 

IFC E*:[j2-4 . *A#C) 5,7,7 

A=EC 1, 1)>MSINCTHETA))**2+E<2, 1) *( GOSC THETA) ) **2 
B=2.#X*CEC 1, 1 ) *C S INC THETA) ) *C GOSC THETA) ) -EC 2 , 1) S INC THETA) ) 

1: *C GOSC THETA) ) ) 

C=X**2*CEC 1, 1 ) #< GOSC THETA) ) **2+EC 2 , 1) *C S INC THETA) ) **2) 

1: ~3 . &EC 1 , 1) JtfEC 2, 1) 

IFC B**2-4 . *A*C) 5,7,7 

Y1=SMC 2) +C-B+SQRTC . *A*C) ) /< 2 . «A) 

Y2= SMC 2) + ( -B-SQRTC B**2-4 . #A*C) )/<2.*A) 

GALL NOTEC XI , Y1 , IDG, - 1) 

CALL NOTEC XI, Y2, IDG,- 1) 

CONTINUE 
• FIN 

SPL0T1=SMC 1) 

SPL0T2=SM( 2) 

GALL NOTEC SPL0T1 , SPL0T2 , IDG , - 1 ) 

RETURN 
END 
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«p «t* -,l> si> »U Ui *1* kl> \1+ »l< *U ii> 0>i]» i'* >l< Ui tl> *V# iU * 1* mAj *L» «.]> >U * l> lU «,*» ilj *U <i> «J-* lb tif «iU *A* »l* v^» l> vU \l* O* ki» kU »V» tl» »|i»b >J> (J/ %t* o* *1* *A* ^i> *> 

* 

DTAPE * 

if. 

GENERAL TAPE DUMPING ROUTINE FOR READING AMD * 

DEBUGGING FOREIGN TAPES * 

if 

SEQUENCE: DTAPE MTU: F if 

:f 

^^^^.^^JfJfJf^^^Jj'.ifjf^^Htitiifil^iifili^^^if^ififilJ^ifiiiifif^^^^^ififififif^if^ififif^il'if^ifififii^.ififiifiliii^^ifif^if 

if 

CREATED AT GEORGIA TECH EES if 

* 

PROGRAMMER: MICHAEL D. FURMAN if 

if 

Cif^:ij:^ifi^:ifif^^^^if:f;f:f^:f^^^i^iJ;i!iif^:^:;l:if^:*;fif:f:fif“!::fiJi;f^:il:i|iil:iiiififii;.‘f.'i:^:^;:fil"fil:^il;if;f^:ifi):if:fif^:ifififif 

INTEGER IDATA( 41 10) , I0UT( 8220) , IFLD( 10) , ISW! 2) 

CALL 0PEN( 1 , "COM. CM" * 1 , I ERR) 

CALL COMARG! 1, IFLD, ISW, I ERR) 

CALL COIIARG! 1, IFLD, ISW, I ERR) 

100 FORMAT! IX, 2016) 

101 FORMAT! IS, 25 14) 

102 FORMAT! IS, 100A1) 

103 FORMAT! IS, 10! 018, 1IIK) ) 

1 04 FORMAT! IS, 20 ! 0 1 5 , 1HK) ) 

ACCEPT "WORDS? Cl), OR BYTES? CO) OUTPUT “.IP 
ACCEPT "BASE: TEN= 1 , EIGHT= 0 ",IN 
ACCEPT "SKIP HOW MANY RECORDS? \ISKP 
ACCEPT "DUMP HOW MANY BLOCKS OF DATA " , I BLOK 
IF C! IP.EQ.O) .AND.! IN.EQ. 1)) 

TYPE "TYPE: NUMERIC= 0 ; ASCII- 1 ; EBCDIC= 2 tt 
READ! 11) IQ 
. .FIN 

CALL MTOPD! 3 , IFLD , 0 , I ERR) 

IF CISKP.GT.0) 

: DO ! 1= 1 , ISKP) CALL MTDIO!3,0, IDATA, IS, IE, IC) 

: , .FIN 

DO C Kl= 1 , IBLOIO 

DO ! I— 1 ,4110) IDATA! I) -0 
DO C 1=1,3220) I OUT! I ) =0 
CALL HTDIO! 3,0, IDATA, IS, IE, IC) 

TYPE "RECORDS PER BLOCIC= ” , IC 
IF ! IP.EQ.O) 

CALL UPAC3! IDATA, IOUT, IC) 

IF CIQ.EQ.2) CALL EBCDIC! IOUT, IOUT, ID) 

IC= IC:f2 
. .FIN 

WHEN CIQ.NE.0) WRITE! 12, 102) ! IOUT! I) , 1= 1, IC) 

ELSE 

IF C! IN.EQ, 1) .AND.! IP. EQ. 1)> WRITE! 12, 100) C IDATA! I) , 1= 1 , IC) 

IF !! IP.EQ. 1) .AND.! IN.EQ.0)) WRITE! 12, 103M IDATA! I) , 1= 1 , IC) 

IF C ! IP.EQ.O) .AND. ! IN.EQ. 1) ) WRITE! 12 , 10 1) ( IOUT! I> , 1= 1 , IC) 

IF !! IP.EQ.O) .AND.! IN.EQ.0)) WRITE! 12, 104) ! IOUT! I) , 1= 1 , IC) 

. .FIN 
WHITE! 12) 

. .FIN 
STOP 
END 
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V 


G 

C EBCDIC (SUBROUTINE) * 
C * 
C CONVERTS EBCDIC CODE INTO ASCII # 
C X 
C SEQUENCE: EBCDIC ( INPUT ARRAY, OUTPUT ARRAY, ARRAY DIMENSIONS) * 
C * 


C 

c 

c 

c 

c 


CREATED AT GEORGIA TECH EES 
PROGRAMMER: MICHAEL D. FURMAN 


* 

X 

X 

X 

X 


CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


SUBROUTINE EBCDIC( INPUT, IOUT, IPR) 

INTEGER INPUT( IPR) , IOUTC IPR) 

COMKON/DUM1/ ICNTABC 263) 

DATA ICNTAB/ 1H? , IH? , in? , 1H? , 1H? , 1H? , 1H? , IH? , IH? , 1 H? , I H? , 
1 1H? , IH?, IH?, IH?, IH?, IH?, IH?, III?, IH?, IH?, IH?, IH?, IH?, IH?, 
1 IH? , IH?, IH?, IH?, IH?, IH?, III?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, 
I IH? , IH? , IH?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, 
2IH?, IH?, IH?, IH?, IH?, III?, IH?, IH?, IH?, IH?, IH?, IH , IH?, IH?, 
21H?, IH?, IH?, IH?, IH?, IH?, IH? , III?, IH. , IH?, 1H( , IH+, IH?, IH?, 
21H? , IH?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, IHS, IE*, IH) , IH?, 
31H? , 1H-, IH/, IH?, IH?, III?, IH?, IH?, IH?, IH?, IH?, IH?, IH, , IH?, 
31H? , IH?, IH?, IH?, IH?, HI?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, 
3IH?, IH?, IH’ , 1H=, IH?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, 
41H?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, 
41H?, IH?, IH? , IH?, IH?, III?, IH?* IH?, IH?, IH?, IH?, IH?, IH?, IH?, 
41H?, IH?, IH?, IH?, IH?, in?, IH?, III?, IH?, IH?, IH?, IH?, IH?, IH?, 
S1H?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, IH?, 
51HA, IHB, 1HC, I HD, 1HE, 1HF, 1HG, 1HH, 1HI, IH?, IH?, IH?, IH?, IH?, 
51H?, IH?, IHJ, 1HK, IHL, 1HM, 1HN, 1HO, IHP, I HQ, 1HR, in?, IH?, IH?, 
6IH?, IH?, IH?, IH?, IH?, I IIS, 1HT, 1HU, 1HV, 1HW, IHX, 1HY, IHZ, IH?, 
6 IE?, IH?, IH?, IH?, IH?, 1110, 1H1, 1H2, 1H3, 1H4, 1H5, 1H6, IH7, 1H8, 
61E9, IH?, IH?, IH?, IH?, IH?/ 


DO (1=1, IPR) 

: Kl= INPUTC I) + l 

: UHEN (K1.LE.256) 

: : IOUT( I) = IGNTAB( Kl) 

: : . .FIN 

: ELSE IOUT( I) = ICNTABC I) 

: . .FIN 

END 


70 






G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

C 

G 

C 

c 

G 

G 


ED 1ST (SUBROUTINE) 

PURPOSE: 

COMPUTES THE DISTANCE BETWEEN TWO VECTORS 

DESCRIPTION OP PARAMETERS 

CALLING SEQUENCE 
INPUT 

VI, V2 - TWO VECTORS 

ND - DIMENSION OF THE VECTORS 

OUTPUT 

GALLING SEQUENCE 

D “ CONTAINS DISTANCE BETWEEN VI AND V2 


* 

# 

X 

* 

* 

•J* 

u* 

* 

Hi 

* 

* 

* 

* 

* 

* 

* 

* 


G ^ 
C CREATED AT NASA/ JSC (ASTEP) * 
C % 




SUBROUTINE ED 1ST ( VI, V2, ND, D ) 
DIMENSION Vl(ND) ,V2(ND) 

D = 0. 

DO 10 I-l.ND 

10 D - D + (VIC I)-V2( I) )#*2 
D = SQRTCD) 

RETURN 

END 


71 




C 

G 

G 

G 

C 

C 


EIGEN ( SUBROUTINE) 

COMPUTES EIGENVALUES AND EIGENVECTORS OF REAL 
SYMMETRICAL MATRIX 


# 

* 

* 

* 

:j: 

* 


202 


10 


15 

20 


J 


G * 

G CREATED AT NAS A/ JSC CASTEP) # 

C * 

SUBROUTINE E I GEN ( AA, N , HV , A, E » R) 

DIMENSION AA( N, N) ,RCN,N) ,E(N) ,A(N,N> 

FORMATC // " SUBROUTINE EIGEN ERROR RETURN-FINAL NORM HAS NOT 
1BEEN REACHED AFTER 100 ITERATIONS ”/) 

DO (1=1,10 
DO C J= 1 , N) 

: AC I, J)=AA( I, J) 

: . .FIN 
.FIN 

RANGE=1.E“6 
IF(MV-l) 10,25, 10 
DO C 1= 1 , N) 

DO ( J“ 1,N) 

IF( I ,EQ. J) GO TO 15 
R( I, J)=0. 

GO TO 20 
R( I, J) = l. 

CONTINUE 
. .FIN 


• ..FIN 

25 ANORM=0. 

DO C 1=1, N) 

: DO ( J=1,N) 

: : IF( I . EG. J) GO TO 35 

: : ANORM= ANORM+AC I , J) *AC I , J) 

35 : : CONTINUE 

: : . .FIN 

: . .FIN 

IFCANORM) 165, 165,40 
40 ANORM=SQRT( ANORM) 

ANRMX= ANORM*RANGE 

ICNT=0 

IND=0 

THR= ANORM 

45 THR= THR/FLOATC N ) 

IGNT= ICNT+1 

IFC IGNT .EQ.' 100) GO TO 200 
50 MQ=2 

55 MP= 1 

X= * 5*C AC HP , MP) - AC MO, MQ) ) 

62 IFC ABSC ACMP,MQ) )-T5HR) 138,65,65 

65 IND= 1 

X= . 5*C AC MP , KP) -AC MG, MQ) ) 

Y= -AC MP , MG) /SGRTC AC HP , MG) *AC MP , MG) +X*X) 
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70 

75 


120 


125 


isa 

140 

145 

160 

165 

190 


20 © 

210 


ORIGINAL PAGE IS 
OF POOR QUALITY 


IF( X) 70,75 , 75 
Y--Y 

S I NX- Y/SQRTC 2 . *( 1 . 4* C SQRTC 1 . -Y*Y) > ) ) 

SINX2=SINX*SINX 

COSX= SQRTC 1 .-SINX2) 

C0SX2= COSX*COSX 
SINCS=SINX*COSX 
XX= AC HP, TIP) 

YY=A( MQ, MQ) 

ZZ=ACMP,KQ) 

DO C 1= 1,N) 

: X= A( I , HP ) *COSX-AC I , MQ) *S INX 

: AC I , MQ) = AC I , MP ) *S INX+ AC I , MQ) *COSX 

: AC I , MP) =X 

: IFC HV— 1) 120 , 125 , 120 

: X-RC I , MP) ^COSX-RC I , MQ) *S INX 

: RC I , MQ) = RC I , UP) *S IHX+RC I , MQ) *COSX 

: RC I , MP) =X 

: CONTINUE 

: . .FIN 

X=2.*ZZ*SINCS 
Y= C XX*CGSX2) 4*C YY*S INX2) -X 
X= C XX*S I NX2 ) + C YY*C0SX2 ) 4*X 
AC MP , HP) = Y 
AC MQ, MQ) =X 

AC MP , MQ) = C XX-YY) *S INCS4-ZZ*C G0SX2-S INX2) 
ACMQ,MP)=0. 

DO C I- 1,N) 

: AC MP , I ) = AC I , TIP) 

: AC MQ, I) =AC I ,MQ) 

:..FIN 

IFCMP ,NE. C MQ- 1 ) ) GO TO 140 

IFCMQ ,NE. N) GO TO 145 

IFCIND . NE. l)GO TO 160 

IND=0 

GO TO 50 

MP=MP4*1 

GO TO 62 

MQ=KQ4-1 

GO TO 55 

IFCTHR — ANRMX) 165, 165,45 

CONTINUE 

DO C 1= 1 , N) 

: DO CJ=1,N> 

: : IFC I .NE. J) GO TO 190 

: : EC I)=AC I, J) 

i : CONTINUE 

: : ..FIN 

: . .FIN 
GO TO 21© 

NRITEC 10,202) 

RETURN 

END 
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ft *t* 

G EL IP * 

C * 

C DRAWS 2-D ELLIPSES OF COVARIANCE MATRICIES * 

G 

C SEQUENCE: EL IP REQFIL SIGFILC # 

C a: 

G * 

C CREATED AT GEORGIA TECH EES * 

C '!' 

G PROGRAMMER: G. DAVID GENTRY * 

G * 

DIMENSION REGC30) , SIGFILC 30) ,EC2,4) , R( 2, 2) , DC 2, 2) . IORDERC20) 
DIMENSION G0VINC4,4,20> ,SMINC4,20) , NIC 20) ,N2(20) ,N3(20) ,N4C20) 
DIMENSION C0V1C2,2) ,SH1<2> , ISWC2) 

GALL OPENC 1 , "COM. CM" , 1 , I ERR) 

CALL COMARGC 1,REG, ISW, IERR) 

GALL COMARGC 1 , REQ, ISW, IERR) 

CALL COMARGC 1,SIGFIL, ISW, IERR) 

GALL FOPENC 2 , REQ, 11 B " ) 

10=10 

ACCEPT "INPUT NUMBER OF ELIFSES PER PLOT n , NUMS 
I IN= 11 


ND=2 

NSIG=20 

CALL MODEC 1,1. ,.5,-1.) 

MS I G= NS I G/NUMS 
DOC JV= 1 , MSIG) 

: NSIG=NUMS 

: KQ=0 

: GALL GSIGC SHIN, COVIN, NSIG,2, 3, SIGFIL, I0RBER, N1 , 1?2, N3, N4) 

: DO C JJ= 1 , 3) 

: : K= JJ+1 

: : DO CEK=K,4) 

: : : AX1=JJ 

: : : AX2=KIC 

: : : ID=0 

: : : DO CII=1,NSIG) 

: : : : SI-I1C 1) =255 . -SMINC JJ, 1 1) 

: : : : SMI C 2) =255. -SMINC KK, 1 1) 

: : : : COV1C 1, 1) =COVINC JJ, JJ, II) 

: : : : C0V1C 1 , 2) =COVINC JJ, KK, 1 1) 

: : : : COV1C2, 1) =COVINCKK, JJ, II) 

: : : : C0V1C2.2) =GOVINC KK, IOC, II) 

: : : : CALL FACANLCCQV1,SH1,ND,E,R,D) 

: : : : CALL DR&WERCSH1 ,E,C0V1 , AX1, AX2,R, ID) 

: : : : ID=ID+1 

: : : :..FIN 

: : : KO=KO+I 

: : : IFCICO.EQ. 1) 

: : : : CALL DRAWCO. ,0. , 1,0000) 

: : : : CALL DRAWC-.5, 8.5, 1 , 1001) 

: : : J..FIN 

: : : IFCK0.EQ.2) 

: : : : KO=0 

: : : : CALL DRAWCO . , 0. , 1 , 0009) 

: : : : CALL DRAWCO. ,-9.5, 1, 1001) 

.* : : : CALL DRAWC 0 . , 0 . , 1 , 9000) 

: : : : . .FIN 

: : :..FIN 

: : . .FIN 

i . .FIN 

CALL DRAWC 0,0, 0,9999) 

STOP 

END 
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G ' * 

G FACANL (SUBROUTINE) * 

C * 

G USED WITH EL IP TO DO FACTOR ANALYSIS OF COV. MATRIX * 

C * 

C * 

C CREATED BY NASA JSC * 

C * 

SUBROUTINE FACANL( COV, SM, ND , E, R, D) 

D I MENS ION COV( ND , ND) , SMC ND) , EC ND , 4) , RC ND , ND) , DC ND , ND) 

MV"0 


CALL EIGENCCOV,ND,MV,D,E,R) 

ND1=ND-1 

IFCIID1 .Ed. O) RETURN 
DO C J=1 , ND1) 

: EMAX=0. 

: DO C I=J,ND) 

: : IFC EC 1 , 1) .LT. EMAX) GO TO 50 

: : EMAX=EC 1,1) 

: : rs=I 

50 : t CONTINUE 

: :..FIN 

: DUM=ECJ,1) 

: EC J, 1) =E( IS, 1) 

: EC IS, 1) =DUM 

: DO C 1= t , ND) 

: : DUH=R( I , J) 

: : RC I, J)=R( I, IS) 

i : RC I , IS) =DUM 

: J..F1N 

: . . FIN 
DUM= 0 . 

DO C 1= 1 ,ND) 

: DUM-DUH+EC 1,1) 

: . .FIN 

DUM=1./DUM 

EMAX=0. 

DO C 1= 1,ND) 

: EC I,2)=DUM*ECI, 1) 

: EC X,3>=EMAX+E( 1,2) 

i EM_4X=E( 1,3) 

: . .FIN 
DUM-0. 

DO (1=1, ND) 

: BUM- DUM+SMC I ) **2 

: . .FIN 

DUM= SQRTC DUM) 

DO C J= 1 , ND) 

: EC J,4)=0. 

: DO C 1= 1, ND) 

: : EC J,4)=ECJ,4)+SM< I)*R(I, J) 

: : . .FIN 

: EC J,4)=ECJ,4)/DUM 

: E( J, 4) =57. 29578&AC0SCEC J, 4) > 

: . .FIN 
RETURN 
END 


■ 5 * ~ 


C * 

G FILTER! * 

c ^ 

C PERFORMS A HIGH-PASS FILTER OR DISPLAY * 

C * 


G 'i' 

C CREATED AT GEORGIA TECH EES * 

G ^ 

G PROGRAMMER: MICHAEL D. FURMAN" * 

G * 

G »i^»fj »■> *.!> * t* *J» «1* *}-* kl- 4 *T**J#vf**L «.J»*t*«J>*L' vL»*]*«4* *1»»1**'*%;* si* »]* *1**1.# ****£.» 

INTEGER INC512,3) , IMAG1C256) , IKAG2C512) , IMAG3C256) 

INTEGER IN2C 512, 3) , IN3C 512,3) 

ACCEPT "512 OR 256 MODE " , IMODE 
IMOM1- IMODE- 1 
IMOM2” IMODE-2 
IANS=0 


IF < IMODE. EO. 256) 

: ACCEPT "ALL 3 IMAGES AT ONCE? ( 1=YES) ‘MARS 

: IF (IAMS.NE.I) ACCEPT “ IMAGE NUMBERC 0-2) MM2 

:..FIN 
DO C IM=0.2) 

i IF C( IANS, ME. 1) .AND. ( IMODE. HE. 512)) IM= IM2 
: 113=2 

: M4=3 

: DO ( 1=0, 1) 

i : GALL IMREADC IM, I, IKAG1,256) 

: : GALL UPAC8C IMAG1 , INC 1 , C 1+ 1) ) , 256) 

: i . .FIR 

i DO <M5=2, IM0M2) 

: : M5M1=M5-1 

: : CALL IMRE ADC III, M5 , IMAG1 , 256 ) 

: : CALL UPAC3C IMAG1 , INC 1 , M4) , 256) 

: t DO C 1=1, MODE) IMAG2CI)=0 

: : DO C Jl=l,3) 

: : : DO CK=2, IMOM1) INAG2CK) = EMAG2CK) + INC K- 1 , JD + INCK, JD + INCK+l 

: : :..FIH 

: : DO C 1=1, IMODE) 

: : : IMAG2C I) = C INC I, M3) *2) -C C IMAG2C I) /9 . 0) +0.5) 

: : : IF C IHAG2C I) .LT.0) IMAG2C I)=0 

: : : IF ( IMAG2C I) . GT. 255) IMAG2CI)=255 

: : : . .FIN 

: : CALL PAC8C IMAG2 , 1 MAGS ,512) 

: : GALL IMWRITEC III, M5M1 , IMAG3 , 256 ) 

: : FHEN CM4.EGt.3) M4= 1 

: : ELSE M4=M4+1 

: : NIIEN CM3.E0.3) H3= 1 

s : ELSE M3=H3+1 

: : . .FIN 

: IF CC IANS.NE. 1) .OR.C IMODE. ECL 5 12) > STOP 

: . .FIN 
GALL BACK 
END 
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G * 

C GETPOLY * 

c ^ 

C THIS PROGRAM SUMS AREAS OP ALL MAXIMUM LIKLIHOOD GLASS * 

C WITHIN A SPECIFIED POLYCOM * 

G POLYGON DATA IS THEN DUMPED TO TAPE * 

C ?ic 

C SEQUENCE: GETPOLY INPUT MTCRF MTO:F * 

C & 

C 

C CREATED AT GEORGIA TECH EES * 

C * 

C PROGRAMMER: NICICOLAS L. FAUST * 

G * 

DIMENSION COUNTC 60) , JBUF( 3300) , IVXC 201) , IVYC 201) 

DIMENSION LAC 233) ,KEUF( 3300) 

DIMENSION IF ILC 20) , ITAP< 20) , JTAPC 20) , ISWC2) 

I IN= 1 1 
110 = 10 


CALL OPENC 1 , "COM. CM" , 1 , I ERR) 

CALL CQMAR&C 1 1 IFIL, ISW, I ERR) 

CALL COMARGC 1, IFIL, ISW, I ERR) 

CALL COMARGC 1, ITAF, ISW, I ERR) 

CALL COMARGC 1, JTAP, ISW, I ERR) 

CALL FOPENC I IN, IFIL, “B") 

CALL KTOPDC 3 , ITAP , 0 , IE) 

CALL KTOPDC 4 , JTAP , 0 , IE) 

WRITEC 110,400) 

DO C 1=1, 3300) KBUFC I)=0 
DO (1=1, 133) LAC I) =0 

400 FORMATC 2X, " INPUT 1 FOR OUTPUT OF POLYGON TO TAPE") 

do c 1 = 1 . 60 ) comm i) - 0.0 

RE ADC I IN) I OUT 
WRITEC 1 10) IOUT 
ICUT= 10800 
WRITEC 110,401) 

401 FORMATC 2X, "INPUT 1 TO USE A THRESHOLD") 

READC I IN) 12 

WRITEC 1 10) rz 
IFC IZ. EQ. 1) 


: WRITEC 110,402) 

402 : FORMATC 2X, " INPUT THRESHOLD ") 

: READC I IN) ICUT 

: WRITEC 1 10) ICUT 

: . .FIN 

ACCEPT " START LINE OF CLASSIFIED DATA ? " , I START 

IYMIN= 10000 
IXMIN* 10000 
IYMAX=0 
READC I IN) NV 
WRITEC 110) NV 
201 FORMATC 212) 

DO C 1=1 ,NV) 

: READC I IN, 206) IVXC I) , IVYC I) 

: IVYC I ) = IVYC I ) - ISTART4 1 

.* IFC IVYC I) .LT. 1) IVYC I) = 1 
: IFC IVXC I) . LT. 1) IVXC I) = 1 

: WRITEC I TO) IVXC I) , IVYC I) 

: IFC IVYC I) .LT. IYHIN) IYMIN= IVYC I) 

: IFC IVXC I ) . LT. IXMIN) IXMIN= IVXC I) 

: IFC IVYC I) .GT. IYMAX) IYHAX=IVYC I) 

S..FIN 
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206 


2B 


211 


100 


IVXC NV+1) = IVXC 1) 

IVYC NV+ 1) = IVYC 1) 

IDEL= IYMAX- IYMIN+ 1 
NRITEC 1 10) IDEL 
IFC lOUT.Ea. 1) 

: JBUFC I)" IDEL 

: CALL MTDIOC 4, 50000K+33Q0, JBUF , IS, IE) 

: . . F IN 

ISKIP= C IYMIN- 1 ) *2+ L 
IFt ISKIP.GT.0) 

: FORMAT* 214) 

: DO ( K- 1 , ISKIP) 

: : CALL IlTDIOC 3,0, JBUF, IS, IE) 

: : . .FIN 

i . .FIN 
ISKMYIIIN-1 
DO ( J= 1 , IDEL) 

: DO ( 1= 1 , 3300) JBUFC I) =60 

: DO (1 = 1,3300) KBUF( I) =10000 

: CALL MTDIOC3,0, JBUF, IS, IE) 

: CALL MTDIOC 3, 0, KBUF, IS, IE) 

: ISW=2 

: CALL F0LY2C I YHIN+ J- 1 , ISW, I VX, IVY, NV, LA) 

: NSTR=LAC 1) 

: JV= 1 

: WITEC 10)NSTR 

: DO ( NSEO= 1 , NSTR) 

: : LSUB=2*NSEQ 

: : L2=LSUB+ 1 

: : JS=LA( LSUB) 

: *. IFCJS.LT. 1) JS=1 

: : JF-LAC L2) 

: : NB.ITEC II0,28)JS, JF, JV 

: : F0RMATC2X, "JS, JF, JV ",3110) 

: : JS1=JS-1 

: : IFCJSl.NE.O) 

: : : DO CK2=JV, JS1) JBUFC K2)=6S 

: : : DO CK3=JV, JS1)KEUFCK3) = 10000 

: : : . .FIN 

: : JV= JF+ 1 

: : DO ( J1=JS, JF) 

: : : IFCKBUFC Jl) .GT. ICUT) JBUF( Jl)=60 

: : : LS1= JBUFC Jl) 

: : : IFCLS1.GT.60)LS1=60 

: : : COUNTC LS 1 ) = COUNTC LS 1 ) +1 . 

: : :..F1N 

: : . .FIN 

: JF1= JF+1 

: DO ( Jl-JF 1,3300) JBUFC Jl)= 60 

: DO ( J1=JF1, 3300) KBUFCJ1)= 10000 

: IFC IOUT.EO. 1) 

i : CALL MTDIOC 4, 5000OK+330G, JBUF, IS, IE) 

: : CALL MTDIOC 4, 50000K+3300, KBUF, IS, IE) 

: : VRITEC 10,211) J 

: : FORMAT! 2X, "LINE # M5," PROCESSED") 

: : . .FIN 

: . .FIN 

DO (It- 1,60) 

t ACRE= COUNTC K) * 1 , 05 

: miTEC 110,100) IC, COUNTC K) , ACRE 

: FORMAT! 2X ; "CLASS # ",II0," SIZE ",Flt.0,“ ACRES ",F11.0> 

: . . F IN 

CALL MTDIOC4,6000OK, JBUF, IS, IE) 

CALL MTDIOC 4, 60000K, JBUF, IS, IE) 

STOP 

END 
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C w 

G GSIG (SUBROUTINE) ft 

C ft 

C GETS GLASSIFICATION SIGNATURES ft 

C NUN I IS INPUT FILE (NAMES OF SIGNATURES TO BE RETRIEVED) ft 

C NUN2 IS SIGFILE FROM WIIICH SIGNATURES WILL BE TAKEN ft 

C ft 

Cftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftft 
C ft 

G CREATED AT GEORGIA TECH EES * 

C :l: 

G PROGRAMMER: NICICOLAS L. FAUST ft 

C ft 

Cftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftft 
SUBROUTINE GSIGC A,B, NS)IG, ITEMS IG, N1 ,N2,N3,N4, APR) 

DIMENSION AC 4, 60) ,B(4,4,G0) ,SM(4) ,C0VC4,4) , ITEMS IGC 34) 

DIMENSION NIC GO) , N2( 60) ,N3(60> ,N4(60) ,APRC60) 

IP= 12 
NUN 1=2 
NUN2=4 


DO (1=1 , NSIG) 

: READC NUN 1 , 100) NA1 , NA2 , NA3 , NA4 , APRI 

: N!< I) =NA1 

: N2(I)=NA2 

: N3CI)=NA3 

: N4( I)=NA4 

: APR( I ) = APRI 

: CALL FOPEN C NUN2 , I TEMS IG, U B") 

: READ B INARYC NUN2) I TAP 1 , ITAP2 , 1 TAP 3 , ID ATI , IDAT2 , IDAT3 

» I S ET= 0 

: WHILE! ISET.EQ.O) 

: : READ B INARYC NUN2,END= 120)NS1 ,NS2,NS3,NS4 

: : READ B INARYC NUN2) NP , SM, COV 

: : IFCNS1 .EQ.NA1 . AND.NS2.EQ.NA2) 

: : : IFCNS3.EQ. NA3. AND. NS4. EQ. NA4) 

: : : : DO (K= 1,4) 

: : : : : AC K, I) =SM( IO 

: : : : : DO (L=l,4) 

:::::: BCK,L, I)=COVCK,L) 

: : : : : : . .FIN 

J..FIN 

: : : : ISET= 1 

: : : : CALL FGLOSEC NUN2) 

: : : * , .FIN 

: : : . .FIN 

: :..FIN 

: . .FIN 
RETURN 

100 FORMAT! 4A2, F10. 8) 

120 STOP FILE NOT FOUND 
END 


noannoonfioon 


U^r *7* rfv »|> rj^ »f» •{ * rf> if* .[» rT*-f»<1» r>r{»rf> r|»rf.itv<t.«f. rf» <7* 


HISTO (SUBROUTINE) * 

j|j 

THIS SUBROUTINE CALCULATES AND DISPLAYS A HISTOGRAM FOR A * 

SPECIFIED DATA SET - 4 CHANNELS = 4 HISTOGRAMS * 

# 

* 

CREATED AT GEORGIA TECH EES * 

H: 

PROGRAMMER: NICHOLAS L. FAUST * 

% 

SUBROUTINE HISTCNUII) 

DIMENSION I COUNTC 4 , 100) ,HAXC4) ,MINC4> 

COMMON/DTRANS/INODLC 256) , IXC512,4) , IDUMC 512) 

COMMON /HISTO/ ICOUNT 
ND=4 

DO ( IJ=1, 100) 

: DO C JI = 1 , ND> 

: : ICOUNTC J I , I J) =0 

: :..FIH 

: . .FIN 
DO (1=1, ND) 

: DO ( K= 1 , NUTD 

: : J=IXCK, I) 

: : IFC J.ECLO) J=99 

: : IFC J.GT.99) J=99 

: : ICOUNTC I , J) = ICOUNT( I , J) + 1 

: : . .FIN 

: . .FIN 

DO (KCH=1,ND) 

: MIN( KCH) - 1000 

: MAXC KCH) =0 

: DO ( KEC= 1 . 60) 

: : IFC ICOUNTC KCH, KK) . LT.HINCKCH) > MINCKCH) = ICOUNTC KCH, KK) 

: : IFC ICOUNTC KCH, KK) . GT. MAXC KCH) ) MAXCKCH) = ICOUNTC KCH, KK) 

: : . .FIN 

: IFC MANX KCH) . Ed. 0) MAXC KCH) = 1 

: DO CJJ=I,60) 

: : ICOUNTC KCH, JJ) = C ICOUNTC KCH, JJ) *20) /MAXC KCH) 

: S..FIN 

: ..FIN 
DO (L= 1 , 20) 

: LI=256-L 

: DO CIC= 1,512) IDUMCK)=0 

: DO CKCH=1,4) 

: : ISUB= C KCH- 1 ) *60+NJ 

: : IFC ICOUNTC KCH, NJ) . GT. L) IDUMC ISUB) « 140 

: J..FIN 

: . .FIN 

IDUMC 1 ) = 140 

IDUMC 60) = 140 

IDUMC 120) = 140 

IDUMC 180) = 140 

CALL CPACK1C IDUM, IMODL, O) 

TOBENC L . NE . 1 > CALL IMNRITEC 0 , LI , IMODL , 256 ) 

ELSE 

: DO <L2= 1,255) 

: : IMODLC L2) = 1062 14K 

: : . .FIN 

: CALL IMNRITECO, LI, IMODL, 256) 

: . .FIN 

ACCEPT "HISTO ENDED, ENTER 1 TO CONTINUE : MSTOP 

RETURN 

END 
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C 

G HISTPR 

C 

C THIS PROGRAH PRINTS HISTOGRAMS AND STATISTICS 

G FOR ALL SIGNATURES IN A FILE 

G 

G SEQUENCE : HISTPR HISTFILE SIGFILE 

G 


C * 

G CREATED AT GEORGIA TECH EES * 

C tj: 

C PROGRAMMER: NICHOLAS L. FAUST X 

C x 


cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


DIPIENSION I COUNT! 4 , 100) , IFIL! 20) , IPIAX(4) , TCOL! 123) 
DIMENSION JFIL! 20) ,C0VC4,4) ,SPIC4) 

INTEGER Z, ISWC 2) 

COPfflON /LBLANK/IBL 
DATA IBL /III / 

CALL OPEN! 1 , "COM. CM" , 1 , I ERR) 

GALL COEARGC 1 , IFIL, ISW, I ERR) 

CALL GOMARGC 1 , IFIL, ISW, I ERR) 

CALL GOPIARGC 1, JFIL, ISW, I ERR) 

ACCEPT " WHICH OUTPUT DIVICE DO YOU WANT (10 OR 12) ?",Z 
CALL F0PENC2, IFIL, “B") 

CALL F0PENC3, JFIL, "B") 

READ BINARY! 2) IT1 , IT2, ITS, ID1 , ID2, ID3 
READ BINARY(3) IT1, IT2« IT3, ID1, ID2, ID3 
WRITE! Z, 100) ITl , IT2 , ITS 
WRITE! Z, 101) ID1 , ID2, ID3 
FOREVER 


READ B IN ARY! 2 , ENB= 1 20 > NA I , NA2 , NA3 , NA4 

WRI TEC Z , 1 02) NA 1 , NA2 , NA3 , NA4 

READ BINARY! 2) I COUNT 

READ B r NARY! 3) NA 1 . NA2 , NAS , NA4 

READ B 1 NA11Y! 3 ) NP , SM, COV 

WRITE! Z, 102) NA1 , NA2 , NAS , NA4 

WII TEC Z) NP , SH, COV 

BO C 1=1,4) I-PIAXC I)=0 

DO C 1= 1,4) 

: DO C J= 1, 10O) 

: : IF! ICOUNT! I , J) . GT. IMAXC I) ) IPIAXC I) = ICOUNT! I , J) 

: : . .FIR 

: IF! IPIAKC I ) . EQ. 0) I PLAN! I) - 1 

: DO ( JJ= 1 , 60) 

: : ICOUNT! I,JJ)=C ICOUNT! I, JJ)# 1©) /IPIAX! I) 

: : . . FIN 

: . . F IN 
DO (10=1,2) 

: KG=2*K 

: DO ( Kl = 1 , 10) 

: : K2= 1 1-K1 

: : DO (K3= 1 , 60) 

: : : IQ0-K3+2 

t : *. ICOL! K10) “ IBL 

: : : Kll=63+K3 

: : : IGOLCICl 1) = IBL 

: : : IF! ICOUNT! KO, IC3) . GE. IC2) IGOLCICl 1 ) = IH# 

: : : K12=K0-1 

: : : IF! ICOUNT! IC12.K3) .GE.K2) ICOLCK10) = 1H* 

: : : ..FIN 

: : ICOL! 2) “ 1HI 

: : ICOL! 63) = 1HI 

: : WRITE! Z, 104) ICOL 

: : . . FIN 

: DO ( KQ= 1,123) ICOL! KQ) = 1H= 

: WRITE! Z, 104) ICOL 

: WRITE! Z, 105) 

: WRITE! Z, 105) 

: WRITE! Z, 105) 

: WRITE! Z, 105) 

: WRITE! Z, 105) 

: . .FIN 81 


At it At 4> At At A> 


100 

101 

102 

104 
106 

105 
120 


FORMAT! 2X, "TAPE IDENTIFIER — \3I2) 
FORMAT! 2X, "SIGNATURE TAKEN ON - ",12,"/" 
FORMAT! 1X,4A2) 

FORMAT! 123A1) 

FORMAT! 4A2) 

FORMAT!/////) 

STOP 

END 


12, VM2) 


82 


V 


« *• 


ORIGINAL PAGE IS 
C 

INITGL (SUBROUTINE) 


G 
G 
G 
C 
G 
C 

c 
c 
c 
c 
c 

G 
G 

c 

G 
C 
G 
C 
C 
C 
C 

c 

G 
C 
C 
C 
C 

c 
c 
c 
c 

C * 

C CREATED AT NASA/ JSC < ASTEP) V 

G * 

SUBROUTINE INITCLC VM, NVG, NVM, KNO) 

DIMENSION VM( 240) , NVG( 20) 

COMMON/ INOUT/ I OUT, INF 
1040 WRITE! I OUT , 1050) 

1050 FORMAT ( 38H CHOOSE VALUES FOR INITIALIZATION FROM) 


PURPOSE: 

INITIALIZES THE MEAN VECTORS AND WEIGHTS FOR THE CLUSTER 
ING ALGORITHMS 

DESCRIPTION OF PARAMETERS 

INPUT 

SNAMELIST-S ININIT 

VII - VECTOR MEANS 

NVM “ NUMBER OF VECTORS IN VM 

NVG - WEIGHT ASSOCIATED WITH EACH MEAN 

KNO - DIMENSION OF EAGH VECTOR IN VM 

HOLLERITH OPTIONS 

ZERO - SETS INITIAL MEANS AND WEIGHTS EQUAL TO ZERO 
OLD - RETURNS WITH PREVIOUSLY COMPUTED MEANS AND WEIGHTS 
NEW - ALLOW MEANS AND WEIGHTS TO BE INPUT 

OUTPUT 

GALLING SEQUENCE 

VM - VECTOR MEANS 

NVM - NUMBER OF VECTORS IN VM 

NVG - WEIGHT ASSOCIATED WITH EACH MEAN 

PRINT 

DISPLAYS INPUT IF CARD OPTION IS NEW 
REMARKS AND RESTRICTIONS 

AFTER N AMEL IST INPUT A CARD WITH THE WORD YES IS REQUIRE 
INDICATING CORRECT INPUTS IF NOT REPEAT NAMELIST INPUTS 


* 

* 

* 

* 

* 

:fc 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

X 

At 

* 

•i* 

A: 

X- 

* 

X 

!|C 

A: 

* 1 * 

X 

* 

* 


WRITE! I OUT, 1055) 

1055 FQRMATC 15H ZERO OLD NEW) 

RE ADC INP, 1030) IAN 
103© FQRMATC A6 ) 

IFC IAN.EQ. 4HZER0) GOTO 1100 
IFC IAN. EQ. 3 HOLD) GO TO 1200 
IFC IAN.EQ. 3HNEW) GO TO 1150 
WRITE! TOUT, 1060) IAN 

1060 FORMATC IE ,A6,23H IS NOT A VALID CHOICE.) 

GO TO 1040 
1100 CONTINUE 

DO 51 1=1,240 
51 VM(I)=0.O 

DO 52 1=1,20 
52 NVGC I ) =0 
NVM= 1 

GO TO 1200 
1150 CONTINUE 
1170 WRITE! TOUT, 1160) 

1160 FORMAT! 2 1H S ININIT VM, NVG, NVM) 

WRITE! I OUT, 1022) NVM, ( NVGC I) , 1= 1 , NVM) 

1022 FORMATC 711 NVM = , I2,7H NVG = ,2016) 

CALL MATPRT! VM, KNO , KNO , NVM, 5HMEANS) 

WRITE! TOUT, 1020) 

1020 FORMATC 32H TYPE YES IF INPUTS ARE CORRECT. > 
READ! INF, 1030) LAN 
IF! IAN. NE.3KYES) GOTO 1170 
1200 CONTINUE 
RETURN 


END 


f 
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G :;: 

C INPUTS I G * 

G * 

C ALLOWS MANUAL INPUT OF SIGNATURE FILE * 

C ;Jc 

C « 

G CREATED AT GEORGIA TECH EES * 

G :;: 

C PROGRAMMER: NICKOLAS L. FAUST # 

G tf. 

DIM ENSIOM COVC4.4) ,XPI!4) 

GALL OPEN( 3 , " INSG 1 ’ , 0 , IE) 

GALL F0PEN14, M HSIG'\ "B") 

READC3)RSIG 

ID 1-01 

ID2= 22 

ID3=77 

IT1-99 

IT2=99 

IT3=99 

WRITE BINAPlY!4) IT1, IT2, ITS, ID 1, ID2, ID3 
BO ! L- 1 , NSIG) 

WRITE! 10,101) 

RE AIK 3, 10O)NA1,NA2,NA3,NA4 
WRITE! 10, 102) 

READ(3)NP 
WRITE! 10, 103) 

READ! 3) KM 
DO (10=1,4) 

WRITE! 10, 1©4) K 
READ! 3) (COVCIC.J) ,J-1,4) 

..FIN 

WRITE R INARYC 4) NA1 , NA2 , NA3 , NA4 
WRITE BINARY! 4) NP, KM, GOV 
.FIN 

100 FORMAT! 4A2) 

101 FORMAT! 2X, " INPUT NA1 , NA2 , NA3 , NA4 " ) 

102 FORMAT! 2X, ° INPUT NUMBER OF POINTS") 

103 FORMAT! 2X, " INPUT ! XM( I > , I “ X , 4) " > 

104 FORMAT! 2X, " INPUT ROW # ",11," OF COV - CGOVCK, J) , J= 1 , 4) ") 

STOP 

END 
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G * 

G L INBOX (SUBROUTINE) * 

C :■: 

C USED TO CALCULATE EXTREMES FOR A LINEAR CLASSIFICATION X 

C * 

C * 

C CREATED AT GEORGIA TECH EES * 

C :|: 

C PROGRAMMER: NICKOLAS L. FAUST * 

C x 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
SUBROUTINE LINBOXC AMEAN, BCOV, NSIG, ND) 

DIMENSION AMEAN ( 4 , 10) ,BC0V(4,4, 10) 

INTEGER DELL IN 
COMHOR/BGX/BELL INC 4,10,2) 

DO ( IC= 1 ,NSIG) 

DO ( 1=1 ,ND) 

D IF= C SQRTC BCOV( I , I ,IO )/2.) *3 , 

DELL IN ( I , K, 1)=AKEAN( 1,10 -DIF 
DELLINC I , K, 2) = AMEANC I , IO +DIF 
. .FIN 
FIN 
RETURN 


END 
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G * 

G LINGLASS * 

C * 

G TIHS ROUTINE CLASSIFIES POLYGONS WITH A LINEAR CLASSIFIER * 

G AND CREATES THRESHOLD ARRAY * 

C * 

C SEQUENCE: L INCLASS INPUT DATAP SIGFIL LINTAP * 

C * 

C 

C CREATED AT GEORGIA TECH EES 

C 

0 PROGRAMMER: NICKOLAS FAUST 4' 

C * 

DIMENSION THRESHC812) , ISWC2) 

INTEGER DELL IN 

COMMON /A/’' AMEANC4, 10) ,BC0V<4,4. 10) , IVXC 101) , IVY( 101) ,BC4Q0) ,DET( 10) 

1, LAC 133) ,NVG( 10) ,COUNT( 10) , IBUFC 1650) ,11(30) , 12(30) , 13(30) , 14(30) 

2, IORDERC 10) , NAM1 ( 10 > ,NAM2t 10) ,NAM3( 10) ,NAM4( 10) 

CGMMON/B/JBUF( 4, 810) , IDUPK 8 12) 

CQMMQN/BOX/DELL IN( 4 , 10 , 2) 

IP=12 


C 

C FILE STATEMENTS 

C 

CALL OPEN( 1 , "COM. CM" , 1 , I ERR) 

CALL COHARGC 1 , 1 1 , IS W, IERPc) 

CALL COMARGt 1,11, ISW, I ERR) 

CALL C0I1ARG( 1 , 12, ISW, I ERR) 

CALL COMARGt 1,13, ISW, I ERR) 

GALL COMARGt 1,14, ISW, I ERR) 

C 

GALL 0?EN(2, 11,0, IE) 

GALL MTOPD ( 5 , 14 , O , IE) 

C 

C OPEN TAPE FILE 

G 

GALL MTGPDC3, 12,0, IE) 

ND=4 

NUN1=2 

NUN2=4 

RE ADC 2,200) NS IG 
WRITEC IP)NSIG 
C 

C GET SIGNATURES FOR CLASSIFICATION 

C 

CALL GS IG( AMEAN , BCOV , NS I G, NUN 1 , NUN2 , 13, I ORDER, NAI11 , NAII2 , NAPI3 , NAM4) 
C 

IYMIN=1 00O 

CALL LINBOXC AMEAN, BCOV, NS1G, ND) 

IXMIN- 100© 

IYPIAX=0 
WRITEC IP, 202) 

READ( 2 , 200) NC 
DO ( 1=1, NC) 

: WRITEC IP, 205) I 

: READC2,206) IVXC I) , IVY( I) 

: WRITEC IP, 206) IVXC I) , IVY( I) 

: IFC IVYC I) . LT. IYMIN) IYMIN- IVY( I) 

: IFC IVXC I) .LT. IXMIN) IXMIN= IVXC I) 

: IF( IVYC I ) , GT. IYMAX) IYMAX- IVYC I ) 

: . .FIN 

IVXC NC+ 1 ) - IVXC 1 ) 

IVYCNG+ 1) = IVYC 1) 

NV-NC 

IDEL= IYMAX- IYMIN 
ISKIP= IYMIN+1 
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g OF POOR QUALITY 

G SKIP RECORDS 

G 

CALL MTDIOC 3 , 30000K+ ISKIP , IBUF, IS, I EE) 

DO (LL= 1,812) 

: IDUTK LL) =0 

* T? TN 

IDUIK 1 ) = IDEL 

IDUM! 2) = IXMIN 

FRITE! IP) IDUM! 1) , IDUM! 2) 

CALL HTDIOC 5 , 50O00K-i*812 , IDUM, IS, IER,NF> 

WRITE! IP) NW 
DO ( J= 1 , IDEL) 

: GALL MTD 10(3,0, IBUF, IS, IE, NW) 

; ISW=2 

i CALL“P0LY2( I SKIP* J- 1 , ISW, IVX, IVY, NV, LA) 

: JS=LA(2) 

: JF=LA( 3) 

: CALL UNPAC4! IBUF, JBUF, JS,JF, 810, 1650) 

: JO= JF-JS+ 1 

C : 

C : 

G : 

: NN=NSIG 

: GALL LINEAR! JBUF , NN, IG , ND , NN, AUEAN , BGOV, JO, IDUM, NVG, THRESH) 

G l 

C : 

G : 

! WHITE! 10,207) J 
: JM0VE=JS-1 

: DO CLL- 1 ,810) 

: : LV=813-LL 

: : I DUIK LV) ~ IDUIK LV-2 ) 

: : THRESH! LV)= r rHRESH( LV-2) 

: : FIN 

: IDUHC 1 ) = JIIOVE 

: IDUIK 2) = JQ. 

: CALL MTDIOC 5 , 5O000K+S12, IDUII, IS, IER) 

: GALL MTDIOC 5, 500O0K+8I2, THRESH, IS, IER) 

: TOT=TOT+JQ. 

: DO ( K- 1 , JO) 

: : LSUB= IDUIK K+2) 

: : COUNT! LSUB) = COUNT! LSUB) + 1 

: :..FIN 

t . .FIN 

DO C IO 1 , NS T G) 

; PER= C COUNT! K)*1©0) /TOT 

: WRITE! 12 , 102) K, COUNT! K) , PER, KAMI ! K) , NAM2C K) , NAII3CK) , NAH4C K) 

:..FIN 

GALL MTD I Q( 5, 6000OK, IDUM, IS, IER) 

GALL MTDIOC 5, 60000K, IDUII, IS, IER) 

102 FQRMATC2X, "GLASS ",14," NUMBER ",F7.0," PERCENT " ,F4. 0,2X,4A2) 

103 FORMAT! 2X, 4A2) 

200 FORMAT! 12) 

202 FORMAT! 2X, “INPUT * OF CORNERS 11 ) 

203 FORMAT! 2X, “INPUT CORNER ^ , 13 , " J, I “) 

206 FOPiMAT! 214) 

207 FORMAT! 2X, “LINE * ",I5 , “ PROCESSED " ) 

STOP 

END 
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LINEAR C SUBROUTINE) 

CLASSIFIES LANDS AT DATA ¥ITH LINEAR DISTANCE 

kU 

CREATED AT GEORGIA TECH EES 
PROGRAMMER: NICHOLAS L. FAUST 


SUBROUTINE LINEAR! V, NVM, I COUNT , ND , NS IG, AMEAN , BCOV , NL, IDUM,NVG, 

1 THRESH) 

D I MENS I ON AMEANC ND , NVM) , BCOVC ND , ND , NVM) 

INTEGER DELL IN 
DIMENSION V! ND,NL) 

DIMENSION IDUM! 812) , THRESH! 812) ,NVG! 10) 

INTEGER V 

COMMON/BOX/DELLIN! 4 , 10,2) 

ICOUNT= 1 

DO ! LL= 1 , NL) IDUM! LL) =9 
DO ! LL“ 1 , NL) THRESH! LL) ~ 0 . 

JPT=0 

FOREVER 

JPT-JPT+1 

IF! JPT . GT . NL) RETURN 
IMIN= 1 
DIST= 1 . E+6 
DO ! KM, NVM) 

ESUM=0 
SUM= 0 . 0 
DO ( L= 1 , ND) 

IF!V!L, JPT) . GT. DELLIN! L, K, 1) . AND. V! L, JPT) . LT. DELLIN! L,X, 2) ) 
SUM=SUM+ABS! AMEAN! L, ID ~ V! L , JPT) ) 

KSUM= KSUM-f 1 
..FIN 
. .FIN 

IFIKSUM.EQ.4) 

IF! SUM. LT. DIST) 

D I ST- SUM 

imin=k: 

. .FIN 
. .FIN 
. .FIN 

THRESH! JPT) =D 1ST 
IDUM! JPT) = IMIN 
NVG! IMIN) “NVGC IMIN)+1 
..FIN 

100 FORMAT! 2X, "ENTERED LINEAR M6," TIMES") 

101 FORMAT! 2X, "EXITED LINEAR") 

RETURN 

END 
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C * 

G LISTALL * 

C * 

G PRINTS ALL SIGNATURES , VERTICILS , OR HISTOGRAMS * 

G FOR A SPECIFIC FILE * 

C * 

C SEQUENCE: LISTALL FILE * 

C * 

C * 

C CREATED AT GEORGIA TECII EES * 

c ^ 

C PROGRAMMER: NICICOLAS L. FAUST * 

C * 

DIMENSION IFILC20) ,SK(4) ,C0V!4,4) , IC0UNTC4, IOO) 

DIMENSION IVX! 101) , IVY( 101) 

DIMENSION ISWSC2) , FIELD! 17) 

CALL OPEN( I , "COM. CHM , IE) 

CALL COMARG! 1 , FIELD , ISWS, IE) 

CALL COMARG! 1, FIELD, ISWS, IE) 

CALL F0PENC2, FIELD, "B") 

WRITE! 10,555) 

READ C 1 1 ) I OUT 

£55 FORMAT! 2X, " INPUT OUTPUT DEVICE # ( 10, 12) ") 

WRITEC 10.401) 

READ! 11) ICON 
WRITEC 10,403) 

READ C 1 1 ) ITOT 

READ BINARY! 2) ITAP1, ITAP2, ITAP3, IDAT1 , IDAT2, IDAT3 
WRITEC I OUT, 100) ITAPI, ITAP2, ITAP3, IDAT1, IDAT2, IDAT3 
FOREVER 

READ B INARYC 2 , END= 1 20 ) NS 1 , NS2 , NS3 , NS4 
IF C ICON . EQ. 1 ) READ B INARYC 2 ) NP , SII, GOV 
IFC ICON. EQ. 2) 

READ BINARY! 2) IC, IYHIN, IYTL4X 
DO (10=1,10 READ BINARY! 2) K2, IVXCK1) , IVYCIC1) 

..FIN 

I F C I CON . EQ. 3 ) READ B I NARY! 2 ) I COUNT 
WRITE! IOUT, 10 1 ) NS 1 , NS 2 , NS3 , NS4 
IFC ICON. EQ. 1) 

J IF! ITOT.EQ. 1) WRITEC IOUT) NP, SII, COV 
2 p . F IN" 

IF! ICON. EQ. 2) 

IFC ITOT.EQ. 1) 

WRITEC IOUT) K. IYM.IN, IYMAX 
DO C Kl= 1 , K> WRITE! IOUT) IC1 , IVX! Id ) , IVYC IC1) 

-FIN 
.FIN 

IFC ICON. EQ. 3) 

: IFC ITOT.EQ. 1) WRITEC IOUT) ICOUNT 

: . .FIN 
..FIN 

FORMAT! IX, 6 12) 

FORMAT! 1X,4A2/) 

FORMAT! 2X, " INPUT SWITCH , 1 -SIG , 2 “VER , 3 -HIS") 

FORMAT! 4 1 10) 


100 

101 

401 

402 

403 
120 


FORMAT! 2X, " 
STOP 
END 


TOTAL LIST - 1 , OR NAME LIST - 0 ?'') 
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G * 

G LISTS IG * 

C * 

C PRINTS ALL SIGNATURES FOR A GIVEN FILE * 

C * 

C SEQUENCE: LISTS IG FILE * 

C * 

G * 

G DESIGNED AT GEORGIA TECH EES * 

G # 

G PROGRAMMER: NIGKOLAS L. FAUST * 

G * 

DIMENSION IFIL( 20) ,SH(4) ,C0V(4,4) , ISWC2) 

GALL OPENC 1 , "COIL CM" , 1 , I ERR) 

GALL CGEARG! 1 , IFIL, ISU, I ERR) 

CALL GOIIARGC 1 , IFIL, ISW, I ERR) 

CALL F OPEN 1 2 , IF IL , n B 11 ) 

READ BINARY! 2) I TAP I , ITAP2, JTAP3, I BAT I , IDAT2, IDAT3 
WRITE! 12, 100) ITAP1 , ITAP2, ITAF3, IDAT1 , IDAT2, IDAT3 
FOREVER 


READ B INARY( 2 , ENB= 1 20 ) NS 1 , NS2 , NS 3 , NS 4 
READ B INARYC 2) NP , SM, GOV 
WcITEC 12, 101)NSl,NS2,HSS f NS4 
■WRITE! 12) NP , SH, GOV 


: . .FIN 

100 FORHATC IS, 6 12) 

101 FORMAT! 1K.4A2) 
120 STOP 

END 
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*1» »J-» *!> *!» -fct> 4fji.r> Oj «hi Ui *L» vi» tU til <i« «U.I. «t_» vL. ^L> Oi fcli >i< (i« til kt« >!• <4> ip «tj lUll* klj «b »U >]» tlj ttj >1* »t. « !■» »U .L.»U«I> 


c ^ 

G MATIN' (SUBROUTINE) * 

G 

G MATRIX INVERSION PROGRAM * 

G * 




G 

G 

G 

C'J” 


c 

c 


G 

c 

c 


CREATED AT NAS A/ JSC ( ASTEP) 


* 

* 

at 


atataiatatatatatatataiaiatatatatatatatatattitfftataiaiat^atatatattitaiatatHtut^atatatat^iftatataiatatatiitatatataiatHtatatatatatatiitaiataiatiitat 

SUBROUTINE HATINC A, N,B,K,ICEY, DETERM) 

DIMENSION A( N,N) ,B(N> , IPIVOTC 20) , INDEX(20,2) 

DOUBLE PRECISION PIV0TC20), T, SNAP , DETERM, AMAX, ZERO, A, B 


INITIALIZATION 
3 KEY = N 
10 DETERM= 1.D+0 
15 DO 20 J= 1 ,N 
20 IPIVOTC J) =0 
30 DO 550 1=1, N 

SEARCH FOR PIVOT ELEMENT 


40 AMAX= Q.D4-0 
45 DO 105 J= 1,N 
50 TF ( IPIVOTC J)-l) 60,105,60 
60 DO 100 IC= 1 , N 

70 IF (IPIVOT(IO-l) 80,100,740 
80 IF CDABS(AMAX)-DABS(A( J,IO)) 85,100,100 
85 I RON= J 
90 ICOLUM=K 
95 AIIAX=A(J,IO 
100 CONTINUE 
105 CONTINUE 

ZERO = 1.0-16 

1 07 IF( DABSC AMAX) - ZERO? 745 ,745,1 10 
110 IP I VOTC ICOLUH) = IP I VOT( ICOLUH) + 1 


C 

G INTERCHANGE ROWS TO PUT PIVOT ELEMENT ON DIAGONAL 

C 

130 IF ( IRON- ICOLUH) 140,260,140 
140 DETERM=-BETERM 
150 DO 200 L=1,N 
160 SWAP=AC IRON, L> 

170 AC IRON, L) =A( ICOLUM, D 
200 AC ICOLUM, L) =SNAP 
205 IFCH) 260,260,210 
210 SNAP = EC IRON) 

230 BC IRQ Tv) = BC ICOLUM) 

250 BC ICOLUM) = SNAP 

260 INDEXC 1,0 = IRON 

270 I NDEXC I , 2 ) = I COLUM 

310 PI VOTC I)=AC IGOLUIL ICOLUM) 


91 


c 

C DIVIDE PIVOT HOW BY PIVOT ELEMENT 

G 

330 AC ICOLUM, ICOLUM) = . ID-f 1 
340 DO 350 L= 1 , N 

350 AC ICOLUM, L) =AC ICOLUM, L) /PIVOTC I) 

355 IF(M) 380,380,360 

360 BC ICOLUM) = BC ICOLUM) /PIVOTC I) 

C 

G REDUCE NON-PIVOT ROWS 

C 

380 BO 550 Ll= 1, N 

39Q I F C L 1 — I COLUH) 400,550,400 

400 T= AC LI, ICOLUM) 

420 ACL1 , IGOLUII) = .0D+0 
430 DO 450 L“ 1 , N 

450 AC LI , L) =A( LI , L) -AC ICOLUM, L) *T 
455 IF CM) 550,550,460 
460 BCLI) = BC LI) - BCICOLUM)*T 
550 CONTINUE 
C 

C INTERCHANGE COLUMNS 
C 

600 DO 710 1*1, N 
610 L=N+1-I 

620 IF C INDEXCL, 1>-INDEXCL,2) ) 630,710,630 
630 JROW* INDEXCL, 1) 

640 JCOLUM* INDEXCL, 2) 

650 BO 705 K= 1 , N 
660 SWAP* AC IC, JROW) 

670 AC K, JROW) = AC K, JCOLUM) 

700 AC K, JCOLUM) *SWAP 
705 CONTINUE 
710 CONTINUE 

DO 800 1*1, N 
J=N+1-I 

80O DETERM* BETERM*P IVOTC J) 

74© RETURN 

745 DETERM = 0.B+0 

746 KEY =1-1 
750 RETURN 

END 
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‘.HMM 


ORIGINAL RAGS IS 
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P '}' ‘4^ *> *)* *1» •> ili O* *1# »U *. 1 * »b »u»Jj <L* « v ^ ^ 4 1 *i* »X- \V 4* •i' 4/ 4» 1 I 4 Uji *i» U* tl* U< *1^ <u *1* *J* »!• il> *1* iJ» «}<■ *4/ *,i* sl> »>• *A- * 1 * »A* *1* 

G Jft 

G MATPRT * 

G ;jc 

C PRINTS A REAL MATRIX * 

C X 

G jjj 

C GREATER AT RASA/ JSC CASTEP) # 

G X 


SUBROUTINE MATPRTC MAT, NR, RD , NC , NAM) 
INTEGER RD,P,G 
REAL MATC KB, NO 
100© FORMATC IH , X 2 , 6F 10 . 3) 

1010 FORMATC 1H ,6119) 

KOUT= 12 

WRITEt KOUT, 1020) NAM, NR, NG 
1020 FORMATC 1H , 15X, A2, 16, 3H BY, 13) 

DO 100 1=1, NR 
DO 100 J= 1 ,NG 
IFC MATC I , J)>150 , 100 , 150 
100 CONTINUE 

WIITEC KOUT, 1030) 

1030 FORMATC 1H , 15X, 12S ALL ZEROES. ) 

GO TO 30 
150 CONTINUE 
P=0 
G=~5 

10 P=F+6 

Q=Q.+6 

IF(NG.LT.P)P=NG 


miTECicouT, ioioxj, j=a,P) 
mi TEC KOUT, 1000) 

DO 20 1=1, NR 

WI TEC KOUT, 1000) I, C MATC I, J) ,J=a,P> 
20 CONTINUE 

IFCNG.GT.P)GO TO 10 
30 CONTINUE 

WRITECKOUT, 1000) 

RETURN 

END 
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G * 
C MESS «i* 
C ^ 
C UNPACKS AND DISFLAYS MULT I SPECTRAL 12 CHANNEL * 
C DAEDALEUS SCANNER DATA AND DISPOSES REFORMATED DATA * 
C TO TAPE * 
C * 
C SEQUENCE: MESS MTU:F MTU:F * 

c * 




C 

c 

c 

c 

c 


CREATED AT GEORGIA TECH EES 
PROGRAMMER: MICHAEL D. FURMAN 


* 

# 

# 

* 

5k 


C5S^5k^5k5k5lci!:5k5k^5fC5}:5k5k^^:^5k^#5iC^^55}:5|:5!:5k5ic^:5K5k^:^C5k*#^:5}:^}i::kH55it5}:5t;^^:i:<t:^5k^5l:5t:^5i:5i;S5^5l:5|;:!:5!i:}:^;^:5lJ5i;57:5}: 


INTEGER INPUT! 700) , E0UTC850) , ITAPE1! 10) , ITAPE2! 10) , IMAG2C 1700) 
INTEGER IMAG3! 1100) , ISW!2) 

CALL OPEN! 1 , "COM. CM" , 1 , I ERR) 

CALL COMARG! 1 , ITAPE1 , ISW, I ERR) 

CALL COMARGC 1, ITAPE1, ISW, I ERR) 

GALL COtlARGC 1, ITAPE2, ISW, I ERR) 

CALL MTOPD! 2 , ITAPE 1,0, IE) 

CALL MTOPD( 3 , ITAPE2 , O , IE) 

ACCEPT "START WITH WHAT ELEMENT ", IEL 
FOREVER 

: CALL MTDI0C2, 0, INPUT, IS, IE, IC) 


L= 1 


DO ! Jl- 1 , 565,3) 

: M0= IANDC ISHFTC INPUT! Jl) , 2) , 17740010 

: Ml- IOR! M0, I AND! ISHFT! INPUT! Jl) , 6) , 30010 ) 

: I OUT! L) - IOR! Ml , ISHFT! INPUT! J 1+ 1) , - 10) ) 

: M2- IOR! ISHFT! INPUT! J 1+1) , 10) , I AND! ISHFT! INPUT! J 1+2) ,-6) , 140010 

: I OUT! L+ 1 ) = IOR! M2 , I AND! ISHFT! INPUT! J 1+2) ,-2) ,377K) ) 

: L=L+2 


: . .FIN 

CALL UPAC8! IOUT, IMAG2,850) 

DO 1 1- 1 ,40) 

: IF ! ! IMAG2! I) .EQ.255) .AND.! IMAG2! 1+1) .LE. 10)) IPASl=I+4 

: . .FIN 


IPAS2-0 

DO ! I=IPAS1,512) 

: IF ! ! IMAG2! I) . EQ. 255) .AND. ! IMAG2! 1+1) ,LE. 10) ) IPAS2= 1+4 

: . . FIN 
L30- 1 

WHEN ! IPAS2.GT.0) 

: DO ! I-IPAS2,756) 

: : IHAG3CL30) = IMAG2! I) 

: : L30=L30+1 

: J..FIN 

: . . FIN 
ELSE 
: L31=l 

: DO ! 1= IPAS1, 756) 

: : IMAG3! L3 1 ) = IMAG2! I) 

: : L31=L31+1 

: : . .FIN 

:..FIN 

CALL PACO! IMAG3! IEL) , IOUT, 512) 

CALL RIKWRITE! 0 , 0 , IOUT, 256) 

CALL MTDIO! 3, 50000K+850, IMAG3, IS, IE, IC) 

..FIN 


STOP 


END 
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C 

c 

G 

G 

C 

G 

C 

C 

c 

G 

G 

G 

C 

C 

c 

G 

C 

G 

G 

G 

G 


MODIFY 

PURPOSE 


(SUBROUTINE) 


* 

+U 
*»♦ 

* 

* 

* 
•j- 

<J> 
* 
m 
* 
* 
at 
at 

•j* 
4l> 

«v< 
* 
* 
sk 
* 
* 
sfs 

G ^ 

C GREATED AT NASA/ JSC CASTEP) * 

G jj: 

SUBROUTINE MODIFY ( VI, V2, Nl, N2, ND ) 

D I MENS ION VI ( ND) , V2 ( ND) 

REAL Nl 
XN1 = CN1) 

XN2 = FLOAT(K2) 

SN1R2 ~ 1 . 0/( N I+FLOATC N2) ) 

BO 10 I=I,ND 

10 VI { I) = ( XN1*V1(I) + XN2#V2( I) )*XN1N2 
Nl * Nl + FL0ATCR2) 

RETURN 
END 


COMPUTES WEIGHTED AVERAGE OF TWO MEAN VECTORS 
DESCRIPTION OF PARAMETERS 
INPUT 

CALLING SEQUENCE 

VI - 1 ST VECTOR MEAN 
V2 - 2 ND VECTOR MEAN 

Nl — NUMBER OF VECTORS USED TO COMPUTE VI 
N2 - NUMBER OF VECTORS USED TO COMPUTE V2 
ND - DIMENSION OF VI AND V2 

OUTPUT 

GALLING SEQUENCE 

VI - LIGHTED AVERAGE OF INPUT VI AND V2 

Nl - NUMBER OF VECTORS USED TO COMPUTE OUTPUT VI 
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dC - 

‘•a 






G 

C 

G 

G 

G 

C 

c 

c 

c 

G 

G 

C 

C 

C 

C 

c 

G 

G 

G 

C 


PACK (SUBROUTINE) 

PURPOSE 

PAGES A STORAGE ARRAY TO ELIMINATE A VACATED SLOT, MOVES 
ALL VECTORS WITH INDEX GREATER THAN INDEX OF VACATED 
SLOT DOWN ONE POSITION IN THE ARRAY 

DESCRIPTION OF PARAMETERS 

INPUT 

CALLING SEQUENCE 

V - DATA ARRAY 

ND - DIMENSION OF EACH VECTOR IN V 

NV - NUMBER OF VECTORS IN V 

IND - INDEX IN V OF VACATED SLOT 

OUTPUT 

CALLING SEQUENCE 

V - PACKED ARRAY 


* 

* 

* 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

*1* 

X 

X 

X 


cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

c * 

C CREATED AT NASA/ JSC (ASTEP) * 

c * 




SUBROUTINE PACKC V,ND,NV, IND) 
DIMENSION VCND.NV) 

IF( IND.EQ.NV) RETURN 
11= IND+1 
DO 20 1= 1 1 ,NV 
12“ I— 1 

DO I© J=1,ND 
1© VC J, 12) =V< J, I) 

20 CONTINUE 
RETURN 
END 
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G * 
G P0LY2 (SUBROUTINE) * 
C * 
G CALCULATES POLYGON INTERSECTIONS FOR LINE AND POLYGON * 
C * 

C * 
G CREATED AT NASA/ERL * 
G * 
C PROGRAMMER: RONNIE PIERSON * 
C # 


SUBROUTINE P0LY2C JY, ISW, IX, IY,N,LA) 

C VERTEX -C IXC I) , IYC I) ) = (ELEMENT, SCAN) FOR ITH VERTEX. 

DIMENSION IXC 101) , IYC 101) ,DC 100) ,FC 101) ,LAC 133) ,SC 101) 

IFC ISW. EO. 1) GO TO 1 
C SLOPE STORE BLOCK X/Y 

DO 20 1=1, N 

IFC IYC I) . EQl. IYC 1+ 1 ) ) GO TO 20 

DC I) = FLOAT! IXC 1+1) -IXC I) ) /FLOAT! IYC I+D-IYC I) ) 

20 CONTINUE 

1 Y= FLOAT! JY) 

M=0 

NOV=0 

DO 40 1=1, N 

C LOCATION OF INTERVALS SUCH THAT IYC I) .LE. JY.LT. IYC 1+ 1) . SCAN=JY 

IFC IYC I) , EO. JY) GO TO 30 
IFC IYC I) .LT. JY)GO TO 22 
IFC IYC 1+1) .LT. JY)GO TO 28 
GO TO 40 

22 IFC IYCI+1) .GT, JY)GO TO 28 

GO TO 40 

C POLYGON BOUNDARY COUNTER AND FUNCTIONAL VALUES. -NON VERTICES- 

23 H-M+ 1 

FC M) =D< I)*C Y-FLOATC IYC I) ) ) +FLOAT( IXC I) ) 

GO TO 40 

C POLYGON BOUNDARY COUNTER AND FUNCTIONAL VALUES. -VERTIGIES- 

C NO V= NUMBER OF VERTICES 

30 M=M+1 

NOV= NOV+ 1 

FC M) = FLOAT! IXC I) ) 

LAC NOV) = IXC I) 

40 CONTINUE 

K=M 

C BOUNDARY SORT - INCREASING ORDER. 

DO 50 J=1,K 
SC J)=FC 1) 

IND= 1 

DO 48 L= 1 , M 

IFCFCL) .GE.SCJ))GO TO 48 
SC J) =FCL) 

IND=L 

48 CONTINUE 

FC IND) = 10000 
50 CONTINUE 

G FOR CN0V.LT.2) INTERVAL SC J) ,SC J+l) MEMBERSHIP GAN BE DETERMINED 

IFCNOV. EO. ©)G0 TO 500 
IFC NOV. EO. 1) GO TO 490 

G DATA SHIFT FOR CONSISTENCE AFTER MEMBERSHIP IS DETERMINED 

DO 52 K=1,M 
52 FCIO=S(K) 

NO I =11- 1 
M=G 
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c 

G INTERVAL FCK),FCK+1> MEMBERSHIP SECTION FOR CN0V.GE.2). 

C 

DO 200 K=l,NOI 

G A IS THE POINT CHECKED TO DETERMINE INTERVAL F(K) ,F(K+1) MEMBERSHI 

A=(F(K) +FCK+1) )/2. 

54 L=0 

J=0 
LV=0 
JV=0 
IN=0 

DO 122 1=1, N 

C LOCATION OF INTERVALS SUCH THAT IXt I) .LE.A.LT. IXC 1+1) 

CX= FLOAT! IX( I) ) 

DX= FLOAT! IXC I+D) 

IFC CX. Ed. DX) GO TO 122 
IFC GX. Ed. A) GO TO 12B 
IFCCX.LT. A) GO TO 124 
IFC DX. LT. A) GO TO 128 
GO TO 122 

124 IFC A. GE.DX) GO TO 122 

128 CY= FLOAT! IYC I)) 

DY= FLOAT! IYC I+D) 

DM= ( DY-CY) / C DX-CX) 

C BOUNDARY POINTS < A t FC A) ) FOR GIVEN INTERVAL IX! *1) .LE.A.LT. IX! 1+11 

FX=DM&! A-GX) +CY 
C IN= 1 INDICATES F!A)=JY 

IF(FX.Ed.Y) IN= 1 
IFC FX. LT. Y) GO TO 130 
C J INDICATES FCAl.GT.JY. 

C 

J=J+1 

IF C FX. RE . C Y) GO TO 122 
IF! A. NE. GX) GO TO 122 
C 

C JV INDICATES Ft A) = IYC I) .GT. JY, FOR SOME I. 

G 

JV=JV+1 
GO TO 122 
C 

C L INDICATES FCAl.LT.JY. 

C 

130 L=L+1 

IFC FX. NE. CY) GO TO 122 
IFCA.NE.CXIGO TO 122 
C 

G LV INDICATES F( A) = IYC I) . LT. JY. FOR SOME I. 

C 

LV=LV+ 1 

122 CONTINUE 
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G 

C INTERVAL FCK) ,FCK+1) MEMBERSHIP DETERMINATION. 

G 

IFC IN.EO. 1>G0 TO 150 
IFC JV.ECL0)GQ TO 140 
J=L 

IF(LV.EQ.0)GG TO 140 
G 

G IF NO DETERMINATION GAN BE MADE, ADD ANOTHER INCREMENT AND TRY AGP 

G 

A=A+.0t 
GO TO 54 
140 L-2*CJ/2) 

G 

C IFC J.NE.L) THE INTERVAL FCK) ,FCK+1) IS IN POLYGON. 

G 

IFC J. NE.L) GO TO 150 
G 

G CHECK FOR VERTEX AT CFCK) ; JY) 

G 

DO 146 J=l,NOV 
CX^FLOATCLAl J) ) 

IFCFCK) .NE.GX)GO TO 146 
G 

G VERTEX C FC K) , J Y) I NCLUDED . 

G 

M=M+1 

SCM)=F(K) 

M-M+l 

SCM)=FCK) 

146 CONTINUE 

IFCK. NE.NO I) GO TO 200 
FX=FCK+1) 

DO 14B J= 1 , NOV 
GX=FLOATC LAC J) ) 

IFCFX.NE.CX)GO TO 148 

M-M+l 

SCH) =FX 

M-M+l 

SC H) =FX 

148 CONTINUE 

GO TO 200 
C 

G INTERVAL FCK) ,FCK+1) INCLUDED 

G 

ISO M=M+ 1 

SC M) -FCK) 

M=M+1 

SC M) -FCIC+1) 

200 CONTINUE 

GO TO 500 
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V 


c 

G ALL INTERVALS FORTHCOMING ARE BOUNDARY TO BOUNDARY ON JY 

C 

C CHECK FOR INTERVAL THAT CONTAINS ONLY A VERTEX 

C 

490 L=2*(M/2) 

IFCL. EQ. H) GO TO 500 
VAV= FLOAT! LA( 1) ) 

M=H+1 

K-M 

492 K=K-1 

IF ( S ( K) . EQ. VAV) GOTO 494 
SCK+l)=SCK) 

GO TO 492 
C 

C STOP INSERT FOR VERTEX LOOP 

G 

494 SC KM )= VAV 

G 

G IFCM.EQ.0)NO POINTS ON SCAN JY ARE INCUDED 

C 

500 IFCM.EQ. 0) GO TO 555 

1=0 

C 

G LOOP START AND STOP ROUND OFF TO MIN AND MAX INTEGER VALUES INCLUB 

C 

501 1=1+1 

LA< I)=S( I) + .999999 
1 = 1+1 

LACD=S!I) 

IF(I.LT.M) GO TO 501 
1 =-’ 1 
G 

G I NTERVAL C START. GT. STOP ) COMPRESS 

C 

504 1= 1+2 

505 IF(LACI) .LE.LACI+1DGO TO 510 

M=II-2 

IFtJL EQ.. 0)GO TO 555 
IF! I. GT. ID GO TO 511 
DO 506 K=I,M 

506 LAC IQ =LAC K+2) 

GO TO 505 

510 IFC C 1+1) ,LT. M) GO TO 504 

511 1=0 
G 

G CLOQPSTOP ON INTERVAL! J> . EQ. LOOP START ON INTERVAL! J+ 1 ) ) COMPRESS 

C 

512 1=1+2 

510 IF C I . EQ.. ID GO TO 555 

514 IFC LA! I) .NE.LAC 1+1) )GO TO 512 

M=M-2 

DO 516 K=I,Ml 
516 LAC K) =LAC K+2) 

GO TO 513 
G 

C LAC 1) = NUMBER OF LOOPS RETURNED 

G 

555 IFCM.EQ.0) GO TO 560 

DO 557 1=1, M 
J=M-I+1 

557 LA! J+1)=LAC J) 

560 LA!l)=M/2 

RETURN 
END 
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G # 
G POLYGL * 
G # 
C PERFORMS CLUSTERING OPERATION * 
c ^ 

G SEQUENCE: POLYCL INPUT OUTPUT * 
C * 


P*fr^^* 0 # ^*^ ^^^*^^^^^*^^^^^*f t *J**i**P^^ l f* , J**i* t P^* t J 4 *i , *l <l *f**^ l J**I < ^**t ,l> t* < ^**^ t J* l J**^*t > *J*^**^'^l**I* l, ^*t**« > *^*t**t 1|l i**^*i* , i»*^«J»*^*^^ |t ^*^l | J i *i**y**J* , J* 


C 

G CREATED AT NAS k / JSC (ASTEP) 

C 


* 

* 

# 




COHKON/URTNUM/ IMGUNT, DATUNT 
COMPION/BLANIC/ICNUM 
INTEGER CAR,BUF, ISWC2) 

REAL NVG 
COMMON/COL/ JS , JF 

COMUON/BOUND/YHIN , XMIN, YMAJI, XMAX 
INTEGER YIIIN , YHAX, XMI N , XHAX 
DIMENSION JPXC 101) , JPYC 101) , LAC 133) 
COmON/NVEC/NVM 
COMMON/BUFFER/BUFC 3300) 

COHMON/A1/ I0UTCB12) 

INTEGER DATUNT 
COMMON/ INOUT/NOUT, NIN 
COMKON/CHARAC/CAR 
COMMON/D IST/ ID 1ST 
INTEGER PL IST 

DIMENSION VMC240) , ADC20,20) 

DIMENSION NVGC20) ,FLIST(20) ,CARC20> , IACOPC 4) 
DIMENSION RMC20) t RV(20) 

DIMENSION VHPC60) , VARC60) 

DIMENSION ACRE! 20) , 11(30) , 12(30) 


C=3, 


CALL OPENC 1 , "COM. CM" , 1 , 1 ERR) 
CALL COMARGC 1,11, ISW, IERR) 
CALL GOMARG( 1,11, ISW, IERR) 
CALL COMARGC 1,12, ISW, IERR) 
CALL HTOPD( 4 , I 1 , © , IER> 

CALL MT0PD(3, 12,0, IER) 

RP=3 


S=l. 

Rl=20. 

R2=20. 

NVHHAX= 15 
NPT= 100 
NET=500 
NMT= 100 
IP=0 

IACOPC 1) =2HME 

IACOPC 2) =2HSI 

IACOPC 3) =2IIAN 

IACOPC 4) =2HQU 

ND=4 

NRT= 1 

IMGUNT=3 

DATUNT=4 

NOSCAL=l 

NBUFSZ= 326 0 

IQ-NBUFSZ/CND+2) 

I BUF2= NBUFSZ- IQ* 1 
IBUF1= IBUF2-IQ 
NIN- 11 
NOUT= 10 


101 


CAR! 1) = 1HA 
CARC2) = 1H, 

CAR!3) = 1H: 

CAR( 4) = 1H- 
CAR!5) = lH, 

CAR(6) = 1H/ 

CAR! 7) “ 1H+ 

GAIK 8) = 1110 
GAR( 9 ) = 1HX 
GAIK 10) = 1H# 

GAIK 11) = IBS 
CAR( 12) = 1H% 

CAR! 13) = 1HB 
CAR! 14) = 2H@ 

CAR( 15)=1H* 

C USER INPUTS 

G 

1 CONTINUE 

NRITECNOUT, 1003) 

1008 FORMAT! 2X, * INPUT VALUES FOR G, RP, Rl, R2 AND NVMMAX’ ) 
READ! N IN) C , KP , R1 , R2 , NVMMAX 
NRITE! NOUT) C , RF , Rl , R2 , NVMMAX 
miTE( NOUT, 102) 

102 FORMAT! 22H TYPE YES IF INPUTS OIO 
READCNIN, 104) IAN 
104 FORMAT! A6) 

IF! IAN.NE.3HYES) GOTO 1 
IF ( NVMMAX. GT . 20 ) NVMMAX* 20 
C 

C FIRST PASS INITIALIZATION 

C 

GALL INITCL! VM, NVG, NVH, ND) 

NRITE! 10, 199) 

199 FORMAT! 2X, "INPUT 1 FOR POLY , 0 - NORMAL") 

READ! 11) I POLY 

IF! IPOLY.NE. 1)G0 TO 201 

WRITE! 10,202) 

202 FORMAT! 2X, "INPUT # OF CORNERS") 

READ! 11) NV 
DO 204 1=1, NV 
WRITE! 10,205) I 

205 FORMAT! 2X, "INPUT CORNER 13, " J, I ") 

204 READ! 1 1 ) JPX! I ) , JP YC I ) 

JPX(NV+1)=J?XC 1) 

JFY( NV+ 1) = JPYC 1 ) 

201 CONTINUE 

JMIN= 10000 
IMAX=0 
IMIN= 10000 
DO 405 K1=1,NV 

IF! JPX! Kl) . LT. JMIN) JHIN= JPX! K1 ) 

IF! JPY! Kl ) . LT. I MIN) IMIN= JPYC IC1 ) 

405 IF! JPY! IC1 ) . GT. I MAX) IMAX= JPYC Kl ) 

NRT= IM4X- IHIN+ 1 

DO 406 K2= 1,812 

406 IOUT! K2) =0 

I OUT! 1) =NRT 
IOUT! 2) = JMIN 
LINIT= IMIN 
LIN= IMIN+1 

GALL MTD I O C 3 , 50000K+8 12 , IOUT, ISl, IER) 

CALL MTDI0(4, 30000K+LIN, IOUT, ISl, IER) 

INRT=0 

NEC*© 

NHC=0 
NPG=0 
IPASS= 1 
JPTP=0 
RMC1)=0. 

RV! 1) =0. 

DO 8 1=1, ND 
8 VMPt I)=0. 

DO 10 1=1, NVM 
10 PL 1ST! I) = I 


102 


G 

G 

C 


FIRST PASS PROCESSING 


G 

G 

C 


20 


INRT= INRT+1 

IFC INRT.GT.KRT) GO TO 50 
IL=LIN3T+INRT~ 1 
IS¥= IL+ I 

GALL P0LY2C IL, ISW, JPX, JPY,NV,LA) 
JS=LA( 2) 

JF=LA(3) 

IF( JS . Eft. 0) JS" 1 

WRITEC NOUT) JS , JF , INRT, " 1 » 

GALL UNPAG1 
NX=JF-JS+1 
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GALL CLUSTAC BUF C I > , VM, ND , NX, NVM, NVMMAX, NVG, C , 5 , RP , R1 , NPC , 
*NPT,PLIST,NEC,NET,NHIN,NMC,NMT, IPASS,BUF< IBUFl) ,BUFC IEUF2) 
s!rIP , JPTP , RM, RV , VMP , VAR) 

JPTP= JPTP+NX 
GO TO 20 


» 


SPECIAL ELIMINATION AND MERGER TESTS 


50 KX=0 
JPTP=0 
NEG=RET 
NMG=NMT 

CALL CLUSTAC BUF CD, VM, ND , NX, NVTI, NVMMAX, NVG, G , S , RP , Rl , NPG , NPT 
*PLIST, NEC, NET, RUIN, NMC,NKI\ IPASS,BUF( IBUFl) , 

*BUFC IBUF2) , IP , JPTP , RM, RV, VMP , VAR) 

C 

C SECOND PASS INITIALIZATION 

G 

INRT=0 

NBACIC= L I N+NRT 

GALL LTD 10(4, 40000IONBACK, BUF , IS, IE) 

FRITEC 10,222)NBACK 

222 FORMATC 2X, “ BACKSPACED ”, 110, " LINES") 

NEC=0 

IFC LIN.LE.0) GO TO 444 

CALL MTDIO(4,30000K+LIN,BUF, IS, IE) 

444 CONTINUE 
NHC=0 
NPG=0 
IPASS=2 
DO 60 1=1, NVH 
PLISTC I)=I 
66 NVGC I) =0 


G 

G SECOND PASS PROCESSING 

C 

62 INRT= INRT+1 

IFC INRT. GT.NRT) GO TO 70 

IL=LINIT+INRT-1 

ISW= IL+ 1 

CALL F0LY2C IL, ISW, JPX, JPY, NV, LA) 

JS=LAC2) 

JF=LA( 3) 

JQ= JF-JS+1 
DO 402 LL- 1,812 
402 IOUTCLL) =0 

IOUTC 1) = JS 
I0UT(2)=Ja 

TvRITEC NOUT) JS , JF, INRT, "2 " 
rFCJS.EO-.O) JS=1 
GALL UNPAG1 
NX= JF-JS+1 

GALL CLUSTACBUFC 1) , VM, ND , NX, NVM, NVMMAX, NVG, C,S,RP,R2, NPG, NPT 
*PLIST, NEC , NET, NXIN, NMG , NMT, IPASS, BUFC IBUFl) , 

*BUFC IBUF2) , IP , JPTP , RM, RV, VMP , VAR) 


JPTP= JPTP+NX 

DO 403 LL= 1,810 

L2= IBUF2+LL— 1 

IFC L2 . GT . NBUFSZ) GO TO 404 

403 IOUTC LL+2) =BUFC L2) 

404 CONTINUE 

GALL MTDIOC 3 ,5000GK+812, IOUT, ISt, IER) 


G PRINT RESULT SUMMARY 

C 

70 WRITE! ROUT, 106) 

106 FORMAT( IX, ’CLUSTER SYMBOL SIZE R MEAN R SIGMA 
KS AVE= 0 . 0 
DO 72 1=1, RVM 
RV! I)=SQRT!RV< I)) 

pp- i 

IFC NOSCAL.EQ.. 0) GC= 1.53046 
ACRE( I) =NVG! I)*CC 
IFCRVGC I) .LT.NSAVE)GO TO 2000 
NSAVE=RVG! I) 

KNUM= I 

2060 CORTIRUE 

72 WRITE! ROUT, 108) I, CAR! I) ,RVGC I) ,RM( I) ,RV! I) , ACRE! I) 
108 FORMAT! 16 , 7S, Al , F 10 . 0 , 2X, 2F7 . 2 , F 10 . 2) 

I=RD*NVM 
DO 74 J=1,I 
74 VM!J)=VHP(J) 

CALL MTDIOC 3, 6O000K, I OUT, IS1, IER) 

GALL MTDIO! 3, 6O000K, IOUT, IS1 , IER) 

G 

G USER OPTION SELECTION 

C 

76 WRITE! NOUT, 110) 

110 FORMAT! 19H CHOOSE OPTION FROM) 

WRITE! ROUT, 1 12) IAGOP 
112 FORMAT! 1H,4A8) 

READ! NIN, 104) IAN 
IF! IAN, EGL5HMEANS) GO TO 80 
IF! IAH.EGL6HSIGM4S) GO TO 85 
IF! IAN.EGL6HANGDIS) GO TO 90 
IF! IAR. EGL4HQDTT) STOP 
WRITE! NOUT, 114) IAR 

114 FORMAT! 1H ,A6,22H IS NOT A VALID CHOICE) 

GO TO 76 
C 

C MEANS DISPLAY 

C 

80 CALL MATPRT! YM, ND , ND , NVM, 5HMEANS ) 

GO TO 76 
G 

C SIGMAS DISPLAY 

C 

85 I=RD*NVM 
DO 86 J=1,I 

86 VAR! J)= SORT! VARCJ)) 

CALL MATPRT! VAR, ND, ND , NVM, 6HS IGMAS) 

GO TO 76 
C 

C ANGDIS DISPLAY 

C 

90 IDISF=2 

GALL ANGDISCVM, NVM, ND, IDISF, AD) 

CALL MATPRT! AD , NVM, NVM, NVM, 6 HARGD IS) 

GO TO 76 
END 


ACRES ’ ) 


X04 
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c 

C PUTIN 

C 

C INPUT SIGNATURES 

G 

C 

C CREATED AT GEORGIA TECH EES 

C 

G PROGRAMMER: NIGKOLAS L, FAUST 

C 

DIMENSION C0VC4.4) ,SMC4) 

CALL 0PENC3, “IRSC",©, IE) 

CALL FOPENf 4, "HSIG", "B") 

READC 3) NSIG 

ID 1 = 0 1 

ID2=22 

ID3=77 

IT1=99 

IT2=99 

IT3=99 

WRITE B INARYC 4) ITI, IT2, IT3, ID1, ID2, ID3 
DO t L= 1 , NSIG) 

: WRITEC 10, 101) 

: READC 3 , 100) NA1 , NA2 , NA3 , NA4 

: WR.ITEC 10, 102) 

: READC 3) NP 

: WRITEC 10, 103) 

: READC 3) XM 

: DO C IC= 1 , 4) 

: : WRITEC 10 , 104) K 

: : READC 3) C GOVC K, J) , J=l, 4) 

: : , .FIN 

: WRITE B INARYC 4) NA1 , NA2 , NA3 , NA4 

: WRITE B INARYC 4) NP, XM.GOV 

: . .FIN 

100 FORMATC 4A2) 

I © 1 FORMATC 2S , n INPUT NA 1 , NA2 , NA3 , NA4 « ) 

102 FORMATC 2X, "INPUT NUMBER OF POINTS”) 

1 03 FORMATC 2K, " INPUT ( SMC I ) , 1= 1 , 4) " ) 

104 FORMATC 2S, “INPUT ROW # Ml," OF COV - C COVCK, J> , J= 1 , 4) ") 

STOP 

END 


105 


*********** 


V 


f 


C * 
C RAINBOW * 
C »{i 
C DISPOSES EXISTING PSEUDOCOLOR MEMORY TO A DISK FILE * 
G DISPOSES PSEUDOCOLOR MEMORY IN DISK FILE TO DISPLAY X 
C x 


C x 

C CREATED AT GEORGIA TECH EES X 

G 5f£ 

C PROGRAMMER: FRED L- THOMPSON X 

C * 

INTEGER IORANC 94) 

CALL FOPEN ( 5 , "RYGBV. ") 

ACCEPT "TYPE 1 TO SAVE PRESENT PSEUDOCOLOR ", IANS! 

IF ( I ANSI -Ed. I) 

REWIND 5 

DO (J= 1,1000) READ BINARY (5,END=10) IORAN 
10 : DO (1=1,94) IORANC I) =0 

CALL RCM (0, IORAN) 

TYPE "TYPE A 30 CHARACTER DISCREPTION" 

TYPE " . " 


20 


30 


4© 


: READ (11,20) ( IORANC I) , 1=65, 94) 

: FORMAT (30A1) 

: WRITE BINARY (5) IORAN 

: TYPE "PSEUDOCOLOR MEMORY NUMBER ",J 

: . .FIN 

ACCEPT "TYPE 1 TO REPLACE PRESENT PSEUDOCOLOR " , IANS2 
IF (IANS2.EGL 1) 

: ACCEPT “TYPE 1 IF YOU KNOW TEE PSEUDOCOLOR MEMORY NUMBER 

: IF ( IANS3.NE. 1) 

: : PAUSE PUT UP PSEUDOCOLOR TEST PATTERN HIT RETURN 

: : REWIND 5 

: : DO ( J=0 , 1000) 

: : : READ BINARY C5,END=30) IORAN 

: : : WRITE (10,40) < IORANC I) , 1=65, 94) 

: : : FORMAT (1X,30A1) 

: : : CALL WCM ( 0 , IORAN) 

: : : ACCEPT "TYPE 1 FOR NEXT PATTERN " , IANS4 

: : : IF (IANS4.NE.1) STOP NORMAL EXIT 

: : : . .FIN 


IANS3 


60 


: . .FIN 

IF CIANS3.EGL 1) 

*. REWIND 5 

: ACCEPT "PSEUDOCOLOR MEMORY NUMBER ? 

: DO C J= 1, INUM) READ BINARY (5) IORAN 

: WRITE (10,60) ( IORANC I) , 1=65 ,94) 

: FORMAT CIX,30A1> 

: CALL WCM (0, IORAN) 

: . .FIN 


: . .FIN 

STOP NORMAL EXIT 
END 


" , INUM 
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G 

C 

G 

C 

G 

C 

G 

C 

C 

G 


RATIO 

THIS PROGRAM DIVIDES ONE ERTS CHANNEL INTO THE OTHER 
THREE TO LEAVE THREE CHANNELS OF DATA. IT ALSO WILL TAKE 
TWO CHANNELS AND BY DIVIDING ONE INTO THE OTHER TWO, CREATE 
A THIRD CHANNEL, ALL TO BE DISPLAYED. 

SEQUENCE: RATIO MTU: F 


* 

* 

* 

# 

* 

* 

Hi 

# 

H'- 

Hi 


C Hi 

G CREATED AT GEORGIA TECH EES * 

C ^ 

G PROGRAMMER: MICHAEL D. FURMAN * 

C * 

INTEGER INPUTC 1700) , IMP(3400) , IMAG( 0:255) , IMAG2C 0: 255) , ITAFE( 10) 
INTEGER ISW<2> 

CALL OPENC 1 , "COM. CM" , I , IERR) 

CALL COMARGI 1 , ITAPE, ISW, IERR) 

GALL COHARGI 1, ITAPE, ISW, IERR) 

CALL HTOPDC 2 , ITAPE , 0 , IE) 

ACCEPT "TWO C 2) OR FOUR (4) CHANNEL RATIO 'MANSI 

WHEN CIANS1.EQ.2) ACCEPT "INPUT NUM. AND DENOM. CHANNELS “ , INUM, I DEN 

ELSE ACCEPT " INPUT DENOMINATOR CHANNEL " , IDEN 

ACCEPT "SCALE FACTOR 'MSCAL 

ACCEPT "INPUT LINE AND ELEMENT ” ,LINE, N 

IF ( ( ( N/2) *2) . EQ. N) N=N - 1 

ISICP= 1+LINE 

IEL=(N*4>-<6/( l+N-( ( N/2) *2) > > -1 

IEND= IEL+ 10 19 

IDEN2=IBEN*2-1 

DO (1=1, 1SK?) CALL MTDIOC2,0, INPUT, IS, IE, IC) 

DO (K2= 1.255) 

CALL MTD 10(2,0, INPUT, IS, IE, IC) 

CALL UPAC8< INPUT, IMP, 1700) 

M2=0 
M1 = 0 

DO <IC= 1,7,2) 

L=0 

IF ( K.NE. IDEN2) 

DO ( J1=IEL, TEND, 8) 

INAG( L) = ( EHPC J 1+10 /< IMPt IDEN2+J 1) + 1 . 0) ) FISCAL 
IMAG(L+1) = ( IMP( Jl+K+l)/( IMP< IDEN2+J1+1) + 1 .0) ) FISCAL 
L=L+2 
.FIN 

DO (1=0,255) 

IF ( IMAG( I) . GT.255) 

TYPE "ATTEMPTED SCALE TO “ , IMAGC I) 

IMAGC I) =255 

.FIN 

. . FIN 

CALL PAC8( I MAG, IMAG2,256) 

CALL RIMWRITE( Ml , M2 , IMAG2 , 256) 

WHEN ( M2. EQ. 0) M2=255 
ELSE M2=0 
Mi = 111+1 
..FIN 
..FIN 
..FIN 
STOP 
END 
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G <ji 

G RDATA :ft 
G * 
G THIS SUBROUTINE UNPACKS DATA FOR A POLYGON * 
C * 


G 

C CREATED AT GEORGIA TECH EES 

G 

G PROGRAMMER: NIGKOLAS L. FAUST 

G 


SUBROUTINE RDATA! ITOP, IBGT,NV, IX, IY, IND,NP> 
DISSENSION LA( 133) , IXC 101) , IYC 101) , JBUF( 100,4) 
COMMON /BUFFER/ IBUF! 1700) 

COMMON/ I NDEV/ I POS 


* 
£ 
* 
a: 

:i : 


IXD(512,4) 


IC=0 

ISUM1- ITOP- IPOS- 1 

222 FORMAT! 2X, "SKIPPED *’,15,'' RECORDS’*) 

IF( ISUM1 . GT. 0) CALL MTDIOC2,30000K*ISUM1, IBUF, IS, IER) 
IF! ISUM1 * LT. 0) CALL MTDIO(2,40000K~ISUM1, IBUF, IS, IER) 
IDEL=IB0T-IT0P+1 
DO ( IL= 1 , IDEL) 

: CALL MTDIGC 2,0, IBUF , IS , IE) 

: ISW= ITOP+1 

: ILIN= ITOP+ IL- 1 

: CALL P0LY2C ILIN, ISN, IX, IY,NV,LA> 

: NR= 100 

: NSEG= 1 

: MM= 1650 

: DO ! K= 1 , NSEG) 

: : ISUB1=2*K 

: : JS=LAC ISUB1) 

: : JF=LA( ISUB1+1) 

: : JD= JF-JS+ 1 

: : CALL UNPAC3C IBUF, JBUF, JS, JF,NN,MM> 

: : DO (K=1,JD) 

: : : IC-IC+1 


: : : DO (KD-1,4) 

: : : : IXDC IG , KD) = JBUFC K, KD) 

: : : s . .FIN 

: : : . . F IN 

: S..FIN 

: ..FIN 
NF=IC 


IPOS= IBOT 

NRITE! 10,222) IDEL 

TYPE " TOTAL # PTS = u t NP 

RETURN 

END 
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*4* ^ *!■* vi« *1* <V «A» «J* «Jj» »!> *Sm Ui <L. lO* tli kUiAoliflj 

C 

C RDIMG (SUBROUTINE) 

G 

G READ AND DISPLAY IMAGE TO SCREEN 

C 

/^ l l**i»*i* , ^ t ^*f*4**]**l**J*^«*f**l*^*f*^*l»*f>*J*4 , ^«^Ml><*«i*»]^*J*»1«<J«^«*f*<l**f*4«»)«»]*»)*4 t *i*^*f*4**f**f*«l»^*^ > ^^J**l* I l»*l«^»*>^^**)**i**i*«J»^*I*^l>*l**l**l i «J*«l**l* , J*^ < 

c 

G CREATED AT GEORGIA TECH EES 

C 

C PROGRAMMER: ROBERT A. MADDOX 

C 




SUBROUTINE RDIMGt IL, IEL,NCH, IBLUP) 

INTEGER IMAG( 260) , IPIAG2< 0:5 12) ,NCH(0:2) 

GOMMON/BUFFER/INPUTC 1700) ,/INDEV/IP0S 

ISUPI1 = IL— IP0S- 1 

IELP=( ( IEL+254) /%) *4 

IEL2= ( ( IEL-S)/2>*4 

IF ( IBLUP, LE.0) 

: TYPE " ILLEGAL FACTOR" , IBLUP 

: RETURN 


: . .FIN 

IF CIStmi.GT-0) CALL MTDIO(2,30000K+ISUMi, INPUT, IS, IE, IC> 
IF ( ISUPI1.LT. 0) CALL PITDIO(2,40000K-I5UM1, INPUT, IS, IE, IC> 
NCH( 0) = I 


100 


ICC=0 
NCH( 1) =2 
NCH( 2) =4 
PI7=255/ IBLUP 
DO (ia=0,M7) 

: CALL MTDI0C2,©, INPUT, IS, IE, IC) 

: DO C 15=0,2) 

: : Mo-0 

: : IF (I5.EQ,. 1) M5=255 

: : I2=NCHC 15) 

: : 13= IEL2+ 12 

: : 14= IELP+ 12 

: : L= 1 


: : DO (J= IQ, 14,4) 

: : : JJ=J 

: : : IF(JJ. GT. 1700) JJ= 1700 

: : : IMAG(L) = INPUT( JJ) 

: : : L=L+1 

: : : . . FIN 

: : IF C IBLUP. GT. 1) BLOW-UP- IMAGE 

: : DO (1=1, IBLUP) 

: : : ECC=ICC+1 

: : : IF( ICC.GT. 768) GO TO 100 

: : : CALL RIPIWRC 15 , M5 , IPIAG, 256) 

: : :..FIN 

: :..FIN 

: . .FIN 
CONTINUE 
IP0S= IL+P37 




RETURN 


TO BLOW-UP- 1 PLAGE 
: CALL UPAC8C IMAG, IMAG2,256) 

: DO ( K7=0 , 255) 

: : K9=256-K7 

: : X£DB=K9/ IBLUP 

: : IPIAG2(K9) = IPIAG2C KDB) 

• 5 * . F IK 

: GALL PACBC IMAG2 , IPIAG, 512) 

: . .FIN 


END 
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RET I MG (SUBROUTINE) * 

* 


PURPOSE * 

RETRIEVES IMAGE AND THRESHOLD ARRAYS CORRESPONDING TO * 

A DATA RECORD * 

* 

DESCRIPTION OF PARAMETERS * 

$f* 

INPUT . * 

CALLING SEQUENCE * 

I FLAG - FIRST ENTRY OR REWIND AND SKIP HEADING RECORD * 

FLAG * 

NUN IT - UNIT NUMBER * 

K; 

OUTPUT Jfr 

CALLING SEQUENCE » 

IMG - IMAGE ARRAY * 

THR - CORRESPONDING THRESHOLD ARRAY * 

NP - HUMBER OF POINTS IN EACH ARRAY * 

* 




* 


CREATED AT NAS A/ JSC ( ASTEP) 


* 

* 




SUBROUTINE RETIMGC I FLAG, NUN IT, IMG, THR, NP) 
DIMENSION IMG< 112) 

COKMON/COL/ JS , JF 
INTEGER THR 

IFC IFLAG.NE. 1) GO TO IO 
REWIND NUNIT 

10 READ B INARY( NUNIT) JS , NP , C IMGC I ) , 1= 1 , NP) 
RETURN 
END 




.r 


OBIOW^ ? NfJ| 
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% 

SAVIKG (SUBROUTINE) * 

* 

PURPOSE * 

SAVES, ON TAPE, IMAGE AND THRESHOLD ARRAYS CORRESPONDING * 
TO A DATA RECORD * 

& 

DESCRIPTION OF PARAMETERS * 

* 

INPUT * 

CALLING SEQUENCE * 

I FLAG - FIRST ENTRY OR REWIND AND WRITE HEADING RECORD * 

FLAG * 

NUN IT ~ UNIT NUMBER * 

IMG - IMAGE ARRAY TO BE SAVED * 

THR - THRESHOLD ARRAY CORRESPONDING TO IMG * 

NP - NUMBER OF POINTS IN EACH ARRAY % 

5 a 

CREATED AT NAS A/ JSC (ASTEP) # 

SUBROUTINE SAVIKGC I FLAG, NUN IT, IMG, THR, NP) 

DIMENSION IMGCNP) ,THRCNP) 

INTEGER THR 
DIMENSION BLKC166) 

COMMON/COL/ JS , JF 
BLK( n-i. 

IF( IFLAG.NE. 1) GO TO 10 
REWIND NUN IT 

10 WRITE B INARYC NUN IT) JS , NP , < IMG< ID ,K= 1, NP) 

WRITEC 12)NP, IMGC 1) 

RETURN 

END 


111 


i- 


C * 

C SCORECARD - 

C * 

C THIS ROUTINE CLASSIFIES TRAINING FIELDS WITH A * 

G MAXLIK CLASSIFIER * 

G * 

C SEQUENCE: SCORECARD INPUT TAP RAM SIGFIL VERTEX * 

C * 

G * 

C CREATED AT GEORGIA TECH EES * 

C * 

C PROGRAMMER: NICHOLAS L. FAUST # 

C ^ * 

DIMENSION AMEANC4,20) , BC0VC 4, 4, 2©) , IVXC 101) , IVYC 101) ,B(400) ,DET(20> 
1 ,LAC 133) ,NVG(20> , ICOUNTC20) , IBUFC 1650) , 11(30) , 12(30) , 13(30) , 14(36) 
2, IORBERC20) , NAM1C20) ,NAM2(20) ,NAH3(20) ,NAM4(20) , ISWC2) 
COMMON/DTEANS/ IDATC 256) , JBUF( 512,4), IDUMC 5 12) 

IP= 12 


C 

C 


FIELD AND OPEN STATEMENTS 


C 

CALL 0PEN( 1 , "COM. CM 11 , 1 , I ERR) 
CALL COHARGC 1,11, ISW, IERR) 
CALL COMARGC 1 , 1 1 , ISW, IERR) 
CALL COMARGC 1, 12, ISW, IERR) 
GALL COMARGC 1,13, ISW, IERR) 
CALL COMARGC 1,14, ISW, IERR) 

C 

CALL 0PENC2, 11,0, IE) 

CALL F0FENC7, "MAXFLD", "B”) 

C 

C OPEN TAPE FILE 

C 

CALL MTOPDC 3,12,0, IE) 

ND=4 


NUN1-2 

NUN2=4 

RE ADC 2,200) NS IG 
WRITEC IP) NSIG 
C 

C GET SIGNATURES FOR CLASSIFICATION 

C 

CALL GSIGC AMEAN,BGOV, NSIG, NUN1 ,NUN2, 13, IORDER,NAMI , NAM2 , NAM3 , NAMI ) 
GALL ICBC AMEAN , BCOV, ND , NSIG, DET, B , B) 

C 


C 

ISIG—0 
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G 


DO C KB IG= 1 , NSIG) 

: ISIG= ISIG+ 1 

i GALL FOPEN! 5 14 ,: B" ) 

: READ BINARY! 5 ) ITAP 1 , I TAP 2, ITAP3 , IDAT1 , IDAT2 , IDAT3 

: K3= I ORDER! KSIG) 

: DO (K2=1,K3) 

: : READ BINARY! 5 ) NA1 , NA2 , 1? A3 , NA4 

: : READ BINARY! 5) K, IYMIN, IYMAX 

: : DO ( 1=1, K) 

: : i READ BINARY! 5) Kl, IVX! I) , IVY! I) 

: : :..FIN 


: . .FIN 

GALL FCL0SE(5) 

WRITE! IP) ITAP1, ITAF2, ITAP3 
WRITE! IP) IDAT1, IDAT2, IDAT3 
WRITE! IP, 103) NA1 , NA2, NA3, NA4 
WRITE! IP)K, IYMIN, IYHAX 
DO C L3= 1 , IO 

: WRITE! IP) Kl, IVX! L3) , IVY! L3) 


ORIGINAL PAGE IS 
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.FIR 


NV=K 


ISKIP= IYMIN+1 


IVX! K+ 1 ) = rvxc 1) 
IVY! K+l) = IVY! 1) 
IDEL= IYMAX- IYMIN+1 
DO (11=1,20) 

: I COUNT! II) =0 

: . .FIN 


JTOT=0 


G 

C 


G 

G 

C 

G 

G 

G 


102 

103 

200 


/*SKIP RECORDS 

GALL KTDIO! 3.30000K+ISKIP, IBUF, IS, IEE) 
WRITE BINARY! 7) IDEL, ISIG 
DO ! J= 1 , IDEL) 

: GALL MTD 10! 3,0, IBUF , IS , IE, NW) 

: ISW=2 

: CALL P0LY2! ISKIP+J-1, ISW, IVX, IVY, NY, LA) 

: JS=LA( 2) 

: JF=LA<3) 

: GALL TJNPAC3! IBUF, JBUF, JS, JF,512, 1650) 

: JG= JF-JS+1 

: WRITE! IP) J, JS, JF, JO, IDEL, ISKIP 

: DO !L=l,JO) 

: : DO !L1= 1 ,ND) 

: : : ISUB=(L-1)*ND+L1 

: : : IBUF! ISUE) = JBUF! L, Ll> 

s : J..FIN 

: : * .FIN 


CALL MAXLIC IBUF , NB , JO., AHEAN , BGOV, NS IG , DET, IE DM, NVG) 


: : JM0VE=JS-1 

: : WRITE B I NARY! 7 ) JMOVE , JQ, IDUM 

: : JTGT= JTQT+ JQ. 

: : DO !K=l,Ja) 

: : : LSUB= IDUM(K) 

: : : ICOUNT(LSUB) = IC0UNT(LSUB>+1 

: : : . . FIN 

: : . . FIN 

: DO (K=1,NSIG) 

: : JPER= ! TGOUNT! K) * 100) ✓ JTOT 

: : WRITE! 12 , 102) K, IGOUNT! K) , JPER, NAM1 C K) , NAM2C K) , NAII3! K) , NAM4l K) 

: : . .FIN 

: ISUM= IDEL+ISKIP 

: GALL MTD 10(3, 400G0K+ISUM, IBUF, IS, IER) 

: . .FIN 

FORMAT! 2X, "CLASS ",I4 , " NUMBER " , 16 , " PERCENT M3, 2X, 4A2) 
FORMAT! 2X, 4A2) 

FORMAT! 12) 

STOP 
END 
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G X 

C SEQCOV < SUBROUT IKE) * 

C :;: 

C PURPOSE * 

C if: 

G RECURSIVE CALCULATION OF SAMPLE COVARIANCE MATRICIES AND * 

G MEAN VECTORS x 

G x 

G INPUT t. 

C * 

G CALLING SEQUENCE * 

C X = CURRENT SAMPLE VECTOR X 

G N = DIMENSION OF X 

G V = PREVIOUS COVARIANCE MATRIX X 

C M = PREVIOUS MEAN VECTOR * 

C I - PREVIOUS # OF SAMPLES X 

C D = SCRATCH AREA * 

C X 

G OUTPUT x 

G :1c 

G V = CURRENT COVARIANCE MATRIX X 

C M - CURRENT ME AIT VECTOR X 

C I = CURRENT # OF SAMPLES X 

C X 

GXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
G X 

G CREATED AT NAS A/ JSC CASTEP) * 

C X 

cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SUBROUTINE SEGCOVC X,N, V, M, I , D) 

INTEGER X 
DIMENSION XCN) 

DIMENSION VCN,N) 

REAL MC N) 

DIMENSION D( NT, 2) 

"WHENC I . LE.0) 

DO CL=1,N) 

V(IC,L) =0.0 
1=1 

RETURN 
..FIN 
..FIN 
ELSE 

FI=FLOAT( I) 

FIP=1./FI 
1 = 1+1 

FI=FLOATCI) 

FI=1./FI 
DO ( K= 1 , N) 

DCK, 1)=X(K)-M(K) 

DCK,2)=D(K, 1)*FI 
M(IO=M(K)+D(K,2) 

..FIN 

FIP= 1 . -FIP 
DO ( K= 1 , N) 

DO ( L= 1 , N) 

: VC K, L) =F IP*V( K, L) +DC K, 2) *D< L, 1 ) 

: . . F IN 
. .FIN 
. . FIN 
RETURN 
END 
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G 

c 

c 

G 

c 

G 

G 

C 

C 

c 

c 

c 

c 

c 

G 

G 

C 

c 

G 

C 


SEQST 

PURPOSE 


(SUBROUTINE) 


SEQUENTIALLY CALCULATED BEAN VECTOR AND VARIANCES 
DESCRIPTION OF PARAMETERS 
INPUT 

GALLING SEQUENCE 

KM - CURRENT MEAN 
XV - CURRENT VARIANCES 

N - NUMBER OF POINTS USED TO COMPUTE XM AND XV 
ND - DIMENSION OF DATA VECTOR 
X - NEW DATA VECTOR TO BE ADDED TO XM AND XV 

OUTPUT 

CALLING SEQUENCE 

XM, XV, N - UPDATED VALUES 


* 

* 

* 

& 

* 

* 

* 

* 

* 

* 

* 

* 

* 

i|t 

& 

# 

* 

* 

* 


*1* *J* *1* +)+ *£+ *1* *t* *J» *1* **■» *1* *L» *1* +3+ *1* *1* fcl* hi# ^S» »3a *lr *3* fcl* fcj* fcj* •i#*V*t**J* *1* *4^ »!» *1* *1* *1* O* «*> %!*'**« 

C :Ji 

C CREATED AT NASA/ JSC CASTEP) * 

C * 




SUBROUTINE SEQSTC XM, XV, N, ND , X) 
D I MENS ION XM( ND) , XVC ND) , X( ND) 
INTEGER X 

IF(N.GT.O) GO TO 29 
DO 10 1= 1 , ND 
XM( I)=X( I) 

1© XV(I)=0. 

N= X 

RETURN 


20 F1=C 1.-1./FL0ATCN) ) 
N=N+1 


F2= 1 ./FLOATCN) 

DO 30 1=1, ND 
T=X( I)-XM( I) 

XM( I)=XMC I)+F2*T 
30 XVC I>=F1*XVC I)+F2*CT**2> 
RETURN 
END 
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C 

C STRIPLESS 

G 

G PROGRAM TO DESTRIPE A DATA SET GIVEN THE 

C CONVERSION TABLE FROM STRNF2 

C 

C SEQUENCES STRIPLESS MTIN: F MTOUTsF TABLE 

C 

»)j kl< 0>*4**I-k vt* ■!««)•* *1**1** *4* ilitl/tli*Jit|j\}j»l*tJ*d>»lj vi* *1**1* »)>*.■' 1 ^>*^*1* *!'*•]»■■ J> "*» * * 

r t‘ % "v* n- 'i* *r* *1* *!*'»>*?’* *t> *r* *r* *T* ■r* rl* <r* *!'• * *#* * , r»«T» W *cj w, » »fwj v v rj-» *f* *i- 

C 

C PROGRAMMERS: RONNIE PIERSON C NASA/ERL) 

C NICHOLAS FAUST (GEORGIA TECH EES) 

G 


* 

a.}* 

P 'i r 

J# 

* 

* 

■<X 

X 

• V 

* 
X 


4 < vl-» «i* »lj. «J> »E* *V t ‘l > «■)* *1 


DIMENSION IEMC4000) ,LT(4,6, 128) . IZ(2,4,412) 
DIMENSION ITAP( 15) , JTAPC 15) , IFIL( 15) , ISW(2) 
CALL OPEN( 1 , "COM. CM" , 1 , IERIl) 

CALL COMARG( 1 , I TAP, ISW, I ERR) 

CALL COHARGC 1, I TAP , ISW, IERIl) 

CALL COMARG( 1, JTAP, ISW, IERIl) 

CALL COMARGt 1 , IF TL, I SW, IERIl) 

CALL F0PENC3, IFIL, "B") 

CALL MT0PDC2, ITAP,0, I ERR) 

CALL MTQPD( 4 , JTAP , 0 , I ERR) 

NOFTl= 1 

READ BINARYC3)LT 
DO (2=1,2) 

: CALL MTDI 0(2,0, IEM, IS, I).RR,NW) 

: CALL MTD 10(4,5 ©000K+NW , JEM, IS, I ERR) 

: . . FIN 

N=0 

NL=© 


TYPE "HOW MANY LINES OF 6 TO SKIP?" 

REA3X 11) KICK 

TYPE “HOW MANY LINES TO PROCESS?" 

RE ADC 1DNLINE 

KICK6=6*KICK 

IF(KIGK.GT.0) 

: CALL MTD 10(2, 30060K+KICE6 , IEM, IS , IERR) 

: . .FIN 

FOREVER 

: CALL MTD I0( 2,0, IEM, IS , IERR, MW) 

: N=N+1 


: NL= NL+ 1 

: IFCNL .GT. NLIME) WRITE-EOF 

: IFCNW .LE. 50) WRITE-EOF 

: CALL PACKER! 1 , IEM, IZ) 

: DO (1=1,4) 

: : DO (J= 1,407) 

: : : K= IZC 1, I, J)+l 

: : : Kl= IZC 2, I , J)+l 

: : : IF(K .GT. 12S)K=12S 

: : : IFCE1 .GT, 128)Kt=128 

: : : IZ( 2, I , J) =LT( I , N, IC1) 

: : : IZ( 1 , I , J) =LTC I , N, 10 

: : : . „ F IN 

: : . .FIN 

: CALL PACKER! 2, IEM, I Z) 

: CALL MTDIOC 4, 50000K+ 165(? , IEM, IS, IERR) 

: IF(N .EQ. 6)N=0 

: . .FIN 

TO WRITE-EOF 

: CALL MTDIOC 4, 600G0K, IEM- IS, IERR) 

: CALL MIR IOC 4 , 6O0OOK, IEM. IS , IERR) 

: CALL MIR 10(4, 1OOO0K, IEM. IS, IERR) 

: CALL MTDI0C2, 1006OK, IEM. IS, IERR) 

: STOP 

s,.FIH- 
END 
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ft ’X‘ "A X 





•}> ii» lU xi> tLrvljilj •> kU »!, vVj tb «J> <ia >V * 1 ^* K±+ ti< u< «l> »J» <vL» iV <1/ »Jrf ki> >b iU »X* fU ^X. >X> K> >I> vU «1/ «X» *^* O' ‘■V tV «lr -vA» * l< «tj «^» vU \U tU *I> »J» tW «Xr iA» >U«ii «i* 
1 J*T*n*rfi«^ n> *■?•, fT> >7> n»!T» «-T-» n» »-f» »T> tT> .TV w*w »f« «yv«f»<yvrt* IT^ fT- ^«T»it*TS^H»f»yHff>iit>^iT>^fi»f>«f*<f»«f< »T» r r* *f»*V , *P 


C * 

G STRIPNF X 

G ?|c 

C TO ELIMINATE STRIPING IN LANDSAT DATA X 

G * 

C SEQUENCE ; STRNF2 TAPE TABFIL * 

G X 

■'Jn *n^^'T>»(»*t*W»»ivrt\'r , 'r*<r» *!■» *T"*j'**T w r* i *'' '**{'* rtv^>*l’«»T> ■‘t > »T*»T» *1'* 

C sfc 

G CREATED AT NASA/ERL * 

G jjt 

G PROGRAMMER: RON PIERSON * 

C X 


DIMENSION IERC620, 8) ,LTC4,6, 128) , IUNPAKC2,4,3> ,DT(4,8, 129) 
DIMENSION IEMC 1650) ,LIMSC 4, 6, 2) 

D I MENS ION ITAP1C 15) , ITAB( 15) 

COMHON/DUM/ IER, DT, LIES, IUNPAK, ITAP 1 
EQUIVALENCE C IEMC 1) , IER( 1,6)) 

EQUIVALENCE CLTC U 1 , 1) , IER< 1,1)) 

GALL FIELBC2, ITAP1 , S701) 

GALL FIELDC3, ITAB,S70l> 

CALL KTOPDC 1 , ITAP 1,0, I ERR) 

CALL F0PENC3, I TAB, "B") ORIGINAL PAGE 18 

N0FTR=0 OF poor QUALITY 

N0FT1= 1 
NOFT= 1 

TvRITEC 10 , 2222) 

2222 F0RMATC2X, "INPUT NPRNT NLINES") 

READ( 1 1) NPRNT.NL INES 
99 CALL MTDIOC 1,0, IEM, IS, XERR) 

CALL MTDIOC 1,0, IEM, IS, I ERR) 

C SIX SCAN READ 

10 N= 1 

I ERROR" 0 

1 1 CALL MTDIOC 1,0, IEM, IS , IERR, NW) 

IFC NW. LT. 1600) GO TO 400 
IFCNLl.GT. NLINES) GO TO 400 

DO CKV- 1 , 620) IERCKV,N) “ IEM(KV) 

N-N+ 1 
NL1=NLX+1 
IF( N. LT. 7) GO TO 11 
C SIX BY SIX PROCESS 

N- 1 
NUII=0 

DO 299 NE- 1,51 
N=N+ 12 
DO 280 J= 1 , 6 
DO C Ll= 1,3) 

: DO < L2= 1 , 4) 

: ; NL=N+(L1-1)&4+L2-1 

: : IWRD=IERCN1, J) 

: : IUNPAICC 1 , L2, LI) = ISHFT( IWHD,— 8) 

: : IUNFAKC2,L2,L1) = IANDC IWRD,377K> 

: : . .FIN 

:..FIN 

C ELEMENT AVERAGE 

DO 280 1-1,4 
DO 200 ICl-1,3 

NUM=NIRU- IUNPAKC 1,1, El ) + IUNPAICC 2 , 1 , El ) 

200 CONTINUE 

DTC I , J , 1) =FLOAT( NRM) /6 . 

NUM=0 

280 CONTINUE 

DO 284 1=1 4 

C DISCARD OF ABNORMAL DATA VALUES 

SH=DTC 1, 1, 1) 

BG=SM 


117 


DO 279 J=2,6 
RAD=DTC I, J, 1) 

IFfRAD.GT.SM)GO TO 275 

SM-RAD 

GO TO 279 

275 IF ( RAD . LT. BG) GO TO 279 

BG=RAD 

279 CONTINUE 
RAD=BG-SM 

IF( RAD.GT. 10. )GO TO 284 
C DATA AGREGATION 

AVG= 0 . 

DO 282 J=I,6 

282 AVG= DTI I , J , 1 ) + AVG 
MN= AVG/6 . + . 5 

IFt MN. LT. 0) GO TO 284 
IFC HN. GT. 127) GO TO 284 

mn=mn+2 

D= DTI 1,7, UN) +6 . 

DO 283 J= 1 , 6 

DT< I, J,MN) = (DTt I, J,MN)*DT< I , 7, MN) +6 . *DTt I,J, 1)>/D 

283 CONTINUE 

C POINT TOTAL UPDATE 

DTI 1 , 7, MN) =D 

284 CONTINUE 
299 CONTINUE 

GO TO 1© 

40© NOFTR= NOFTR+ 1 

IF I N0FTR. LT . NOFT) GO TO 99 
C DETERMINATION OF USEABLE DATA RANGES 

DO 1099 1=1,4 
J=2 

1013 J=J+1 

IFCBTt 1,7, J> .LT.299. >GO TO 1013 
IERt I,6)=J 
J= 125 
1017 J=J-1 

IFtDTt 1,7, J) .LT,299.)G0 TO 1017 
1099 IERt 1+4,6) =J 
C MEAN ADJUST 

DO 410 1=1,4 
Ll= IERt 1,6) 

Kl= IERt 1+4,6) 

DO 41© NS=L1,K1 
SUM= 0 . 

DO 405 J=i,6 

405 SUH= SUM+DTt I , J , NE) 

AVG= SUM/6. 

A=0. 

L= 1 
IC = 3 

DO 406 J=l,6 
B=DTt I, J,NE)-AVG 
IFtB.LT.0)B=-B 
IF IB. LT. A) GO TO 406 
A=B 
L= J 

406 CONTINUE 

SUM= SUK-DTC I , L , NE) 

AVG=SUM/o . 

IFCL.EO. K)K= 1 
A=0. 

DO 407 J= 1 , 6 

IFt J.EGLDGO TO 407 

B=DTt I , J, NE) - AVG 

IFt B. LT.0 . ) B=— B 

IFtB.LT. A)GO TO 407 

A=B 

K= J 

407 CONTINUE 
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SUH= SUM-DTC I , K, RE) 

410 DTC I,8,NE)=SUM/4. 

C SMOOTHING OF SPARSE DATA OVER USEABLE DATA RANGES 

DO 1199 1=1,4 
J=IER( I,6)+I 
K= IERC 1+4, 6) -1 
DO 1199 L=J,K 

IFC DTC 1 , 7 T L) . GT. 299 . ) GO TO 1199 
N=L 

1113 N=N- 1 

IFC DTC I , 7 , N) . LT. 299 . ) GO TO 1113 
M=L 

1117 M=M+1 

IF C DTC 1 , 7 , M) . LT. 299 - ) GO TO 1117 
RAD=DTC I ,7,L>/300. 

SK= 1 , -RAD 

IF( DTC 1 , 7 , L> . LT. 10 . ) DTI 1 , 8, L) =FLOAT( L-2) 

A=( DTC 1 , 3, L) -DTC I,8,N) )/( DTC 1, 8, M) -DTC I,8 ( N)) 

DO 1121 Ll= 1 , 6 

BG=DTC I, LI, M) -DTC I ,L1 ,N) 

BG=A*BG+DTC I,L1,N) 

1121 DTC I , L 1 , L) = SM*BG+RAD*DTC I , L 1 , L> 

1199 CONTINUE 

IFC NPRNT. EQ,. 0) GO TO 500 
DO 499 NE=2r 63 
MN=NE-2 

NRITEC 12,401)MN 

40 1 FORMATC 9 X, '’BAND 1 11 , 6X, "BAND 2 " , 6X, "BAND 3 " , 6X, "BAND 4 B , 6X, "MEAN " , T2) 
DO 404 J= 1 . 8 
A=DTC 1 , J,'NE) 

B=DT(2, J,NE> 

G=DTC3, J,NE) 

B=DTC4, J,NE) 

NRITEC 12,402) J, A, B,C,D 
4©2 FORMATC 2X, "DET" , IX, 1 1 , 4( 3X, F9 . 2) > 

404 CONTINUE 

499 CONTINUE 

G MID-RANGE TABLE BUILDER 

500 DO 498 1=1,4 
DO 498 NE= 1 , 6 
L= IERC 1,6) 

IC= IERC 1+4,6) 

SM=DTC I,NE,L) 

K=DTC I,NE,K) 

J=SM 

A=FLOATCJ) 

IFC SH. HE. A) J“J+1 
LIMSC I,NE, 1)=J 
LIHSC I,NE,2)=K+2 
DO 498 H=J,K 
A= FLO AT CM) 

K=M+1 

431 IFC DTC I , NE,L) . LE. A)GO TO 435 
L=L-1 
GO TO 431 
435 L2=L+ 1 

IFC DTC I , NE, L2) . GT. A) GO TO 44 1 

L-L+l 

GO TO 431 

441 A2= DTC I , NE , L2) -DTC I , NE , L) 

B=A-DTC I , NE, L) 

B2=DTC I,S,L2)— DTC I,8,L) 

LTCI , RE , N) =DTC 1 , 8, L> +B*B2/A2+ . 5 
498 L=L + t 
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C BE LOW- RANGE TABLE BUILDER 

DO 415 1=1,4 
DO 415 NE=1,6 
K=LIMS! I,NE, 1) 

L=K 

DO 415 J=1,K 
M=LT! I , RE , L+ 1 ) — 1 
IFCM. LT.0)M=0 
LT! I,NE,L)=M 
415 L=L-1 

C ABOVE-RANGE TABLE BUILDER 

DO 425 1=1,4 
DO 425 NE= 1 , 6 
K=LIMS! r,NE,2) 

DO 425 L=K, 128 
ir=LT( I , NE , L— 1 ) + 1 
IFCH.GT. 127) M= 127 
425 LT(I,NE,L)=M 
C TABLE SMOOTHER 

DO 525 1=1,4 
DO 525 J= 1 , 6 
BO 525 K=2, 126 
ri=rc+i 

L=LTC I , J,H) — 1 

IFCL.LE.LT! I,J,K))GO TO 525 
N=K-1 

IF!LT! I, J,R) .LT.LTC I, J,K> ) GO TO 521 
LT! I , J , K) =LT( I , J , K) + 1 
GO TO 525 
521 L=K+2 

IF(LT( I, J,L) .GT.LTC I, J,M))GO TO 525 
LT( I,J,ID=LT! I, J,M)-1 
525 CONTINUE 

DO 541 1=1,4 
DO 541 J= 1 , 6 
541 LT( I , J, 1) =0 

IF( NPRNT. Ed. 0) GO TO 550 
DO 549 K= 1,128 
MN=K-1 

WRITE! 12 , 60 1 ) MN 

601 FORMAT! 6X, "CHNL 1 CHNL 2 GHNL 3 CHNL 4 

DO 549 J= 1 , 6 
L=LT( 1 , J, K) 

L2=LTC2,J,K) 

M=LT(3, J,K) 

w- T t V( 4t T TT1 

WRITE! 12 ’ 603) J , L, L2, II, N 
603 FORMAT! 2X, “ DET " , I 1 , 3X, 4! 13 , 620 ) 

549 CONTINUE 
C 

C DATA CORRECTION 

C 

550 CONTINUE 

WRITE BINARY! 3) LT 
CALL FCLOSE! 3) 

CALL MTDIO! 1 , 10000K, IEM, IS, IERR) 

666 STOP 

701 STOP STRNF2 TAPE TABFIL 
END 


INPUT* 



i . 


i -1 


Opf^ L . PA -G£ is 


P0 °* 

\J "t v r V» T V' 'v <v* *V »T»-'T^ ip *ys«-J» n' r {»»V' , T' *V»* >* T i»*y» , > > I p" <r» •v^n<*l'*T* , T'*'l >T i' , *'I' , i'* , '«'‘ , l' , *^' , I' n'H' »J* *f* *%» »,■ 

G * 

C SUMCAT it 

C :;■ 

G THIS ROUTINE SUMS CLASSIFIED CATEGORIES INTO 

C GENERAL CATEGORIES FOR OUTPUT * 

C * 

G SEQUENCE: SUMCAT IRFIL * 

G «f* 

■ +}* *f- »f* +1* ■..)■» *!> »1* *1* «0> ■> » N •> *t* *J> */»»!' U* »j* iJ* \U*l* *J* iJj *> < i> *> <> ‘J* 

*fv »?» ^ ff> «N ^ ''XT' 'f'n> fi* ^><1^ *t ♦ * 1 » 'k'TT* n» W * «r» *i» H' 4* ^ ^ .Jh /(* ^. <y* <N *■!- 

C * 

C CREATED AT GEORGIA TECH EES * 

C «J* 

C PROGRAMMER: NICKOLAS FAUST * 

G * 

DIMENSION COUNT! 60) , IEQ! 60 * 20) , NUM! 20) , I CAR! 10,20) 

D I PENSION JFIL! 34) ,LSUM!2G> , PER! 60) ,KFIL(34) , ISNSC2) 

GALL OPEN! I , "COM. CM" , 1 , IE) 

CALL COMARG! 1, JFIL, ISNS, IE) 

CALL COMARG! 1, JFIL, IS¥S, IE) 

CALL OPEN! 3 , JF IL , " B " ) 

CALL COMARG! 1,KFIL, ISWS, IE) 

CALL FGPEN ! 2 , KF IL , "B " ) 

READ! 3) I SUM 
DO ! L= 1 , ISUM) 

: READ !3, 101) ! ICAR! J,L) , J=l, 10) 

101 : FORMAT ( 10A2) 

READ! 3) NUM! L) 

NUM2=NUH( L) 

READ! 3) ! IEQC J,L) , J= 1, NUM2) 

. .FIN 

READ! 3) NWQ. 

DO CIC=1,NWQ) 

STOT= 0 . 

DO ! EG* 1,60) COUNT! KG) = © . 

READ! 2, 130) NA1 

130 : FORMAT! 14) 

¥RITE( 12, 131) NA1 

131 : FORMAT! 2S, 11 WQMU * ", 15///) 

DO CKZ=1,60) 

: READ! 2,70) COUNT! KZ) t PER! KZ) 

70 : : FORMAT! 2K, F 12 , 2 , 2X, F 12 . 2) 

: ¥RITE( 10) COUNT! KZ) , PER! ICZ) 

: . .FIN 

DO ! N= 1 , ISUM) 

NUM2-NUMCN) 

82=0 
S1=0 

DO (L= l t NXJM2) 

IQ=IEQCL,N) 

S1=S1+C0UNT! IQ) 

S2= S2+PERC IQ) 

.FIN 

WRITE! 12, 1 ©2.) N , < I CAR! J , N) , J-l, 10) ,S1,S2 
ST0T= STOT+S 1 

102 : : FORMAT! 2X, 12, 2X, 10A2,2X, F10. 1 , " ACRES “ , 2H, F6 , 3 , « PERCER1 

FIN 

WRITS! 12, 103)STOT 

103 : FORMAT! /6X, "TOTAL ", 10X.F12. 1, » ACRES") 

WRITE! 12,200) 

20© s FORMAT! ////////////) 

.FIN 
STOP 
END 


1 
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C sj; 

C SYMINV (SUBROUTINE) * 

G sjc 

C INVERTS A SYMMETRIC MATRIX * 




C * 
C CREATED AT NASA JSC * 
C * 


SUBROUTINE SYMINVC A, AI , DET, N) 

DIMENSION A(N,N) , AHN»N) , Q.CC 12) ,D( 12) ,DIC 12) 
INTEGER R 
DET= 1 . 0 
R=0 

10 R=R+ 1 

IF(R.GT.N)GO TO 19 
aC(R)=1.0 
1=0 

20 1 = 1+1 

IFC I.GT.R-DGO TO 29 

QCC I) =0.0 

J=0 

30 J=J+1 

IFC J.GT.R-DGO TO 20 
QGC I) = OCC I)-AIC I, J)*A< J,R) 

GO TO 30 


29 CONTINUE 

DCR) =A( R,R) 

IF(DCR) ,LE.0.0)GO TO 52 
IC=0 

40 IC=K+1 

IFCK.GT.R-DGO TO 49 
DC R> = DC R>+ AC R, ID *Q.C(K) 

GO TO 40 
49 CONTINUE 

DET=DET*DC R) 

IFC C DC R) /AC R, R) ) . LT. 1 . E-8) GO TO 52 

51 DIC R) = l. 0/DC R> 

GO TO 60 

52 DICR) =0. 0 

IVRITEC 12, 10 10) R, DCR) ,DET 

1010 FORHATC9H AT STEF ,I1,4H D =,E10.5,6H DET = ,E10.5) 

6© I=© 

70 CONTINUE 

1 = 1+1 

IFC I ,GT. R) GO TO 100 
AICR, I)=0. 0 
AI C I , R) = 0 . 0 
J=0 

80 J=J+1 

IFC J , GT. R) GO TO 70 
€fl)=QCC DittD ICR) 

AIC I, J>=AICI,J)+QD*QCC J) 

GO TO 80 
100 CONTINUE 

GO TO 10 
19 CONTINUE 

RETURN 
END 
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^ n» 'r» 'i> *T fc - r f ■* *t* > ' t* *t> *t» -t* U>«T» h* n» , T'n»’v *r» n> r x* n’ n* rv» n» •{' •>» n* ^ n*T»n» «f* • H* 'ir'-f* <^v 

G 

THRDST (SUBROUTINE) 


G 

G 

G 

G 

G 

G 

C 

G 

G 

C 

G 

C 

G 

G 

G 

G 


PURPOSE 


UPDATE THE MEM AND VARIANCE OF THE THRESHOLD STATISTICS 
DESCRIPTION OF PARAMETERS 


INPUT 

GALLING SEQUENCE 

XI, VI, N I - CURRENT MEAN, VARIANCE, AND WEIGHT 
X2,N2 - VALUE AND WIGHT OF NEW SAMPLES 


OUTPUT 

CALLING SEQUENCE 

XI, VI - UPDATED VALUES 


C 

C CREATED AT NAS A/ JSC CASTEP) 

C 

SUBROUTINE THRDSTC XI, VI, N1,X2,N2) 
REAL N3 

IF(N1.GT,0.) GO TO 1© 

X1=X2 

V1=0. 

RETURN 

10 X- ( N1+FLQATCN2) ) 

R1=(N1)/X 

R2=FL0ATCN2)/X 

X=R1*X1+R2*X2 

V1-R1*( V1+X1**2> +R2*( X2**2) -X**2 

xi-x 

RETURN 

END 
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G 

C T0P02 

C 

G READS CYBER CONVERTED TOPO TAPES AND 

C DISPLAYS THEM TO THE COMTAL. 

C 

G SEQUENCE: T0P02 KTU:F 

C 

*J' l P ^ *i* «-p "J+ M J *J/ ^V»>V* * ■■» ■,!» -J> *J » 'J-* »A* '•l* «1* 4* » U «. I .. - U *!*■ * U *. U » L, v t, *.}+ *J> k L, « t-. ^ . U •■ I ■• »>* tla si- O# *1- ^ * f* *J> »i> *J a * U « 

U , v* f 4'*V *N *v» n> •v* »T* *l' 1 »>* *l> *T* *T» T* rf> *T» •!> «V> »{» rf* fp »}» *N *t* *i ■ r - 


:H 


C 

C CREATED AT GEORGIA TECH EES 

G 

C PROGRAMMER: MICHAEL D. FURMAN 

G 

INTEGER IMP( 1810) , IK0KKX512) , IHAG(256) , ITAPE( 10) , ISFC2) 
CALL OPEN! 1 , 11 COM. CM" , 1 , 1 ERR) 

CALL COMARGC 1, I TAPE, IST7, IEPJI) 

CALL COMARGC 1 , 1 TAPE, ISW, IEHR) 


Hi 


( 1=YES,0=N0) 
IFSCAN- 4 


ACCEPT "FAST SCAN? 

NHEN ( IFSCAN.EQ. 1) 

ELSE IFSCAN=1 

ACCEPT "SKIP RECORDS? " , ISEP 
ACCEPT "MAM ACCEPTED PEAK VALUE 
P2= I PEAK/25 5 . 0 


IFSCAN 


(AFTER SUBTRACT TON ) 


IPEAK 


ACCEPT "ASSIGNED PEAK SHADE VALUE? " , ISHADE 
ACCEPT "SUBTRACT NHAT VALUE FROM ALL DATA? " , ISUB 
ACCEPT "BEGIN WITH ELEMENT " , I EL 
IEL= IEL+7 


IEND= XEL+5 1 1 

IF C IFSCAN.EQ. 4) IEND= 180O 
IF C IEND.GT. 1800) IEND=180C) 
CALL MTOPDC 3 , ITAPE , 0 , IE) 
LINE= ISKP 


DO (1=1, ISKP) CALL MTDIO(3,0, IMP, IS, IE, IC) 

DO CHS* 1,256) 

: TYPE "LINE NUMBER ".LINE 

: LINE=LINE+ IFSCAN 

: DO (1=1, IFSCAN) CALL MTDIOC3.0, IMP, IS, IE, IC) 

: K= 1 

: DO ( I=IEL, IEND, IFSCAN) 

: : IWOKKC K) = IMP( I ) 

: : K=IC+1 

: : . .FIN 

: K2=K-1 

: SUBTRACT-CONSTANT- AND-CUT-OFF-PEAKS 

: CALL PAC8( IWORK, I MAG, 5 12) 

: CALL RIMWRITE(0,G, IMAG.256) 

: IP (M5.EQ.2o6) ACCEPT "TYPE A *1* TO CONTINUE MANS 

: IF ( IANS. EQ. 1) M5= 1 

: . .FIN 
STOP 

TO SUBTRACT-CONSTANT- AND-CUT-OFF-PEAKS 
: DO (M6=1,K2) 

: : WHEN ( I TvORK( M6 ) . LT. I SUB) I NORK( M6 ) = 0 

: : ELSE I WORK( M6 > = I WORK( M6 ) - 1 SUB 

: : IF ( I¥0KK(M6) .GT. IPEAK) IW0RKCM6) = IPEAK 

: : I ¥0RK( M6 ) = IWOKKX M6 ) /P2 

: :..FIN 


: . .FIN 
END 
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'J* ‘l* *1* *i**l'*J' *4» M*»V \U «j> *1> *1* »^» ^ t]> »V- tji »J» *•> *|* >i> »l»«l»«Jj>ii*J**Ljj »X» »U iL* ■J* *5^ *L»kA* vl' 

*V* 'T' *p T' "T* *v «T* r \ tl *4^ ffv if* #|> ft rf* r]» rf* /[ w 7« f]ini q**o 

C si* 

*T* 

c 

G 
G 
C 
G 
C 
G 
G 


TRAIN3 

THIS PROGRAM IS DESIGNED FOR TRAINING FIELD SELECTION 
AND STATISTICS GENERATION FOR THREE IMAGES 

SEQUENCE: TRAIN3 MTU:F 


* 

* 
Hi 
Hi 
Hi 
Hi 

C Hi 

C DESIGNED AT GEORGIA TECH EES * 

'i* 
Hi 
Hi 
H: 
H: 


PROGRAMMERS: 


NICKOLAS L. FAUST 
ROBERT A. MADDOX 
MICHAEL D. FURMAN 


G 
G 
G 
G 
C 

DIMENSION V( 4, 4) ,XM(4) , IVC 101) , IVX( 101) , IVY! 101) 

DIMENSION D( 4,2) ,HXC4) , IDEVC30) 

INTEGER FIELDC 17) , ISWSC2) 

COMMON /DTRANS/ IDATC256) , IXD(512,4) , IDUMCU12) 

COMMON /HISTX/ IG0UNTC4, 100) _ 

COMMON /BUFFER/ IBUFC 1700) ORIGINAL PAGE i£> 

£ OF POOR QUALITY 


G 

c 


OPEN-FILES 
I NIT- VARIABLES 
FOREVER 

ASK-QUESTIONS 
CONDITIONAL 
C IFUN.EQ. 1) 

( IFUN.EQ. 2) 

C IFUN.EQ. 3) 

C IFUN.EQ. 4) 

C IFUN.EQ. 5) 

( IFUN.EQ. 6) 

C IFUN.EQ. 7) 

C IFUN.EQ. 3) 

( IFUN.EQ. 9) 

C IFUN.EQ. 10) 

C IFUN.EQ. 11) 

C IFUN.EQ. 12) 

C IFUN.EQ. 98) 

( IFUN.EQ. 99) 

C OTHERWISE) 

. TYPE "PLEASE 
: . .FIN 
.FIN 
.FIN 


READ- IMAGE 
ALARM-SCREEN 

MAGN I FY-SUBSET-OF- SCREEN 
TAKE- TRAIN ING-SAMPLE 
TAKE-GCP 

D IGITIZE-BOTJNDARY 
GALL FSWAPC "SGRIPTO. SV") 
GALL FSWAPC "INK. SV") 

GALL FS WAP ( " WFMBOB . S V " ) 
CALL FSWxAPC "GHCOLR.SV") 
CALL FS WAP ("FI LTER2 . S V " ) 
GALL FSWAPC "CLYDES. SV 11 ) 
BLANK-GRAPHICS 
CLOSE-FILES 


ENTER ONE OF THE NUMBERS GIVEN" 


96 


TO OPEN-FILES 

: CALL OPENC 1 , "COM. CM" , 1 , IFER) 

: CALL COMARGC 1 , FIELD , ISW3 . IFER) 

: GALL COMARGC 1, FIELD, ISWS, IFER) 

: IF C FIELDC 1) .NE. "MT") 

: : TYPE "PARAMETER ERROR TRAIK3 MTU:F" 

: : TYPE "TAPE NAME PROBABLY NOT ENTERED" 

: : STOP PARAMETER ERROR IN TRAINS KTU:F ERROR HALT 

: : . .FIN 

: GALL MTOPDC 2 , FIELD , 0 » IFER) 

: CALL FOPENC 3, "VERTEX", "B") 

: GALL FOPENC 4 , "HISTS " , »B " ) 

: CALL FOPENC 5 , "SAVS IG" , "B") 

: TYPE "INPUT TAPE NUMBER IN 612 FORMAT 11 

: REAIK 11,96) IT1 , IT2, IT3, ID1 , ID2 , ID3 

: FORMAT! 6 12) 

: WRITE BINARY! 3) IT1 , IT2, IT2, ID1, ID2, ID3 

: WRITE BINARY! 4) IT1 , IT2, IT3, ID1 , ID2, ID3 

: WRITE BINARY! 5) IT1 , IT2, IT3, 1D1 , ID2, ID3 

: CALL KTD I DC 2 , 0 , IBUF , IS , IFER) 

: CALL 3\1TD IOC 2,0, IBUF, IS, IFER) 

: . . Ft t-i 125 


c 

c 


4 - 


C 

C 


C 

G 


TO IN IT- VARIABLES 
i MODE=0 
: MODE 1=2 

: MAG= 1 

: 1GH= 1 

: ND=4 

: . . F IN 


TO ASK- QUEST I ONS 


TYPE 

n 

11 




TYPE 

TYPE 

11 THE 

U If 

FOLLOWING FUNCTIONS ARE 

TYPE 

It 

1 

— 

READ IMAGE 

9 - 

TYPE 

II 

2 

— 

ALARM 

10 - 

TYPE 

11 

3 

— 

MAGNIFY 

11 - 

TYPE 

11 

4 

— 

TRAINING 

12 - 

TYPE 

11 

5 

— 

GCP " 


TYPE 

II 

6 

— 

DIGITIZE" 


TYPE 

II 

7 

— 

USE SCRIPTO 

98 - 

TYPE 

TYPE 

11 

II 

8 

II 

— 

USE INK 

99 - 

ACCEPT 

FIN 

"FUNCTION NUMBER ? 

" , IFUN 


SELECT BY NUMBER" 


USE 

USE 

USE 

USE 


WM" 
GUCOLR 11 
FILTER2" 
CLYDES " 


BLANK GRAPHICS" 
STOP" 


TO TAKE- TRAINING-SAMPLE 
: K=1 

: I GLCSE= 1 

: PAUSE POSITION CURSOR FOR TRAINING FIELD THEN RETURN 

: CALL GETXYC IXS, IYS) 

: IXS=IAND( IXS, 37710 

: I YS= I AND (IYS, 377K) 

: IVXC 1 ) = IRS 

: IVYC 1) = IYS 

: IXS 1= IXS 

: IYS1= IYS 

. TYPE u CURSOR COORDS « f IXS , IYS 
: TYPE * POSITION CURSOR AT NEXT VERTEX" 

■ ACCEPT " ENTER 1 TO READ POSITION, 2 TO CLOSE : MCLOSE 
: REPEAT UNTIL ( ICLOSE. NE. 1) 

: : CALL GETXYC IX, IY) 

: : CHECK"* XY 

: : CALL VECTOR( 1, IXS, IYS, IX, IY, I, 1) ; DRAW ON GRAPHICS 1 

: : TYPE n END PTS OF SEGMENT" , IXS, IYS, IX, I Y 

: : K=K+1 

: : IVXC K> = IX 

: : IVYCK) = IY 

: : IXS= IX 

: : IYS= IY 

: : ACCEPT " ENTER 1 TO READ POSITION , 2 TO CLOSE : MCLOSE 

: : . . FIN 

: CALL VECTORC 1, IXS, IYS, IXS1 , IYS1, 1,1) 

: IYMAX=0 

: IXMAX=0 

t IXFIIN= 100OO 
: IYIIIN= 10000 

: ivx(ic+i) = rvx( i) 

: IVYC K+l) = IVYC 1) 

: KP=K+1 

: DO (IR=1,KF) 

i : IVXC IK) = IVXC IR) /EL4G+ ICOL 

: : I VYC IR) = IVYC IR) /MAG+ IL 

: : IFC IVYC IR) . GT. IYPLAX) IYMAX= IVYC IR) 

: : IFC IVYC IR) , LT. IYMIN) IYMIN= IVYC IR) 

* : IFC IVXC IR) . GT. IXMAX) IXMAX= IVXC IR) 

: : IFC IVXC IR) .LT. IXMIN) IXMIN= IVXC IR) 

: :..FIN 
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25 


107 


G 

G 


C 

c 






PAGE IS 
0P QUALITY 


ISTART= IYMIN 
IEND= IYMAX 

TYPE " EOS LIMITS 11 , ISTART, IEND , IXMIN, IXMAX 
GALL RDATAC ISTART, IEND.K, I VS, IVY, IXD,NP) 

GALL HIST3B(NP,3) ;POT HISTOGRAM ON GRAPHICS 3 

ACCEPT " INPUT 1 FOR STATS : ", ISTAT 

CONTINUE 

IF( ISTAT, EQ. 1) 

: DO ( Jl= 1 , 4) 

: : 5iMCJl)=0.0 

i : DO ( J2= 1 , 4) 

r : : VCJ1,J2)=0.O 

: : : . .FIN 

: : . . F IN 

: 1=0 

: DO (LK= 1 , NP) 

: : DO (LJ= 1,4) 

i : : MXC LJ> = IXDC LIC, LJ) 

: : i . .FIN 

: : TYPE LK 

: : GALL SEQCOVC MX, 4 , V, XM, I , D) 

: i . .FIN 

: GALL MATPRTC XM, ND , ND , 1 , 5HHEANS) 

: GALL MATPRTC V, ND , ND , ND , 3HC0V) 

,* TYPE 11 STATS CALCULATED FOR ".NP," POINTS" 

: ACCEPT " DO U NISH TO SAVE THIS S I GNATURE? I = YES : « , IS IG 

; IFC ISIG.EGt. 1) 

: : TYPE " INPUT 8 CHARACTER NAME FOR SIGNATURE" 

: : EEADC 1 1 , 107) NA1 , NA2 , NA3 , NA4 

t F0RMATC4A2) 

*. : WRITE BINARYC5)NA1,NA2,NA3,NA4 

: : WRITE BINARYC 5) NP, XM, V 

: : WRITE BINARYC 3? NA1 , NA2 , NA3 , RA4 

t : WRITE BINARYC 3) SC, IYMIN, I YSIAX 

: : DO (Kl= 1 ,K) WRITE BINARYC 3) K1 , IVXC Kl) , IVYC Kl) 

: : WRITE BINARYC 4) NA1 ,NA2,NA3,NA4 

: : WRITE BINARYC 4) I COUNT 

: ..FIN 
: . .FIN 

ACCEPT 11 HAS SIGNATURE BEEN TREATED AS U WISH? 1= YES : ", I SAT 
IFC ISAT.NE. 1) 

: ISTAT= 1 

; GO TO 25 
: . .FIN 


TO CLOSE-FILES 
: CALL FCLOSEC 3) 

: CALL FCLOSEC 5) 

: CALL FCLOSEC 4) 

: STOP TRAINS MTU:F NORMAL EXIT 

: . .FIN 


TO CHECK-MAG-FACTOR 

: WHILE CKAG.NE. 1.AND.MAG.KE.2.AND.MAG.NE.4. AND.KAG.NE.8) 

: : TYPE "INCORRECT MAG FACTOR" 

: : ACCEPT "ENTER MAG FACTOR <1,2,4, OR 8) " , MAG 

: :..FIN 

i . .FIN 
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■ =» 


V 


G 

C 

TO READ- IMAGE 

: ILT= IL 

: ICQLT=ICOL 

: TYPE "INPUT TOP LEFT COORDS FOR IMAGE C COL, LINE) 1 ' 

i ACCEPT " ENTER -1,-1 TO KEEP OLD COORDS " , ICOL, IL 
: IF (IL.EQ.-l) 

i : IL= ILT 

: : ICOL= ICOLT 

: .‘..FIN 

: ILS= IL 

: IES= ICOL 

: ACCEPT "MAGNIFICATION FACTOR (1,2,4, OR 3) ? " , MAG 

: C HECK- HAG- FACTOR 

: CALL RDIMGC IL, ICOL, I CH, MAG) 

: . .FIN 
C 

c 

TO ALARM-SCREEN 

: ACCEPT "ALARM TO WHICH OVERLAY (0-3) " , IOV 

: CALL ALRM2C IOV) 

: . .FIN 
C 
C 

TO P1AGN IF Y-SUESET-OF-SCREEN 

: PAUSE POSITION CURSOR AT UPPER LEFT AND HIT RETURN 

: CALL GETXYt IX, IY) 

i CHECK- XY 

: ICOL= IX/MAG+ IES 

: IL= I Y/MAG+ ILS 

: TYPE "SCREEN COORDS", IX, IY 

: TYPE " TAPE COORDS", ICOL, IL 

: ACCEPT "MAG FACTOR (1,2,4, OR 8) ? " , MAG 

: CHECK-MAG-FACTOR 

: CALL RDIMGC IL , ICOL, ICH, NAG) 

: * . FIN 
C 
C 

TO TAKE-GCP 

: ACCEPT "INPUT 1 TO READ CURSOR POSITION ",IG 

i REPEAT UNTILC IG.NE. I) 

: : CALL GETXY( IX, IY) 

: : CHECK- XY 

: .* TYPE "SCREEN COORDS ",IX, IY 

: : ICOLKT= IX/MAG+ICOL 

: : ILET= IY/MAG+ IL 

: .* TYPE " TAPE COORDS", ICOLKT, ILKT 

: : ACCEPT "INPUT 1 TO READ CURSOR POSITION ",IG 

: : . .FIN 

.‘..FIN 
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C 

C 

TO DIGITIZE-BOUNDARY 

: ACCEPT "DRAW BOUNDARY ON WHICH OVERLAY (0-3) ? 'MOV 

: PAUSE POSITION CURSOR AT START AND HIT RETURN 

: CALL GETXYC IX, IY) 

: CHECK-XY 

; type "SCREEN COORDS ",IX, I Y 
: ICOLS= IX/NAG+ I COL 

: ILS2= IY/MAG+ IL 

: TYPE " TAPE COORDS M COLS, XLS2 

: 1X1= IX 

: I Yl= IY 

: IXS=IX 

: IYS= I Y 

: ACCEPT "ENTER: 1 TO READ POSITION, 2 TO CLOSE, 3 TO STOP “ , ICLOSL 

: REPEAT UNTIL! ICLOSE.NE. 1) 

: : CALL GETXYC IX, IY) 

: : CHECK-XY 

: : CALL VECTOR! IOV, 1X3 , IYS , IX, I Y, 1 , 1) 

: : TYPE “LAST VERTEX AT” 

i ; TYPE " SCREEN COORDS ",IX, I Y 
: : ICOLS= IX/HAG+ I COL 

: : ILS2= IY/ILAG+IL 

; : TYPE " TAPE COORDS", I COLS, ILS2 

: : IXS=IX 

: : IYS= IY 

: : ACCEPT "ENTER: 1 TO READ POSITION, 2 TO CLOSE, 3 TO STOP ",ICL 

: : . .PIN 

: IF ( ICLOSE.Ea.2) CALL VECTOR! IOV, IXS, IYS, 1X1 , IY1 , 1 , 1) 

: . .FIN 
C 

c 

TO CHECK- XY 
: IX= I AND! IX, 37710 

: IY= I AND! IY, 377K) 

: . .FIN 
C 
C 

TO BLANK-GRAPHICS 

: ACCEPT "WHICH OVERLAY TO BLANK ? C0-33 " , IGRN 

: DO (1=1,16) IDUMI)=0 

: DO (1=0,255) CALL GWR! IGRN, I , IDUM, 16) 

: . .FIN 
C 
C 

STOP 

END 
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G ^ 
G UHP AG 1 (SUBROUTINE) * 
G ^ 
C UNPACKS LANDSAT GODDARD FORMAT * 

G if: 


c ^ 

C "CREATED AT NASA/ JSC ( ASTEP) $ 

C * 

C********************#***********************************:,':**;f::!c;fc:{i:;:i|::f:;i:if:il:;|:ilc;}: 
SUBROUTINE UNPAC1 
INTEGER BUF 

GOMKON/BUFFER/BUFC 3300) 

COMMON/COL/ JS , JF 

CALL MDIO( 4,0, BUF , ISTAT, TERR, NW) 

IF( NW. GT. 1750) miTEC 10 , 1) 

WRITEC 10,2)NW 

2 FORMAT! 2X, "# OF WORDS READ = ",I5> 

1 FORMAT ( 2X, "ERROR IN UNPAC 1 ° ) 

IM= 1630 


JC=0 

DO 20 1=1, NW 
NM= 1+IM 

20 BUFC NM) =BUF( I) 

DO 3 J= JS, JF 
I=J-1 

Ll=-16+C I+C I/2)*6)*8 

JC=JC+1 

DO 4 1=1,4 

IL=L 1+1*16 

11= IL/16 

12= IABSC IL— I 1*16) 

11 = 11+1 

IF( I2.EQ., 8)G0 TO 10 
11= ISHFTC BUF ( 1 1+ IM) , “8) 
GO TO 5 

10 CONTINUE 

II=IAND(BUF( II+IM) ,377K) 
5 CONTINUE 

KK= I+( JC-1) *4 
BUF( KK) =255—11 
4 CONTINUE 

3 CONTINUE 

RETURN 
END 


/ 
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V 


*> 


ofp^t L page js 
* POOR QUALITY 


GXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
G X 
C UNFAC4 (SUBROUTINE) * 
G X 
C UNPACKS PART OF LINE INSIDE POLYGON X 
G X 

cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
c X 
G CREATED AT GEORGIA TECH EES * 
G *}£ 
G PROGRAMMER: NIGKOLAS L. FAUST * 
G * 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


SUBROUTINE UNPAC4C IBUF , JBUF , JL , JR, N , M) 
D I MENS ION JBUF C 4 , N) , IBUF ( II) 
rF< JR. GT. 810) JR=810 
DO CK-1,4) 

: L= C ( JL*~ 1 ) /2) #4+\K 

: LAST= L+ ( JR-JL+ 1 ) *2 

: Jl=l 


: DO (I=L,LAST,4) 

: : JBUFCK, Jl) = ISHFTf IBUFC I) ,-8) 

: : JBUFCK, J1+1) = IAND( IBUFC I) ,377K) 

: : J 1- Jl+2 

: : . . FIN 

: . .FEN 

RETURN 

END 
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C * 

UPPLT ( SUBROUTINE) 


PURPOSE 

GENERATES PRIORITY LIST GIVER WEIGHTS 
DESCRIPTION OF PARAMETERS 
INPUT 

CALLING SEQUENCE 

NVG - WEIGHTS, I-TH VALUE IS NUMBER OF POINTS IN I-TH 
CLUSTER 

NVM - NUMBER OF ENTRIES IN NVG 
OUTPUT 

CALLING SEQUENCE 

PLIST - PRIORITY LIST, J-TH VALUE IS INDEX OF THE J-TH 
LARGEST ENTRY IN NVG 


s): 
* 
* 
X 
:Ji 
# 
X 
* 
X 
X 
5ft 
X 
* 
* 
X 
X 
X 
# 

C ^ 

C CREATED AT NAS A/ JSC CASTEP) * 

c X 

SUBROUTINE UPPLTC PLIST, NVG, NVM) 

REAL NVG,N 
INTEGER PLIST, TLIST 

DIMENSION FLISTCNVM) ,NVG(NVM) , TLISTC20) 

DO 10 1=1,20 
10 TLISTC I) =0 
DO 30 L= 1 , NVM 
N=— 1 

DO 20 1=1, NVM 
IFCTLISTC I) .EQ. I) GO TO 20 
IFCNVGt I) .LE.N) GO TO 20 
N=NVG( I) 

J=I 


20 CONTINUE 
TLISTC J) = l 
30 PLISTCL) = J 
RETURN 
END 
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G * 
C ZCOUNTY * 
C * 
G CALCULATES PIXEL COORDINATES FOR ANY COUNTY PARTIALLY * 
C OR WHOLLY CONTAINED IN A SCENE AND WRITES RESULTS TO DISK. * 
C * 
C SEQUENCE: NCOUNTY TAPE FILE * 
C «f* 

C # 
C CREATED AT GEORGIA TECH EES * 
C * 
C PROGRAMMER: NICKOLAS L. FAUST * 
C * 


*J* * 4 * 0 > •!« tit ^ 1 . tli tl. ii> <u %u 4 * vU U j tl< J/U. ^ il* > 1 > ^i/U. * 1 # tU >> ti/ 4> . 1 * U> «i> < 1 > U. ^ ■!> <t» 4j« ii> d. Ai> +±r •[. 


DIMENSION NGODEC 50) 
COMMON/DUM/ I AC 4 100) 

CALL OPENC I , "COM. CM" , 1 , IERR) 
CALL COMARGC 1, FIELD, ISWS, IER) 
CALL C0MARG< 1, FIELD, ISWS, IER) 
CALL MTGPD( 2 , F IELD , 0 , IE) 

CALL COMARGC 1, FIELD, ISWS, IER) 
CALL FOFENC3,FIELD,0, IE) 

I SET* 1 


ITEMPC 2, 300) , FIELDC 17) , ISWSC2) 


ORIGINAL PAGE IS 
OP POOR QUALITY 


I0U=3 


NCOD=0 

ACCEPT " A1 , A2 , A3 , B 1 , B2 , B3 ",A1, A2, A3,B1,B2,B3 
ACCEPT "SKIP RECORDS? n ,ISK 
IF CISK.GT.0) 

: DO £ I~ 1 , ISK) CALL MTD IOC 2,0, I A, IS, IE, IC) 

: . .FIN 
K=0 


IV=0 

IC=10 


200 

300 


UNTIL CIC.LT.S) 

: CALL MTDI0C2, 0, I A, IS, IE, IC) 

: TYPE " DATA 11 , C IAC KK) ,KK- 1, IC) 

: NCOUNTY=IAC 1) 

: YI= I AC 2) * 10O00 . + 1 AC 3) 

: XI* IAC 4) *10000. + IAC 5) 

: IFCXI.LE. .001. AND. YI.LE. .001) STOP 

: IV=IV+1 

: IFC NCOUNTY. NE.NCOD) 

: : ISET= 1 

: : JSET= 0 

: : KSET=0 

: : IFCK.NE.0) 

: : : NIV=IV-1 

: : : DO CJ=1,NIV) 

: : : : IFC ITEMPC 1, J) .GT. 0.AND. ITEMPC 1, J) . LT. 4000) JSET* 1 

: : : : IFC ITEMPC 2, J) .GT.0. AND. ITEMPC 2, J) .LT, 4000) KSET= 1 

: : : : IFC JSET.EQ. 1 . AND.KSET.EQ. 1) ISET=0 

: : : : JSET=0 

: : : : KSET=0 

: : *• : . .FIN 

: : : IFC ISET.EQ. 0) 

: : : : NIV*IV-1 

: : : : IV= 1 

: : : : WRITEC IOU)NCOD 

: : : : WRITEC I OU)NIV 

: : : : DO CKl= 1 ,NIV) WRITEC IOU, 200) ITEMPC 1,K1) , ITEMPC 2, Kl) 

: : : : FORHATC IX, 215) 

: : : : WRITEC I QU, 300) NCOD 

: : : : F0RMATC2X, " END OF PIXELS FOR COUNTY = M5> 

: : : :..FIN 

: : : . . FIN 


: : IV* 1 

i : NCOD= NCOUNTY 

: : K=1 

: : ISET= 1 

: : . .FIN 

: PY=A1+A2WXI+A3*YI 

: PX=Bl-t-B2*XI+B3*YI 

: ’ ITEMPC 1, IV)=PX 

: ITEMPC 2, IV) =PY 

: . .FIN 

STOP 

END 
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APPENDIX C 

MAIN PROGRAM FLOWCHARTS 
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ADDSIGM 



SSW ® 5 


INPUT - 1) SIGNATURE NAME FILE 
2) SIGNATURE FILE 
OUTPUT - NEW BINARY 

SIGNATURE FILE 
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CLUSTER 


INITIALIZATION 


INPUT 

CLUSTER 

BUILDING 

PARAMETERS 


/ READ IN 
f VERTICES 
OF POLYGON TOJ 
BE CLUSTERED/ 


READ 
LINE OF 
DATA 


INPUT - LANDSAT RAW DATA TAPE 
OUTPUT - CLUSTERED DATA SET 



LAST LINE? 


REWIND 

TAPE 
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(next page) 



dt - 


CLUSTER (Continued) 
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LINCLASS 



INPUT - LANDSAT RAW DATA TAPE 
OUTPUT - CLASSIFIED DATA TAPE 
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M60CL (Continued) 
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STRIFNF (Continued) 
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TRAIN3 



INPUT - LANDSAT RAW 
DATA TAPE 
OUTPUT - FILES FOR 

a) SIGNATURES 

b) VERTICES 

c) HISTOGRAMS 
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TRAIN3 (Continued) 
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V 


TRAIN3 D (Continued) 
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TRAIN 3 (Continued) 


2 


READ 

CURSOR 

POSITION 


* 


COMPUTE 

SCREEN AND TAPE 
COORDINATES 


TYPE OUT 
COORDINATES 









■ ** 
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TRAIN3 (Continued) 
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TRAIN 3 (Continued) 
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TRAIN 3 (Continued) 






TRAIN3 (Continued) 
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